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On the Problem of Transformation Bieaetien 
of the Spin Pseudo-Vector 


er 


Summary 


Polarization 4- vector of Dirac’s particles is introduced by the covariant 
method, and its transformation properties are investigated. The results ob- 
tained are used to discuss the phenomena related to parity non-conservation. 


$ 1. Introduction 
The questions of the spin properties of the particles have become recently 
of special interest due to the discovery of parity non-conservation!). We 
consider the case of free fermions with positive energy. We use the con- 
servation of the total angular momentum ?)’) as the basis for determining the 
” 


— {Mine Styria} —(. 


is the “orbital” angular momentum tensor, and aor 


Dura 
is the “internal” angular momentum tensor. salt 
T he | öy Cyt ) 


is the canonical energy tensor, x, = 0, 0, isthe velocity matrix, while the 
quantity 


ji Ü 01 On 03 0, 03 04 when all suffixes are different 
| 0 when at least two suffixes coincide ie = 


is the antisymmetrical tensor of the third order determining the spin com- 
ponents. The pseudo-vector spin matrix o, is defined by the relation 


') T.D. Lee, u. €. N. Yang, Physic. Rev. 104, 154 (1956); 105, 1671 (1957). 

*) D. Iwanenko u. A. Sokolow, „Klassische Feldtheorie‘“, Akad. Verlag, Berlin 
1953. 

3) A.A. Sokolov, “Introduction to the Quantum Electrodynamics” Fizmatgiz, 
M.-L. (Russ.) 1958, see also A. Sokolow, ‘Quantenelectrodynamik’’, Akad. Verlag, 
Berlin 1957. 
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In these equations 9, 9,, 03 are the well-known Dirac matrices 


0 010 0 ‘0 0 0 
vo0oıp 0 0 
0010 00:0 
1 0 0 0 0010 
0-10 0 0001 
000-1, 0100 
00 0 —1 
i’ being the second imaginary unit (see 3)), i’ = 72 — — 1 while i+ = —i, i’+ =’, 
et $ 2. The polarization vector 
er The longitudinal polarisation of Dirac particles may be described by the 
operator S = V + ei This is the integral of motion, its eigenvalues s 
iV_p2 \ 


being equal to the doubled projection of the spin pseudo-vector on the direction 


of the momentum k*). The transverse and “time” components not being the 
integrals of motion only their mean values may be given 
= x, {u = 1, 2, 3, (2) 
However, due to the volume element under the integral sign such a definition 
of polarization will have the other covariance properties than those of the 
quantity y* o, y. Therefore we take for the polarization vector 
Ko, = (3) 
This quantity transforms according to the same law as the quantity y* o, y. 
Here K = | k* + kg and k, characterize respectively the energy and rest 
mass of the particle. In case of the free particle with positive energy the y- 


function has the form?) er 
YC, 0, (é)exp(—ie Kt 4 i kr). (4) 


f (—K) Sin ©, e 
4) In ref.5) second quantization of the free-particle Dirac equation was carried 
out taking account of the subsidition Sy = sy. This made it possible to take account 
of the spinor properties of the particles. The method developed in ref. ®) was used in °). 
5) A. Sokolow, Journ. of Phys. USSR 9, 363 (1945). 
6) D. Iwanenko u. A. Sokolow, Sow. Phys. 11, 590 (19387), 


Sokolo: 


and © 
coordi 
it is e: 
zation 


Here 


differ 
Fı 


in the 


an 
84, be 
at re: 


It 


one n 
ted fr 


then 
yt 
nian. 
O 
dimer 


eigen‘ 
toniaı 


of zeı 
bd, | 
by n 
the 
= 
parti 
while 


= 
+6 
u 
> 
4 
“a 
i 
Alik 
i. Fis. ( 


rest 


(5) 


ried 
yunt 


Sokolov, Ternov a. Loskutov: Transformation Properties of the Spin Pseudo-Vector 243 


and ©, g are the spherical coordinates of the vector k. Choosing an auxiliary 
coordinate system so that the momentum would be parallel to the z-axis, 
it is easy to obtain the following expressions for the components of the polari- 
zation vector ¢,,: 
= NK (Ci 
(Ct C, Ct = K 8, = — 88 Cos = 2% 
6 

iN ky (C71 C, — Cf = Ks =k — §Sind=k, 88 
ka. 
Here N = (Cf C, + 0%, C_,)"' is the normalisation coefficient’), ö the phase 
difference between the amplitudes C_, and C,. 

From Eq. (6) ist is clear that both in case of zero rest mass particles and 

k 
in the ultra-relativistic case (Bı =r> 1), the spin pseudo-vector components, 
\ 

s,, and s,, perpendicular to vector k, tend to zero, while the time component, 


8,, becomes equal to the longitudinal one, ss. On the other hand for the particle 
at rest, the time component becomes equal to zero. 


It is easy to show that 


Cn be kg 


7) Analogous results can also be obtained from paper of Tolhoec’). We note that 
one may take advantage of the covariant expression for spin matrix as explicitly separa- 
ted from the orbital angular momentum 


> „1 
aa =, Kaa 


then we obtain once more the eq. (6) for positive energy particles. We have also: 
K o,y=yt I’, y, where the matrix I’, = ky 0304, commutes with the Hamilto- 
nian. 

On the contrary, the longitudinal s§ and transversal s}, components of the three 
dimensional unit vector possesses different transformation properties; s and 8?, are 


eigenvalues of the operators (* i) and Qs |x le PI, these commuting with the Hamil- 


k 


tonian for the free particle case B- Er It is seen from Eq. (6) that for particles 


of zero rest mass only longitudinal spin component will be different from zero (0, = 
ks, 01,9 = 0). In this case it is more convenient to describe the. transverse polarisation 
by means of the fourdimensional polarisation tensor (Mjs = Hg, My = TE, etc.). For 
the wagers of the particles’ magnetic and electric intrinsic moments we have: 
My = Ko; Ms = ky 883 Ei, = + ksta; & = 0. From the last Eq. it it seen that for 
particles of zero rest mass only the transverse components are different from zero, 
while the longitudinal ones tend to zero. 


8) H. A. Tolhoec, Rev. mod. Physics 28, 277 (1956). 
®) Ju. M. Shirokov, Zhurn. Eksp. Teor. Fis. (Russ.), 35, 1005 (1958); A. I. 


Alikhanov, G.P. Eliseev, V. A. Lyubimov u. V. V. Ershler, Zhurn. Eksp. Teor. 
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i. ©. m polarization vector ¢, is space-like, while the BEER vector 


ky = K 


is time-like, 


In addition it may be noted, that the vectors k, and ¢, turn out to be mutually 


orthogonal: 


For an unpolarised beam of electrons s, = 0, while the phase 6 is a rapidly 
on the average, both Cos 6 
The partial polarization is also possible: for example 0- 


varying quantity, 


so that, 


since 


ky k, = — 


a rapidly varying quantity. 


are strictly fixed. 


the k-axis. 
k 8g, 
mined by quantity 8. 


From the Eq. (6) it is not difficult to obtain 


ky = 0. 


Then we will have: 3 = Ks, & = ky V1 — 
whence it is seen that all the spin pseudo-vector components 


3-dimensional spin pseudo-vector 


[a] = | st + + 33 
For s, # + 1 it is less than unity. 


— 


kg. 


— 83,6, 0, 


me 


— 


3 + - 83). 


$ 3. Transformation of the spin components ER ST 


In this section we shall try to find the transformation law of the spin pseudo- 
vector under the Lorentz rotations. For this purpose we shall assume, that 


> 
in the original (unprimed) coordinate system, the momentum & lies in the 


plane (i.e. the angle 


bitrary orientation. 


We note that these conditions do not restrict the generallity of the dis 
In this case, the y-function in the original system will be of the form: 


cussion. 


ı Sin 
p= 
V2 Cos 
Sin 
where 
_Vı 


In order to obtain the y -function in the primed coordinate sy eee mov ing- 
relatively to the unprimed one along the z-axis with the velocity cB, it is 
necessary to substitute all the unprimed quantities for the primed ones. 


(6, Cos 50, 8in $) 


2 


2 


- 0) while the spin pseudo-vector s has an ar- 


exp(—icKt+ik?), 


and Sin 6 are zero. 

83| < 1, 6 being 
For the total polarisation the quantities 8 and 6 
In this case one of the transverse polarisation vector com- 
ponents may be made equal to zero by rotating the coordinate system around 
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addition it is known (see, for example ref .*)), that under the Loren tz’s trans- 
formations the w-function transforms according to the law 


x 8h 2) y, (12) 
where 
1 B 
chy=-——— and skhy= and = 0, Gs. 


yı- 


Taking accou nt of the relations (11) and (12) and having in mind that 


—icKt+ikr=-—icK't'+ik'r, we obtain the following simultaneous 
equations: 


(af Cy Cos (apf. + by f_) Cy Sin = (Ci Cos 1, 8in$), 
+ b9f-) C, Sin + (dof — by f_) ö(c + 


by f.) C, Cos (ag + bo f+) C_, Sin f_o (ci Cos — + C! Sin 


2 


2 2. 
(13) 
The factor 6 = (L’/L)": on the right hand side, coming out as a result of the 
variation of volume ZL? will not be present in the final results, as well as the 


normalisation coefficient; likewise the quantity Z? is always absent in the \ 
final result (see®)). Here we have taken the following notations: We 


Y Y 


Defining Cj, Cj, K’ and ©’ from the Eq. (13) we can find the spin com- 
ponents in the unprimed system from Eq. (6) where the unprimed quantities 
must be replaced by primed ones: 

(A, — BC os @) 8, + Sin © it (14a) 
V1 — Bp, Cos — (1 (1 — Bi)’ 
[(B, — B Cos 9) 8, — B (1 Sin — (14b) 
(1 — BB, Cos ©) V(1 — BB, Cos 9% — (1 — BF) (1 — fi) 


8, V1— (14c) 


1 — BB, CosO’ 
(Br — B Cos @) + Bs, Sin@ 


In such a way we obtain for the absolute value BL ch >, 


+ sf + = (Bi — B ‘os ) BP) 


(14d) 
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It is seen from (14a) and (17), that the tangent of angle x between the 3-di- 
mensional vectors k’ and 8’ is given by the expresion: 


tga = — 


We assume that in the original coordinate system the vectors $ and & were 
parallel (s, = 1). Then in the primed system the quantity s3 w ill be less than 
unity according to the formula (14a oe 

8 
V (1 — 6 B, Cos @)? — (1 — PP) (1 — BF) 
It means that the perpendicular spin component sj; appears with the longi- 
tudinal one (see (14b)). 


The angle ketween vectors s’ and k’ can be found from the relation 2): 


(1-2) Sin @ 
= — BR, Cos ©) (B, — B Cos ©)" (20) 
i 


This Eq. may be used to investigate the ~-meson polarization. 
As it is known when the z-meson at rest (in the centre of mass system) 
disintegrates into u-meson and neutrino the spin of u-meson must be directed 


10) In arriving at the last eq., we have taken into account that under the Lorentz’s 
transformations, the quantities K’ and k’ transform according to the law: 
K’ = — : (1-Pß,C0os0), k’ = : ya — B B, Cos 9? — (1 — BR) (1 — BP). 

Note, that the same results for the transformation law of the spin pseudo-vector could be 
obtained from the relations (7), (8) and (9) taking account of the transformation law of the 
anti-symmetrical tensor under the Lorentz rotations. 

To demonstrate it we consider the quantity sj, which characterises the doubled spin 
projection on the momentum direction. 

It follows from (9) 


In addition it is known that se in a 
where 
and 
= K {s,Cos © — s, Sin ©}. (16b) 


Substituting (16b) into (16a) and taking account of Eq. (15), we can easily obtain the 
formula (14a). Other formulae (14b), (14c), (14d) may be obtained in the same way. 
The transformation law for the spin pseudo-vector may be pues > — 


Ik k, i’ kt (for 
\ 
details see ")). 


11) A. A. Sokolov, I. M. Ternov and Ju. M. Loskutov, Zhurn. EksP- Teor. Fis. 
(Russ.) 36, 930 (1959). 

12) The analogous , value for tg « = tg (9; — O;) can be obtained from the formulae for 
Cos ©; and Cos ©, given in the ref.™), where ©; and ©; are respectively the angles of 
the momentum and spin vectors with respect to the z-axis in the primed system. 
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along its momentum. On the other hand in the laboratory system, where the 
a-meson is moving. the spin of the #-meson will turn relativ ely to the momen- 
tum by the angle determined by the equation 


ß, /1 — p2 (1 — 62) Sin®, of: 
Pu os (Py, + B, Cos @,,) 
where the angle ©, = © gives the flight direction of the u-meson with velo- 
m2 — m2 
city ce B, =¢ py = Ca = = in the centre of mass system, and cf, is the 


velocity of the z-meson in the laboratory system, equal to minus the velocity 
of the laboratory system relative to the centre of mass one (Bx. = — B). 

The problem of the appearance of u-meson transverse spin component 
in the laboratory coordinate system was investigated also by Askoli!*). Ho- 
wever, the formula given by Askoli for the angle «, 


Sin ©, 
tg, = Ps Pi (22) 
is different frome - equation by the factor 


1+ 8, B,, Cos ©, 2, 
that takes account of the relativistic contraction of the transverse spin com- 
ponents 1). 

In fact, the spin axial vector characterises the circular polarisation in the 
plane perpendicular to its direction. For the longitudinal spin component 
this plane is perpendicular to the particle velocity and therefore the polari- 
zation remains unchanged. On the other hand for transverse spin component 
s the particle velocity will lie in this plane, and therefore the polarization 
will change. 

We note, that according to eq. (20) for the particles with zero rest mass, 
the velocity 6, = 1 and therefore the angle x becomes zero: i.e. the axial 


vector s being parallel to the momentum i (85 = s = 1) in a certain inertial 
system, these vectors will remain parallel (s3 = s’ = 1) in all other inertial 
systems. 

This circumstance allowed one of the authors (see ref., ref.15)16)) to use 
the results of®) to discriminate the neutrino and the anti-neutrino, their 
rest mass being ‚oa to zero, by the different values of the quantities s 
(respectiv ely s = — 1 and s = 1). 


_ 38) G. Ascoli, Z. Physik 150, 407 (1958); see also J. H. D. Jensen, H. Overis, 
Kgl. norske vid selskabs forhandl. 31, N 6 (1958). 

4) Kerimov, who has drawn our attention to the paper of Askoli’), has particip- 
pated also in the disucssion of the section. In ref. 13) rotation of the three dimensional 
unit polarisation vector is described by the formulae: 

83° = 83 Cos xy — 8] Sin a 
= Cosa, + 8$Sin a, 


= 
These eqs. are easily obtained from eqs. (6) and (14). 
15) A. A. Sokolov, Zhurn. Eksp. Teor. Fis. (Russ.) 33, 794 (1957). 
1) A.A. Sokolov and B. K. Kerimov, Ann. Physik 2, 46 (1958). 
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We note i in conclusion that if the polarization of neutrino is characterised 
by its spirality, i. e. by the sense of the rotation of the 6 y component, per. 


pendicular to the momentum &, then in passing from ‘‘right’’ coordinate system 
to the left one, this rotation remains unchanged as was shown in!”), i. e. in 
passing form the “right” coordinate system to the left one, the neutrino is 
not transformed into an anti-neutrino. 

In fact, for a passage, the momentum (being a polar vector) reverses its 
orientation (hight = — kıeit), while the spin axial vector remains unchanged 
(d.aı — A: However, a given direction of an axial vector characterizes dif- 
ferent senses of rotations in the right and the left coordinate systems respecti- 
vely. Therefore, in case of the space inversion, the spirality (i. e. the sence 
of the polarization vector rotation relative to the momentum) will be con- 
served. 


17) A. A. Sokolov, Nuclear Physics 9, 420 (1959), see also A. A. Sokolov, ..Max- 
Planck-Festschrift 1958“, p. 309; A. A. Sokolov and Ju. M. Loskutov, Ann. Physik 
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Bemerkungen 
zur Theorie der Konzentrationsdepolarisation Pr. 
von fluoreszierenden Lösungen 4 


Es wird ein einfaches Verfahren zur Berechnung der Konstanten = an- 


gegeben, die in der Theorie von Wawilow über den Einfluß der Konzen- 
tration auf die Fluoreszenz von Lösungen auftritt. Ferner wird gezeigt, daß 
der Wert dieser Konstanten dem Kugelvolumen vom Halbmesser R, gleich 
ist, wo R, den in der Theorie von Förster über die Konzentrationsdepolari- 
sation der Lumineszenz auftretenden kritischen Abstand zwischen benach- 
barten fluoreszierenden Molekülen darstellt. 


In allen Theorien über die Selbstdepolarisation der Photolumineszenz von 
Lésungen!~*) treten gewisse Konstanten auf, die durch Vergleich mit dem 
Experiment zu bestimmen sind. In einigen Fällen ist eine solche Bestimmung 
mit recht großen Schwierigkeiten verbunden, da hierbei öfter zusätzliche 
Messungen anderer physikalischer Größen notwendig werden. Dies trifft z. B. 
für den Fall der Konstante R,°*) zu, welche in der Theorie von Förster auf- 
tritt, und zu deren Berechnung wir unter anderem die mittlere natürliche 
Fluoreszenzabklingdauer sowie die Uberdeckung von Absorptions- und 
Fluoreszenzspektren kennen müssen. Bei der Mehrzahl von zähflüssigen und 
festen fluoreszierenden Lösungen fehlen jedoch die entsprechenden Angaben ®). 


Kürzlich wurde ein einfaches Verfahren?) zur Bestimmung der Kon- 
stanten R, angegeben, bei dem aus der Kurve der Konzentrationsdepolari- 
sation die Anzahl der lumineszierenden Moleküle in 1 em? abgelesen wird, 
für welche der Polarisationgrad 


3.P. 
0 
6— Py 
1) S. I. Wawilow, J. Physics URSS 7, 141 (1943). 
?) Th. Förster, Ann. Physik 2, 55 (1948). 


3) G. Weber, Trans. Faraday Soc. 50, 552 (1954). 

*) A. Jablonski, Acta Phys. Pol. 14, 295 (1955); 17, 481 (1958). 

5) Die Konstante R, bestimmt denjenigen kritischen Abstand zwischen zwei gleichen 
fluoreszierenden Molekülen, für welchen die Ausstrahlungswahrsche inlichkeit der Über- 
tragungswahrscheinlichkeit von Anregungsenergie gleich ist. 

*) Th. Förster, Fluoreszenz Organise her Verbindungen, Stuttgart 1951, S. 172—180. 

*) C. Bojarski u. A. Kawski, Ann. Physik im Druck. 
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beträgt, wobei P, die Grundpolarisation des fluoreszierenden Lichtes be. Der 


zeichnet®). Unsere Kenntnis von n’ erlaubt es, Ry aus der Beziehung R, = (1) geg 
3 


zu berechnen. Den Ausdruck (1) erhielten wir bei Voraussetzung Au: 


2 


| 4an’ 
gleichmaBiger Verteilung fiir diejenige Konzentration der fluoreszierenden 
Buben bei welcher der Abstand zwischen benachbarten Molekülen gerade 


R, beträgt. Bei der Berechnung des Ausdrucks (1) gingen wir unter anderem Ta 
I von der Voraussetzung aus, daß für Moleküle i-ter Ordnung?) die Ausstrah- “4 
Ps lungswahrscheinlichkeit nach allen Richtungen gegeben sei durch: schein] 
1 Da di 
Nunmehr wollen wir beweisen, daß die obige Voraussetzung in einfacher | lässige 
a Ww eise aus der allgemeinen Gleichung hervorgeht, die Waw ilow}) in seiner Wi 
e Theorie über den Einfluß der Konzentration auf die Fluoreszenz von Lösungen 
= Nach Wawilow gilt hierbei folgender Ausdruck fiir die Ausstrah- 

lungswahrscheinlichkeit Er 


. nt + Ha NT + Ho > 
Konst: 


t bezeichnet hierbei die Fluoreszenzabklingdauer der Molekiile. Sie ist 
von n abhängig und nimmt bei starken Konzentrationen mit wachsendem n Di 
wie folgt ab: 


‘ Fluore 
x 
(4) 125- 
. 

angibt 

Die in den Ausdrücken (3) und (4) erscheinenden Größen A), x, und % Au 


bezeichnen gewisse Konstanten, die von den Eigenschaften des lumines- é 
zierenden Moleküls und von der Umgebung abhängig sind, n ist die Zahl der f — ein 
Farbstoffmoleküle pro em? und z, die mittlere natürliche Fluoreszenzabkling- den k 
dauer für x — 0. Die Erfahrung lehrt, daß in vielen Fällen 2, = 0 ist. 


Unter Berücksichtigung dieser Tatsache erhalten wir bei Einsetzung des Be ‘me 
. x 
Ausdruckes (4) in (3) 
1 1 i 


Wenn wir numehr mit n’ diejenige Konzentration der lumineszierenden Mole- 
bezeichnen, fiir welche 


ist, dann geht der Ausdruck (5) in die von uns vorher angenommene Formel (2) 


über. 


8) Unter „Grundpolarisation‘ ist die Polarisation des Fluoreszenzlichtes der Lösungen 
gemeint, die in Abwesenheit irgendwelcher depolarisierender Wirkungen senkrecht zur 
Schwingungsrichtung der erregenden Strahlung zu beobachten ist. 

%) Die angeregten fluoreszierenden Moleküle, welche unmittelbar von der primären 
Lichtstrahlung ihre Elektrone nanregungsenergie erhalten haben, bezeichnen wir al 
Moleküle nullter Ordnung, und diejenigen, zu denen die Anregungsene art ae von don Mole- 
külen (¢ — 1)-ter Ordnung überging, als Moleküle i-ter Ordnung. | 


r 
"4 
+: 
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ne Der Polarisationsgrad, der dieser Konzentration n’ entspricht, wird durch 
0= f (1) gegeben. 

tzung Aus dem Ausdruck (6) geht hervor: 
erade 
lerem = und — bezeichnen hier die resp. Verhältnisse der Übertragungswahr- 


scheinlic chkeiten mit und ohne Löse "hung zur Ausstrahlungswahrscheinlichkeit. 
Da die Depolarisationswahrscheinlichkeit bedeutend größer als die Lösch- 


2 


° . AR e . T, % 
(2) wahrscheinlichkeit ist, können wir in erster Näherung — gegen — vernach- 
x x 


acher | lässigen. 


einer Wir können also annehmen, daß 
ingen 
trah- 


Ermitteln wir daher aus der Kurve der Konzentrationsdepolarisation die 
(3) Konzentration n’, welche zum Werte P laut (1) gehört, so finden wir die 
T 
Konstante —. 
e ist “2 
mn Die auf diese Weise berechneten Werte von — betragen für Lösungen von 
Ho 

Fluorescein 7,6 - 10-19 cm? und fiir Rhodamin B in wasserfreiem Glycerin !) 
(4) § 12,5 - 10-19 em3, während Wawilow dafür 8,15 - 10-19 em? bzw. 12 - 10-1? cm? 
angibt). 

T, 
id % Aus der Berechnung der Konstanten — geht hervor, daß der Wert von 


I der $ — einen Kugelvolumen vom Halbmesser R, gleich ist, wobei Ry nach Förster 


den kristischen Abstand zwischen benachbarten Molekiilen darstellt Ry = 


; 1) Nach den Messungen von P. P. Pheofilow u. B.J.Sweshnikow, J. Physics 
(5) J URSS8, 493 (1940). 


Gdansk (Polen), Technische Hochschule, I. Physikalisches Institut. fi 


[ole- 
5 Bei der Redaktion eingegangen am 25. Mai 1959. pl 
4 
(6) 


= 
1 
ling- | 
des 
laren : 
r als 
Mole- 
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= 
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Verstärl 


werden 
eine Ga 


Gasentl: 
Theorie der Diffusionswellen. II deren Z 


(Amplifikationsbedingungen in Niederdruckentladungen) | Entladı 


einmal 
Von H. Rother dureh 7 


außerde 
Inhaltsiibersicht peratur: 

Im Rahmen einer vorher!) entwickelten Theorie wird gezeigt, daß diffu. | hältnis 
sionsbestimmte Amplifikationsprozesse (laufende Schichten usw.) in Gas. | Wert 1 
entladungen nur in Strombereichen auftreten, in denen Trägererzeugung oder f bereich 
Wandverluste nicht mehr proportional zur Elektronenkonzentration zunehmen. f Edelgas 
Abweichungen dieser Art in der Trägererzeugung sind möglich, wenn Stufen. f stromst: 
ionisaticn einen merklichen Anteil der Trägerbilanz ausmacht und noch | werden 
keine Sättigung der Konzentration angeregter Atome eingetreten ist. Hier. | regter A 

für kommen Stromstärkebereiche zwischen & 1mA und etwa 10 A in Frage, | hohen 
je nach Gasdruck und -art in Übereinstimmung mit dem experimentellen Schicht 
Befund. Bei sehr kleinen Stromstärken (< 1 mA) treten laufende Schichten Konzen 
. auf, sobald die subnormale Entladungsform », S n realisiert ist, für die nach f 8° gerii 
ticker?) die Wand-Diffusion wegen Abweichung von der Quasineutralität Theorie 
genügen 
Stromb 
Eine nähere Diskussion von Dispersionskurve, Resonanzkurve und Exi- f Gasart ı 
stenzgrenzen ermöglicht Aussagen über eine Reihe von experimentell bisher J stärkeb« 
nur schwer zugänglichen Bestimmungsstücken des Plasmas. existiert 
Bei hohen Drucken ( < 1) existieren keine diffusionsbestimmten 

sti 
Amplifikationsprozesse, wenn ausschließlich Elektronen und positive Ionen fi regter A 
für die Trägerbilanz eine Rolle spielen. Verfasse 
Ein möglicher diffusionsbestimmter Mechanismus, der auch hier erheb-f jn Plası 
liche Feldüberhöhungen erzeugen kann und eventuell für die Erniedrigung f hohe Qı 
der Durchschlagsspannung bei hohen Drucken verantwortlich zu machen ist, f bar wir 


. 
genähert proportional n, (1 Damp» verlaufen sollte. 
e 


wurde kürzlich diskutiert!). lonisier 

Bahnqu 

A. Einleitung merklic 

wait : _. laufend 
In einer vorhergehenden Arbeit!) konnte der lange Jahre hindurch weit- 

gehend ungeklärte Mechanismus der laufenden Schichten qualitativ un Wi 

quantitativ befriedigend dargestellt werden. Ausbreitungs- und Rückkopp- 

lungseigenschaften und die Lage der Resonanzfrequenz ergaben sich in guter oy 
Übereinstimmung mit dem Experiment, sofern die Existenz von nahezu unge ‘ K 
dämpften oder aber amplifizierten Lösungen als realisiert angenommen deuten d 

ER. 4 angeregt, 
1) H. Rother, Ann. Physik (7) 4, 373 (1959) sowie Z. Phys. (im Druck). ir 


2) G. Ecker, Proc. physic. Soc. (London) B 67, 485 (1954). 
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eine Gasentladung erfiillen muB, damit solche Lésungen existieren, daB also 
eine beliebige lokale Störung sich unter ständiger (mitunter erheblicher) 
Verstärkung zur Anode (oder Kathode) hin bewegt. Wir werden zeigen, daß 
hierzu die Existenz eines Energiespeichers erforderlich ist — für normale 
Gasentladungen in Form von angeregten Atomen (oder Fremdzusätzen) 
deren Zahl mit der Elektronenkonzentration anwächst oder für subnormale 
Entladungen in Form eines Diffusionspotentials, das dafür sorgt, daß eine 
einmal aufgetretene Erhöhung der Elektronenkonzentration nicht mehr 
durch Trägerverluste kompensiert werden kann. Damit bei Vorhandensein 
eines solchen Energiereservoirs tatsächlich eine Amplifikation stattfindet, ist 
außerdem erforderlich, daß das Verhältnis von Trägererzeugung zur Tem- 
oH 
aT, 
hältnis von Elektronentemperatur zu Feldstärke x Radius einen bestimmten 
Wert nicht unterschreiten. Diese Bedingungen beschränken den Druck- 
tereich für die Existenz laufender Schichten bei Niederdruckentladungen in 
Edelgasen im allgemeinen auf etwa 0,1—10 Torr. Die Puppschen Grenz- 
stromstärken (/grenz: P const) werden realisiert durch das Unwirksam- 
werden dieses Energiespeichers infolge Sättigung der Konzentration ange- 
regter Atome durch Rück wärtsprozesse (Stöße zweiter Art usw.) bei hinreichend 
hohen Elektronenkonzentrationen. Das Verschwinden (Dämpfung) der 
Schichten bei hinreichend kleinen Stromstärken?) tritt dann ein, wenn die 
Konzentration angeregter Atome gegenüber den Atomen im Grundzustand 
so gering ist, daß Stufenionisation gegenüber der in der Schottkyschen 
Theorie als maßgeblich angenommenen lonisierung aus dem Grundzustand 
genügend klein wird. Falls diese untere ,,Grenzstromstarke bereits in den 
Strombereich der subnormalen Entladungsform fällt (das hängt von Druck, 
Gasart und Rohrradius ab) so sollten laufende Schichten im gesamten Strom- 
stärkebereich von null bis herauf zu den Puppschen Grenzstromstärken 
existieren. Bemerkenswert ist die Tatsache, daß in der Existenz der laufenden 
Schichten erstmals ein deutlicher, experimenteller Hinweis gegeben ist, für 
ein starkes Anwachsen des Wirkungsquerschnittes der Stoßionisierung ange- 
regter Atome gegenüber Atomen im Grundzustand®), das schon früher?) vom 
Verfasser im Zusammenhang mit der Erniedrigung der Ionisierungsspannung 
in Plasmen vermutet wurde. (Es wurde dort angenommen, daß für nicht zu 
hohe Quantenzahlen — etwa bis die Laufzeit des stoßenden Elektrons vergleich- 
bar wird mit der klassischen Umlaufzeit des gebundenen Elektrons — der 
lonisierungsquerschnitt angeregter Atome proportional ist dem Bohrschen 
Bahnquerschnitt.) — Es ist ganz ausgeschlossen, daß ohne die Existenz einer 
merklichen Stufenionisation in einer normalen reinen Edelgasentladung 
laufende Schichten auftreten. Schätzt man aber die Wahrscheinlichkeit 
(vgl. Abschnitt C) für die Stufenionisierung ab unter Zugrundelegung eines 
Wirkungsquerschnittes etwa gleich groß dem des Grundzustandes, so erhält 
man um den Faktor 10 oder mehr zu geringe Wahrscheinlichkeiten. — 

5) Für Neon vgl. etwa H. Achterberg u. J. Michel, Ann. Physik (7) 2, 365 (1959). 

‘) Kürzlich erhaltene experimentelle Befunde von H. Burck (noch unveröffentlicht) 
deuten darauf, daß auch der Stoßquerschnitt für elastisch 


angeregten Atomen ähnliches Verhalten zeigt. 
5) H. Rother, Ann. Physik (7) 2, 326 (1958). 


peraturabhangigkeit der elektronischen Energieverluste sowie das Ver- 
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werden konnte. Im folgenden soll berechnet werden, welche Bedingung: 
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Das vielfach beobachtete Auftreten wilder Instabilitäten in ruhig und 
störungsfrei brennenden Entladungen nach Beimischung von (z. B. organi- 
schen) Fremdzusätzen, sowie Bestrahlungs- und Tilgungseffekte (Penning- 
gemisch) muß auf eine Veränderung der Stufenionisation durch diese Ein- 
flüsse zurückgeführt wurden. Die Löschung der laufenden Schichten in Edel. 
gasen bei Quecksilberzusatz findet ihre einfache Erklärung darin, daß die 
angeregten Edelgasatome nicht mehr durch Elektronenstoß ionisiert werden, 
sondern statt dessen ihrerseits die Hg-Atome ionisieren. 


E? 
Vernachlässigung der Terme von der Ordnung InenU, gegen eins. Der Fall 
TEN, Ue i ‘ 


4nen, U, 


E2 
und Feldstärken (Vorstadium des elektrischen Durchschlages bei hohen 
Drucken). Es seien als Träger nur Elektronen (n,) und positive Ionen (n,) 
wesentlich. Dann läßt sich die Trägerbilanz schreiben in der Form: 


“<€1 gestattet Aussagen über die Verhältnisse bei hohen Drucken 


ön, 


An. Wee) 


Cx 


reichend kleinem Ionisierungsgrad — für Plasmen, in denen Diffusionswellen 
eine Rolle spielen, ist der Ionisationsgrad immer < 1 — die Schwankung der 
Ionentemperatur 7, zu vernachlässigen. Effekte infolge der radialen Ab- 
nahme der Gastemperatur sind im allgemeinen vernachlässigbar und werden 
nicht berücksichtigt. 

Die Energiebilanz der Elektronen können wir dann in der Form ansetzen: 


und Poissongleich 


Hierin 


Es ist: b, = Elektronbeweglichkeit, b, — Ionenbeweglichkeit. E = Elek- 
trische Längsfeldstärke, 7, = Elektronentemp.; k = Boltzmann konstante; 


fy i 
tronen, 
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B. Die Dispersionsgleichung der Diffusionswellen 
\ Wir die Dispersionsgleichung der Wellen ab 1) gemachte 5 
W <opplu eng tralgas i hin- u 
MM wird gesetzt: 
5, U, on 3\, eu, 
b, U, on 3 ou, 
i 


(1) 
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e = elektrische Elementarladung, m, = Elektronenmasse; 


e m. f 
m, | IKT, J Auftinwide (w V IKT, 
oo oo 
Jf dw J WA dw 
0 0 


f dw 

d 43 


fo ist die (normierte) Geschwindigkeitsverteilungsfunktion der Elek- 
tronen, die für große Elektronendichten in die Maxwellverteilung, für 
kleine Elektronendichten in die verallgemeinerte Druyvesteyn-Davy- 
dovsche Verteilung®) übergeht. 

Für Edelgase wird die Trägerbilanzfunktion Z für den zylindersymmetri- 
schen Fall (vgl.')): 


n 


Z= 24 +—2mi— (=) b, U, (1 (Z = 0 für ungestörte Säule). (3) 


Ny 


Die Energieverluste H’ setzen sich außer den Wandverlusten zusammen!) 
oo 
aus den elastischen Verlusten ~ U," za fo 42 w* dw, den unelastischen Ver- 
6 


Na 


lusten 2.» U,» und dem negativen Anteil von den Stößen II. Art pro- 
0 


a 
portional 2%. Die drei Anteile unterscheiden sich wesentlich durch die re- 
lativen Werte ihrer Ableitungen nach T,,: (> 1-2; 8 5; zT) U, = 

e ‘e 
lonisierungsspannung, U, = mittlere Anregungsspannung. Nach kleinen 
m > oO > 
Drucken hin nehmen die unelastischen Verluste gegeniiber den elastischen 


Verlusten wegen anwachsender Elektronentemperatur stark zu. 


Die fiir das Folgende wichtige relative Konzentration = angeregter Atome 


ist mit den in !) verwendeten Definitionen: 


No Im * + mi ‘ 
wo die Volumenrekombination x”, im allgemeinen vernachlässigt werden 
dürfte. 

8) Vgl. S.Chapman u. T.G. Cowling. Mathematical Theory of Non-Uniform 
gases. Cambridge 1952. Experimentelle Anhaltspunkte fürden Übergang zwischen beiden 

] gang 

Verteilungen in Neonentladungen scheinen Davies und Cowcher [(Austr. J. Phys. 8, 
108 (1955)] aus Rauschmessungen gefunden zu haben — s. weiter unten Abb. 4. 

*) Für die Wandverluste ist als Näherung die von Ecker?) unter vereinfachenden 


Voraussetzungen gewonnene Beziehung Pyana (1 ( bU, (1 + Dam» 
n, n, 
4 3eU, U, 
zugrunde gelegt. Ecker gibt für 7 den Wert 7 = > (m, — Masse 


8 Bed, 
der Ionen) an, dessen Ableitung wir uns nicht anschließen möchten. "Fir Edelgase zu- 
mindest scheint dieser Wert nicht brauchbar zu sein (siehe weiter unten). 
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erhalten wir für Wellen der Form », o, v m ei(Kz-ot) mit reellem w nach zustan 
gleichen Methoden wie in!) für v, 0, v die komplexe Dispersionsgleic hung: 
i K® +- [(i +7 ) Oo + +19, + K® 
+ [—k, — (dg + &) w] K? + [—(2 + (5) 
+2a—ih,o = 0*) 
¢ 
Mit den Abkiirzungen F=f 
H=H--2; H, 
2 b, E? ön, b, E? Zur 
n.U: 02 n,U, 02 sch 
T b, E? aU,’ b, E? @n, (6) in ein 
nehmen die Koeffizienten von (5) die Gestalt an: (Die Einführung der dimen- ket 
sionslosen Größen (6) ist die natürlichste, Diffusionsproblemen dieser Art _ 
unmittelbar -angepaBte Darstellung): 
A, = 46,2, + ex 1) H, — (H 4 
1 2 4 
2 2 Durel 
Ae 3 Hy — 37’ Verla 


é;, ¢;, die proportional sind zu o = 


4nen,U, 
&= a (04+ 364 1); & = + +7 6H, —+(20- \) 


> 4 2 Zum 
Um eine Vorstellung von der relativen Größe der einzelnen Terme zu geben, ' Elek 
sei eine Abse ‚hätzung ihrer Zahlenwerte angefügt: Für Maxwellverteilung porti 
und geschwindigkeitsunabhängiges A ist 6 = 2; 6* — 3 (für reale Edelgase 
u *) Terme in w? sind weggelassen. Sie spielen für die normalen positiven Schichten — 
Al 
+ 


: 
. 
= 
i 
4 
= 


uch 


(5) 
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sind 6 und ö* kleiner: für N, wird 09 < + ö, 2); Zp >1(10 Zp < 200), 


i 3 Hi 3 Nm Im 
wenn 1 Torr < p< 10 Torr) Z, m 4g: H, = 7; 4, s = ne wo q,, 
der StoBquerschnitt angeregter a qo der Stoßquerschnitt im Grund- 
zustand vgl.t). 


€. Existenzbereiche laufender Schichten in 


Für Niederdruckentladungen ist generell o < 1 unc ı Die funda- 


mentale Dispersionsgleichung (5) hat dann bei Vernachlässigung der Terme 
&. ¢, die Gestalt: 
F=f—itg = 6, Kt —i 6, K®— (,+id,0) K 
+ 2a-thw=0. (9) 

Zur Auffindung der Amplifikationsbereiche ist zunächst zu bemerken (wie 
schon in!) erwähnt), daß Amplifikation für das Gleichungssystem I—IV nur 
in einem endlichen Frequenzband auftreten kann, bei Amplifikation für sämt- 
liche Wellenlängen von A = 4, bis A = oo z.B. wäre keine quasistationäre 
l.ösung mehr möglich. 

Mit K = k + idergibt sich die Bedingung für Amplifikation eines Wellen- 
paketes längs seiner Fortpflanzungsrichtung: 


CM 
10 
ek (10) 
In Bereichen. in denen — sein Vorzeichen nicht wechselt, ist das Fre- 


ok 

quenzband amplifizierter Störungen begrenzt durch zwei Punkte @, und w,, 
für die d = 0 ist. d = 0 bedeutet aber für (9), daß K reell ist, daß also Real- 
teil / und Imaginärteil g von F einzeln verschwinden müssen, d. h.: 
= 6, K*—h, K?—-2wK+2a=0 (9a) 

g = 0, + 6,0 K*—(a+b,)K+hw= 0. (9b) 
Durch Entwicklung von F an der Stelle w,, K, (oder w,, Ky) ergibt sich der 
Verlauf der Resonanzkurve im Amplifikationsband in 1. Näherung: 


AK = K — K, = Ak+ id= 
of rit 
An 


Zum erstenmal wurden Betrachtungen ähnlicher Art für laufende Schichten 
durchgeführt von Wojaczek®) im Anschluß an Dispersionsmessungen in 
stromstarken Argonentladungen. Wojaczek erhielt aus einem System stark 
vereinfachter Bilanzgleichungen mit den Zusatzhypothesen, daß die ambi- 
polare Diffusion träge sei gegen die Schichtfrequenz und die Störung der 
Elektronenkonzentration An, konstant über den Querschnitt, anstatt pro- 
portional zur ungestörten Dichte n,, die Dispersionsgleichung: 


= 


8) K. Ann. Physik VII, 3, 37 (1959). 
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Aus dieser Dispersionsgleichung leitet er die 
= ( b, E ab, die eine erstaunlich gute bereinstimmung mit dem Beeiche 
von ihm gemessenen Teil der Dispersionskurve liefert, obwohl die gemachten 
Zusatzannahmen nicht im Einklang stehen mit anderweitigen experimentellen 
Befunden?). 

Eine Aufklärung dieses Widerspruches liefert die Diskussion der Dämp- 
fung von (9*). mit 

Man erkennt sofort, daß Real- und Imaginärteil nur gleichzeitig verschwin- 
den können für reelles k, wenn die Gleichung Kun 

2,4\2 2,4\2 E? U* = Eine 
R UT, < < 
mindestens eine reelle Lösung hat. Das ist nur die Fall, wenn “ep 
reelle Lö 
2,4-E-R 
2, 
(9*b) 


Der scheinbare Erfolg der Gl. (9*) rührt daher, daß für die ausgemessenen Ä 
Entladungen die linke Seite von (9*b) mit geeignetem 6+ gerade = 1 war. J} Die Bedi 
Es folgt aber aus (9*b) das paradoxe Ergebnis, daß für kleinere „;--Werte als 5 
die der Entladung von 2A, d.h. für noch größere Stromstärken, amplifi- 
zierte Schichten (mit unendlichem Amplifikationsband!) existieren müßten, 
für kleinere Stromstärken aber keine Schichten möglich sind — der experi- Jen mul 
mentelle Befund ist gerade umgekehrt. 
Dieser Umstand mußte Wojaczek im wesentlichen verborgen bleiben, 
da er zur Lösung seiner Gl. (9*) schon die — experimentell offensichtlich er- 


Da mit 
forderliche — Existenz einer ungedämpften Lösung voraussetzte. Die Über- 
einstimmung des Realteiles von (9*) mit den experimentellen Dispersionskurv en herab, s 
gerade für die speziell benutzte Entladung muß daher als ein Spiel des Zu- durch n 
falls — bzw. der zufälligen Realisierungsmöglichkeit durch geeignete Zusatz- fälle leic 
hypothesen (s. oben) — angesehen werden. Die Arbeit läßt daher wie alle Eine nu 
vorhergehenden keinen Schluß von den Eigenschaften der jeweiligen Ent- 
ladung auf die Existenz diffusionsbestimmter Instabilitäten (laufende Schich- 
ten) zu; und die hier durchgeführte Lösung des vollständigen Gleichungs- ai 
systems I—IV ist unumgänglich. hal 

zeigt sc 

D. Amplifikationsbereiche der Diffusionsgleichungen bereich. 

a) Existenzbedingungen der laufenden Schichten (wegen 

Zur Diskussion der Gl. (9) wird zunächst eine simultane Lösung von Drucker 

(9a und 9b) durchgefiihrt. Alle reellen Simultanlésungen von f = 0, g = 0 lation g 
haben (als Lésung von g (k) = 0) k-Werte auf der Kurve: 

2 cos 1+ —w|=Va, W (12) 
~~ 9) Der Ef Der Effusionsstrom in einer von Schichten durchlaufenen Gasentladung schwankt Hier 
mit der Frequenz der Schichten. one Atom d 


= 
4 
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— yw? 
a + by, do 
cos = — = ; w=— —w. (12a) 
(1 + 36, 3d,Ya, 

Bereiche E< 0) liegen zwischen den beiden (falls vor- 


handen) reellen Wurzeln W,, W, der rr 


flat = BW? ( —w w) = 
mit ; 66 h 
p=? 


Eine zweite Bedingung für die Existenz eines Amplifikationsbandes 
u <w< Wy, ist: |cos p (w,)| < 1. 

Da nach (12) und (12a) stets w W < 1 und w > 0, so folgt: (13) hat nur 
reelle Lösungen (und zwar genau zwei) w, und w, wenn: 


sa, > 0 (by > 0). 


Die Bedingung (14) ist die grundlegende Verunntouny fiir die Existenz lau- 
fender Schichten in Niederdruckentladungen. Sie besagt, daB 


4 
by Z, > (Hy + 1-8) (14’) 
sein muß, oder genähert: i : 
OZ >%» b, (14’’) 


Da mit (14’) gleichzeitig b, > 0, so wai b„ u Schwellenwert für Z, etwas 


herab, spielt aber wegen ~ = S 7 keine entscheidende Rolle, wie man sich 


durch numerische Lösung von (12) und (13) fiir die entsprechenden Grenz- 
fälle leicht überzeugt. (14) stellt natürlich nur eine notwendige Bedingung dar. 
Kine numerische Behandlung von (13) ergibt als hinreichende Bedingung: 


B>e>0, (14) 


wo e je nach Größe von a, im allgemeinen zwischen 0,1 und 1 liegt. (14°) 
zeigt schon, daß bei vorgegebenem Anteil der Stufenionisation der Druck- 
&HU,) 
ou, 
(wegen steigendem Anteil an oe wate Energieverlusten), nach groBen 


Drucken hin durch Abnahme von z , da mit (6), (3) und (4) für Z, die Re- 


lation gilt (unter Vernachlässigung von x n,) 


Ben 


bereich begrenzt ist, nach kleinen Drucken durch das Anwachsen von 


Nm Zmi 


(iz) (9o/9m Zm 1) ( Mp 2; +1) + (15) 


Hierin stellt nach !) &/n, die mittlere Wahrscheinlichkeit pro angeregtes 
Atom dar, ohne Mitwirkung eines Elektrons (effektiv!) in den Grundzustand 


he 
4 
a 
. 4 
n: 
m 
on SEEN 
(13) 
p- 
») 
}) 
n 
Is 
1, 
1 
| | 
4 
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9? 


überzugehen, — — z* ist die Wahrscheinlichkeit für Stöße II. Art und n ol Ng 2, 


No mi 
die Wahrse heinlichkeit für ein angeregtes Atom, durch Elektronenstoß ioni. 


Nm 2 
siert zu werden. —"“"* ist das Verhältnis der Stufenionisationsprozesg 


zur Ionisierung aus dem Grundzustand. (14) ist im allgemeinen erfüllt, wenn: 
i 
Abb. 1 gibt den Verlauf von Z, für Neon schematisch wieder; aufgetzuil 
7 (R\ J 
ist Z, (5) über z (p = Gasdruck). In Abb. 2 ist die Druckabhangig. 
keit der von Achterberg und Michel 1. ec. angegebenen unteren Existenz. 


& 


bff 


5 


2 


0" 


g 


Abb. 1. Existenzbereiche von laufenden Schichten in Niederdruckentladungen ak 


Funktion der Elektronenkonzentration Nn.m di (für Neon). ------- Abfall von Z, 
bei Maxwellverteilung, — — — ß = 0,3 (siehe (14)) 


. 
grenzen der laufenden Schichten dargestellt. Aufgetragen ist 


const-n, über p. Die ausgezogene Kurve stellt mit (20) den theoretischen 


Verlauf dar fir ae 1. Bei Druyvesteynverteilung der Elektronen 


2 
Um+Umi Ne A 
. z 
ist wegen e wp? der Quotient’ pP mit y 1,5. 
N, 0“i 
Die gestrichelte Kurve würde den Existenzgrenzen bei Max well verteilung 
Nm Zmi 


entsprechen, für die nahezu unabhängig sein sollte vom Druck. Der 


experimentelle Befund spricht deutlich für die Druyvesteyn verteilung 
(Die Kurven wurden für p = 3 Torr zur Deckung gebracht.) 
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Die in ?) angegebene Druckabhängigkeit für die Existenzgrenzen des sub- 
{ U, 
. „Enormalen Entladungsbereiches mit dem Eckerschen Wert 7 = = a Reb 
ioni- p #60, 
(s. Fußnote *)) ergibt experimentell gerade den entgegengesetzten Verlauf 
(= p-’, ¢ > 0) als theoretisch zu erwarten. Man muß jedoch berücksichtigen, 


daß dieser Wert n nur für Hg-Ent- 


enn: ~ 
ladungen (nach Rechnungen von Me [mA Torr mv] 
Beker 1. e.) zu korrekten Ergeb- JP 
16a) nissen führt. Bei Edelgasen liegt 7 
nicht in der richtigen GréBenord- ER 
L6b) | Pung- In den Eckerschen Wert 
fürn geht die (unbewiesene) Hypo- 
these n, ~ n, ein und außerdem eine 0% Moxwellverteilung 
Randbedingung für den Diffusions- 05 
igig- | strom, gegen die wir einige Bedenken 
enz.f haben. Eine‘ nähere theoretische ost 
Untersuchung dieser beiden Punkte ( 
erscheint uns erforderlich, um die auf- He 
Nach großen Stromstärken hin or 2 
gilt die von Pupp’) gefundene 1 2 3 10 P [Torr] 
Gesetzmäßigkeit: Abb. 2. Niederstromgrenze der laufenden 
- p® = const, (17) Schichten als Funktion des Druckes. « Ex- 
die perimentelle Werte nach Achterberg u. 
woe ale erte hat: Theoretische Kurve fiir 
He Ne A Kr Druyvesteynverteilung, ~-- Theoreti- 
4 1 0,7 05 Havel sche Kurve für Maxwellverteilung : 
Nach (16a) folgt mit (17) für die „Löschungswahrscheinlichkeit“ — . 
0 » 
pro angeregtes Atom die Relation: j 
E2 
cons 
2U... He Ne A Kr 
7, == mi a — 
woy=e—1- u, ' y hat die Werte y = 04 04 02 01° 
Die ~ — erfolgende Strahlungsdiffusion hat bei leichteren Atomen offen- 
D . . r 
-®§ bar einen größeren Anteil an den ‚‚resonanten‘“ Verlusten x als bei schweren, 


hen} wo sie nach (18) kaum eine Rolle spielt. (be: p-E ist zwischen 1 und 10 Torr 
nen 
genähert konstant, U2 nimmt mit p zu.) 
Die von Ac hterberg und Michel ].c. angegebenen Existenzgrenzen 
laufender Schichten nach kleinen Stromstärken hin treten auf, wenn gleich- 


zeitig erfüllt ist (16a und b): 


ung = 
Der a) = <1 und 
mo 
ng. 


) Handbuch d. Physik, Bd. XXII, 8. 140, 141 


| 
H. Rother: Theorie der Diffusionswellen. II | 
b) < (19) 
, Berlin— Göttingen—Heidelberg 1956. 
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Ob schichtenfreie Strombereiche zwischen normaler und inch Ent. 
ladung überhaupt auftreten, hängt davon ab, ob der Eckersche Parameter, 
eine Überlappung des subnormalen Bereiches 7 = x, mit dem Bereich (16a) 
ausschlieBt oder nicht. Die Grenzstromstarken schichtenfreier Nieder. 
strombereiche (für N,:1—10 mA) ergeben aus (3) mit den genähert bekannten 
Koeffizienten 2,, und 29. unter Benutzung von « nach (17) zwei Bestimmung. 
gleichungen für z,,; und «, aus denen sowohl die Strahlungsverluste als auch 
die experimentell" “und theoretisch) noch völlig ungeklärte®) mittlere Ioni. 
sierungswahrscheinrichkeit z,,, abgeschätzt werden können. Leider fehlen 
bei den von Achterberg und Michel durchgeführten Messungen Werte für 
b,. E und U, völlig. Eingehende Messungen dieser Parameter in den Grenz. ee 
strombereichen ließen interessante Aussagen über die atomaren Daten der 


~ 
T 


: Entladungen zu (die anderweitig nur schwer zugänglich sind), zumal Existenz. a 3 
r grenzen und Dispersionskurven mit großer Genauigkeit gemessen werden} Param 
u können, wobei die der linearen Theorie nach (9) entnommenen Dispersions- 1 \ 
u kurven z. B. im Neon wegen der sehr geringen Amplitudenabhängigkeit der («. W) 
E Wellenlänge auch bei starker Amplifikation eine brauchbare Näherung dar- 
2 stellen dürften. Die Untersuchung der Existenzgrenze nach kleinen Strom- mit ¢ 
“ stärken hin läßt zum anderen auch die Diskussion empfindlicher Einflüsse troner 
7 von Fremdzusätzen oder äußeren Einflüssen (Strahlung) auf die Konzen- u 
tration angeregter Atome zu. 
Der experimentellen Bestimmung direkt zugänglich sind: Resonam-} ¢ 
u frequenz, Wellenlänge (die Dispersionskurve in einem schmalen Frequenz- Pye 
m 


bereich in der Umgebung der Resonanzfrequenz ist durch Veränderung der “ 
Rohrlänge meBbar). Von den schwach amplifizierten oder gedimpften Wellen 
in der Umgebung der Grenzstromstärken: Dispersionskurve und Resonanz- ergetis 
kurve in einem relativ großen Frequenzbereich (je nach Steilheit der Re- 
sonanzkurve). 

Da die Relation (12) an den beiden Grenzen des Amplifikationsbereiches vertei 
mit der tatsächlichen Dispersionskurve übereinstimmt, kann sie als brauchbare 


4 


Maxy 


Interpolationsformel benutzt werden. Ihre Korrektion zwischen beiden Pu PE 
Koinzidenzen w, und w, kann der Gl. (21) entnommen werden. Wir schreiben ze. 
(12) in die für eine numerische Auswertung geeignete Form (vy = Frequenz, laufen 
s = Schichtabstand): 
W (w) (1 W (w) Ist u 
ae U, Zr “4 (Bem 
Die Funktion des Parameters w ist für = 0,07 (entspricht 
(w) 09% stark 
p ~1Torr) und Fa = 0,22 (entspricht p ~ 10 Torr) in Abb. 3 dargestellt. me 
(Man erkennt nt ye die qualitative Übereinstimmung mit den in?) angegebenen Resor 
Dispersionskurven — insbesondere Druckabhängigkeit der Neigung und Grok 
Stromstärkeabhängigkeit von Smin fiir » = 0.) ment 
Der Koeffizient x soll den Ubergang erfassen von Maxwellverteilung Die \ 
(x = 1) in eine Verteilung nach Art der Druyvesteyn-Davydovschen strich 


. 
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(x > 1). Dieser Übergang sollte nach den Messungen von Cowcher und 
Davies l. c. im Druckbereich zwischen 1 und 10 Torr zwischen 100 mA und 
1mA stattfinden, entsprechend dem anomalen Verlauf der Elektronentem- 
peratur in diesem Bereich (Abb. 4), der nur dadurch erklärt werden kann, daß 


T° | 
2 5 r 
4} 
3k 
Abb. 3. Dispersionskurven mit p als Abb. 4. Elektronentemperatur in Neonent- 
Parameter in reduzierten Einheiten ladungen nach Davies u. Cowcher l. “il 
(w, jr); (Auf der Abszisse lies w statt = 


mit abnehmender Elektronenwechselwirkung die Zahl der heiBen Elek- 
tronen sehr schnell abnimmt (vgl. auch Abb. 2). wenn die Max well verteilung 
m, 3m 
4M 


e ‘Te übergeht in die Verteilung e ) (Chapman-Cowling]. 
Bei Druyvesteynverteilung ist die Stufenionisation wegen n„2 


mi 


Umi 
e !: gegenüber der Ionisation vom Grundzustand 2; ~e 
. . Nm 2mi 
ergetisch stark bevorzugt. Es wird in (16a) mit ( am ) ==] bei 
2; /Maxwell 


Maxwellverteilung mit denselben Entladungsdaten bei Druyvesteyn- 
Um Umi 


=e % »>1. Bezeichnen wir die zu den 


. Nm 
verteilung : ( = 
N, 2; /Druyvesteyn 
Puppschen Grenzstromstärken gehörige Elektronenkonzentration mit 7%. max; 
so wäre nach (16a) die minimale Elektronenkonzentration 7, min, bei der noch 


laufende Schichten auftreten bei Maxwellverteilung: 
n, 
Ne min = R \2 20) 
(53) 


u ist ein Zahlenfaktor, der sich aus (13) bestimmt; auch in mittleren Bereichen 


ist u < 10. z,, ist aus (3) und (4) als Maximalwert für 2: , entnommen. 
No 


m mt 
(Bemerkenswert ist hieran, daß wegen des großen Wirkungsquerschnitts q,,; 
der angeregten Atome bei Strömen oberhalb der Puppschen Grenzstrom- 
stärken die ‚effektive Ionisierungsspannung‘‘ gemäß dem Boltzmannfaktor 

Um 
e “ von z, praktisch identisch wird mit der Anregungsspannung des 
Resonanzterms.) Nach (20) sollte n.min theoretisch höchstens ein bis zwei 
Größenordnungen unterhalb mmax liegen, nach!) liegt aber memin experi- 
mentell bei Neon z. B. drei bis vier Größenordnungen unterhalb n.max- 
Die Werte n,min nach (20) sind bei den zugehörigen Drucken in Abb. 4 ge- 
strichelt eingezeichnet. Die untere Grenze für n,. bei der Maxwellelektronen 
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noch gerade imstande sind, laufende Schichten anzufachen, fällt beim Neon 
offenbar genähert zusammen mit dem anomalen Anstieg der Elektronentem- 
peratur, der nach Davies und Cowcher l. c. dem Übergang von Maxwell. 
in Druyvesteynverteilung entspricht. Es wäre interessant, experimentell 
zu suchen, ob bei anderen Gasen ein bis zwei Größenordnungen unterhalb der 
Puppschen Grenzstromstärken eventuell ebensolche ,,Homogenitatsfenster“ 
auftreten, wie sie zwischen normaler und subnormaler Entladungsform ge- 
funden werden?). 


b) Resonanzkurve und Dispersionskurve 
Die Entwicklung von (9) an den Stellen &, bzw. w, mit d(@,) = d(@,) = 0 
gibt eine erste Verbesserung der Interpolationsformel (12) K, = k,, wo k, die 
lösung (12) von g (k,, ©) = 0 ist. 
[EF (w,, k)/öw],, _ 


1k+i Ad = = =K (w) — K (w,,) = K (w) — k, = [AF (w, k)/@K | 
K = 


m = 1,2. 


(20) 
Mit den Lösungen k, und k, von f (h,. w) = 0 (siehe Abb. 4) und nt = 


at 

"=< 0: cos y = = 
Va | ck) | \é 


wird: 


a) Ak= (sin + y Aw 


b) Ad = + sin y cos ae). “AY 


Man entnimmt (21a), daß d > 0 wird, sobald die Kurve k, (w) bei «, und 

@, geschnitten wird von der Kurve k,(w), äquivalent mit der Bedingung (13). 
Für © < o, und @ > @g ist d negativ. (21a) besagt, daß die Dispersions- 
kurve s (w) gegenüber der Näherung (12’) eine Wölbung nach oben erfährt und 

zwar wegen sin? y (7) 7 um so stärker, je höher der Druck ist (gestri- 

e e 

chelte Kurven in Abb. 3) — in ausgezeichneter Übereinstimmung mit der 
Gestalt der von Michel und Achterberg I. c. angegebenen Dispersionskurven. 
Das Maximum der Resonanzkurve d (») liegt im Rahmen der in A@ linearen 
Näherung (21) beim Minimum von f(k,(®),®) nach (13). Eine genauere 


(21) 


Form der Resonanzkurve läßt sich im Falle 7< 1 angeben. Wir entwickeln 
(9) in der Gestalt: 


oF 
(K, w)—F (k,, = [(k k,) id] 7 


Wegen F (K,; m) = 0 und F (k,, w) = f (ky, ©) erhalten wir 


d = — f ag 2? 


wobei f; 9 
(23) ist fü 
ladungen 
Maxima ı 
stimmen | 
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wobei f, g usw. mittels (12) alles bekannte Funktionen von @ sind. d(m) nach 
(23) ist für die von Michel und Achterberg untersuchten Niederstroment- 
ladungen in N, in Abb. 4 dargestellt. Die experimentelle Existenz zweier 
Maxima und die relative Verschiebung ihrer Höhe mit zunehmendem Druck 
stimmen qualitativ mit (23) überein. 


Für Entladungen mit großen d-Werten (4 zZ 1) ist (22) nicht brauchbar. 


Hier ist experimentell aber ohnedies nur die Lage des Hauptmaximums (bei 
hinreichend langem 
Rohr praktischidentisch 24 
mit der .,Resonanzfre- aw 
quenz‘‘ vgl.!)) bekannt. N 
Wir indentifizieren das 
Resonanzmaximum hier 
in der Näherung ( 2) a 
mit dem Maximum ay 
—f in (13)*) 

Wir unterscheiden 
für die Lage der Reso- ah a | 
nanzfrequenz », zwei 
Strom bereiche, die stetig 
ineinander übergehen: 
Den Bereich normaler 
Stromstärken (z. B. für 
p=5Torr: J = 20mA), 
in denen die Elektronen 
genihert Maxwellver- 


EN 


€ 


ıI FL; 1 


Abb. 5. Resonanzkurven. a) ein Maximum, b) zwei 
: > > Maxima. Die ausgezogenen Kurven II entsprechen jeweils 
teilung besitzen. Hieı dem höheren Druck (nach (23’)) 

liegt das Maximum 


weitgehend unabhängig von der Stromstärke bei w = 1,5 (W = 0,55) und 


man hat für die 
„—»E U, , 82, 
u,(1 +: ") 
i( Zr 
(24) stimmt gut mit der Puppschen Relation überein!) und liefert für 


Z° 0 (genähert konstante Wurzel) die von Donahue und Dieke?®)!%) 
7 a 

gefundene Proportionalität zwischen Resonanzfrequenz und Feldstärke. 
In s, gehen bei starker Amplifikation nichtlineare Effekte ein. Hier ist im 
. . . . . 
allgemeinen weniger gute Übereinstimmung zu erwarten. Für die in ®) unter- 
suchten Grenzstromstärken ist die Übereinstimmung bis auf einige Prozent 


*) Für beliebiges hd ergibt (9) direkt die allgemeingültige Relation: 


k 
(=) % 0 — (dg hy — 25, — dg hg) kt + (2 (6, —1l)a—2b,—h oh,)k? + 2 
") T.Donahue u. G.H.Dieke, Physic. Rev. 81, 248 (1951). ae lan: va 
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genau. In dem Strombereich unterhalb 20 mA (wieder auf p = 5 Torr bezogen) 
in dem die Elektronentemperatur stark ansteigt (vgl. Abb. 4), verschiebt 
sich das Maximum von — f in (13) immer mehr nach kleineren Werten », 
und », wird mit abnehmender Stromstärke abnehmen (gerade umgekehrt wie 


bei höheren Stromstärken). Da gemäß (12’) » ~ w und s ~ so wird die 


W (w)’ 
zugehörige Schichtweite s, gleichfalls abnehmen in RO A mit dem 
von Michel und Achterberg 1. c. gefundenen Verhalten. Leider sind die 
Angaben dieser Autoren hierzu und die Kenntnis der Größen 6, - E, U, nw 

ungenau, so daß quantitativeVergleiche noch nicht durchgeführt werden können. 

ig Aus (12) folgt für die Phasengeschwindigkeit V der Schichten: 

w 

ic W (wo) 
Diese Relation gilt wegen der Amplitudenabhängigkeit von s, streng nur 

im Bereich der Grenzstromstirken (linearisierte Theorie). 

Als Zahlenbeispiel betrachten wir zum Schluß die in *) untersuchte Neon. 
entladung: 

p= 3Torr; R=0,5cem; J, = J grenze = 3A — nach Abb. 4 muß die 
Maxwellverteilung realisiert sein bis J, ~ 30 mA. 


V=4»=4 b, E. (25) 


or. V 
T,ı © 22000°K; = 2,5—; b,, 2-10 10% em? 
cm sec 


© 25000°K; E, = 5,6 —; by, = 2: 10 © 0,5 - 1011 em) 


cn 
Für J, folgt aus der Relation (3): (—) 6, U, = 2; + 7 mi, Sowie aus 
0 


(4) und (15): x 3- 10sec! em-3; —" 10-%. Für das Verhältnis der| 


Jonisierungsquerschnitte ergibt sich 20 (z%, wurde nach der bekannten 
02 
Relation 


- - U . 
Zo, = 1,35: 108a - p (Torr) - U,YU, e (0.5 + (U,, U, in Volt) 
berechnet mit a, = 0,055 für Neon). 

Aus der Existenz der laufenden Schichten bis zu J = 1mA < 30 mA 
folgt, daß wegen (15) mindestens bis zum Einsetzen der Druyvesteyn- 
verteilung!) bei J, = 30 mA der Quotient — > 1 sein muß. Hieraus 
folgt ar = 20 im Einklang mit den Berechnungen aus Jgrenze = Jı- 
gleich dem Verhältnis der 
F 1 
Quadrate der Bohrschen Bahnen, die wir proportional iD 
Ui 


konsistent angesehen werden. 


Setzen wir zum andern das Verhältnis os 


1 5 
bzw. annehmen’), 


so folgt fiir Neon » 18,5. Die Ergebnisse können als durchaus 


_ 12) Der von Davies und Cowcher 1. c. aus dem anomalen Anstieg von U, postulierte 
Übergang der Verteilungsfunktion in die Druyvesteynverteilung, wird hier völlig un 


und 30 mA. 


abhängig experimentell bestätigt durch die Existenz von Schichten im Bereich zwischen! 
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Die Betrachtung des Falles m *<1 liefert keine physikalisch in- 


teressanten Ergebnisse. 
Aus (5) ergibt sich die Dispersionsgleichung 


F =f +ig=i6)K® + (6, —i6,w) K* + (— 6, +iH,) K® 


—iH,o K*—H,wK =0 
(26-1); ta, — 43; 


09 


9 
H, = 3 H, > 37» 
Amplifikationsbereiche können nur zwischen zwei simultanen, reellen Lösungen 
der Gleichungen: RR 
g = K* (6, K® — 6,0 K? + H, K — H,o) = 0 on 
liegen. Eine gemeinsame asymptotische Lösung K = 0 fiir wm > oo existiert, 
die jedoch für 0< w < oo gedämpft ist, und nur für sehr kurze Wellen- 
längen nahezu ungedämpfte Wellen ergibt. 
Da f = 0 nur reelle Lösungen hat mit K >o, wo für > co: K = 
gilt, so folgt wegen y < 2: 
ow H, — H, K 1 EEE. y H, 
dy (57-3) 0 <—. 
Die reelle Lösung k, von f (k,) = 0 ist stets größer, als die entsprechende 
reelle Lösung k, von g (k,) = 0: SON. > 
Folglich kann mit dem hier behandelten Mechanismus keine diffusions- 
bestimmte Feldüberhöhung zustande kommen. Es pflanzen sich lediglich sehr 


kurzwellige Störungen nahezu ungedämpft mit einer Geschwindigkeit v = 
b,E fort?3). 


k,>k, für w< oo. 


3) Daßsich Fluktuationen von Bereichen kleiner als die Debyelänge nahezu unge- 
dämpft mit » = b,E von der Anode zur Kathode bewegen, ist nicht auf den Fall 
4zen,U, 

EB < 1 beschränkt, sondern folgt aus (5) generell für jede Art von Gasentladung. 


Berlin W 8, Physikalisch-Technisches Institut der Deutschen Akademie 
der Wissenschaften. 


Bei der Redaktion eingegangen am 7. Juli 1959. 


x 


N 


en), 
iebt 
u, 
die 
dem 
die 
nur 
nen. 
(25 
nur 
eon- 
die 
i} 
aus 
mA 
# 
end), 
haus 
‚lierte 


EN | 


Er Free Volume Approach to Polyisobutylene 


: melt Viscosity under Pressure’) 


By Harkrishan Singh 


The Doolittle melt viscosity expression 7 = A e®'! has been applied to 
published data on polyisobutylene to calculate a parameter a, which is de- 
fined as the ratio of relaxation times at any pressure P to those at a reference 
pressure. To check the results, the parameter a, has also been calculated by 
the method of reduced variables. The agreement between the theoretical and 
experimental values of a, is found to be very good. Also the treatment se- 
parates melt viscosity of polyisobutylene into two factors, A dependents 
only on structural parameters, and f, which is determined by the specific 
volume is a function of temperature and pressure. 


The effects of temperature on the flow properties of a high polymer of a 
given molecular species have been described through the variation of a sui- 
tably defined free volume?“). The free volume description has also been 
successfully used to account for the effects of high pressure on the mechanical 
properties of a particular polymer isobutylene®). Further the Doolittle®) 
free volume theory of viscosity has been used to express the temperature 
dependence of mechanical phenomena in the glass transition range of amor- 
phous polymers, in terms of the parameter a, which is defined as the ratio 
of the relaxation times at any temperature T to those at a reference tem- 
perature. 


It is the purpose of the present paper to apply the Doolittle equation of 
viscosity and dilatometric data, on polyisobutylene in order to determine 
the influence of pressure on the parameter a, which is defined as the ratio of 
the relaxation times at any pressure to those at a reference pressure for any 
fixed temperature. 


1) Part of this work was submitted as Dissertation to the Graduate Schcol of the Uni- 
versity of Texas, USA in partial fulfilment of the requirements for the degree of Doctor 
of Philosophy. 

2) T.G. Fox and P. T. Flory, J. Appl. Phys. 21, 581 (1950). 
3) F. Bueche, J. Chem. Physics 21, 1850 (1953). Su lag 
4) Ferry, Grandine, and Fitzgerald, J. Appl. Phys. 24, 911 (1958). 
5) H. Singh and A. W. Nolle, J. Appl. Phys. 30, 337 (1959). 
6) A. K. Doolittle, J. Appl. Phys. 28, 236 (1952). 
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Theory 
According to Doolittle’), the melt viscosity of polymers is given by the 
expression 
n = A eBil (1) 
where A and B are constants independent of temperature, and f, the relative 


free volume, is —y °, where V, is the occupied volume. Following Doo- 


little, the occupied volume is defined as “the limiting volume to which a 
real liquid will contract if it were to continue to behave as a nonassociated 
liquid without change of phase all the way to absolute zero”. Pressure as 
well as temperature will affect /. The effect of temperature has been already 
studied. We will study the effect of pressure. 

In the pressure range above the atmospheric pressure, we will assume 
that at a fixed temperature, the specific volume of polymers is a linear function 
of pressure over a certain range. Then the specific volume of a polymer at 
any given temperature and pressure is given by the expression. 


V=V, + ap — T,) + qp(P — Po) 
and the relative volume is 
where P, is the hii pressure. The expression for viscosity then becomes 
B 
log = log A + (4) 
1 dV 1 dJ 
2.308 | Vo dT (T Te) dP (F I 


This treatment indicates that the pressure dependence of viscosity can 
be related to the pressure dependence of volume independent of molecular 
weight. Thus the following relation will apply to a wide range of materials in 
the liquid state y 

: 
Here V, and V, are the specific EEE at the reference pressure and any 
other pressure respectively, the temperature remaining the same. 

An interesting result that follows from this analysis is that it is possible 
to divide the melt viscosity into two terms; one the constant A, is a function 
of molecular weight and other structural N en — the other the 
exponential term which is a function of free volume. = 3 


log a, = 


Application to Polyisobutylene : 


Williams?) has shown that for polystyrene the vaue of Bis essentially 
constant and is equal to 0.91. As it is found to be independent of temperature 
and molecular weight we may assume it to be independent of pressure and 
also have the same value for polyisobutylene. The occupied volume V, can 
be computed from William’s relation Ris 

0 
7) A.K. Doolittle and D. B. Doolittle, J. Appl. Phys. 28, 901 (1957), 
°) F. Bueche, J. Chem. Physics 20, 25, 599 (1956). 


*)-M. L. Williams, J. Appl. Phys. 29, 1395 (1958). 
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In Equ. (5) v ’ wna V, represent the specific volumes at reference pressure and 
any other pressure respectively. In the present investigation atmospheric 
pressure has been used as the reference pressure. 

Using relation(5) the parameter a, has been computed for various pressures, 
The values of V, and V, have been taken from the work of Singh and Nolle')., 
The density vs. pressure curve in their work may be split up into two parts; 
the first part up to a pressure of few hundred pounds per si shows a very rapid 
rise and the second part indicates a linear increase of density with pressure till 
the pressure becomes very high and the density changes become very small. 
Only the linear portion of this curve has been used in the present investigation. 
It is likely that the apparent initial rapid rise corresponds to a small percen- 
tage of air content in the 
sample. Hence the extra- 
polated value of density 
at atmospheric pressure 
has been used for caleu- 
lating V,. The results of 
Pressure Doolittle Reduced these calculations are 
in Psi Viscosity Equation variables given in Table 1. As an 
independent check on 


Table 1 

Values of log a, for Polyisobutylene from Doolittle, 
Viscosity Equation and from the method of Reduced 
Variables 


B | B- these results, values of 
3 ll 12 14 PR: 
= 0.21 0.23 0:24 the parameter a, have 
9000 0.30 0.33 0.34 also been calculated by 
12000 0.38 | 0.42 0.43 the method of reduced 
15000 0.45 — 0.50 | 051 variables. 


Applieation of the Method of Reduced Variables 


Ferry!)!) on the basis of some simple assumptions regarding stress 
relaxation in polymers introduced the idea of reducing the simultaneous 
variation of rigidity with frequency, temperature and density to a single empiri- 
cal function. 

On the assumption that the behaviour of a linear polymer in small shear 
deformations can be represented by an array of springs coupled by viscous 
elements, we may represent such an array by a number of Maxwell elements 
in parallel in which each spring. @, is relaxed by a dashpot r,@,. Then the 
real part of dynamic rigidity @’, the real part of dynamic viscosity 7’ and 
the steady flow viscosity n are given by 


Now if we make the assumption that all relaxation mechanisms have the 
same temperature dependence so that when the temperature is changed from 
a standard value 7, to 7’, each r, is multiplied by the same factor a, and that 


10) J.D. Ferry, J. Amer. Chem. Soc. 72, 3746 (1950). 4 
u) J.D. Ferry and E. R. Fitzgerald, J. Colloid Sci. 8, 224, (1953), 
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the real parts of bulk and shear moduli and all the relaxation times have the 2 
same pressure dependence such that when pressure is increased froma standard = 
value P, (atmospheric in this case) to P, then @’, x’ and each 1; is 2 
bythe same empirical factor a,, the various unniiiins become 


2,2 
(©) = (a, ap) = 
= Ap ap no (w ap ap) 
(@) = (ay ap) & 
x’ (w) = ap xo (11) 
a 


where Go, 40. 79 and x are the quantities at standard pressure. Here it may be 
mentioned that no correction has been applied to these quantities as the Ar x j 
perature dependence of these is not known. But as the temperature range 5% 
covered is not large, it may be assumed that the correction is not large enough 
to affect the results appreciably. 

The relations (8) and (11) combined ar M'a 


M’ (w) = ap [ +? (w ap 
Thus = when plotted against w a, a, at various temperatures and pressures _ 

p 
should superpose to give one composite curve. This has been done and asa 


result it has been found that polyisobutylene gives one smooth function. 

To determine the values of a,, log M’ is plotted against log w a, at various 
pressures and is shown in fig. 1. These coordinates may be converted to 


27000 
18000 
15000 
12000 
9000 
60 65 70 #5 60 65 90 


log wa, 

Fig. 1. Jog M’ plotted against log wa, for poly ne at various pressures (psi) 

at 4 mc. for empirical determination of a, 
og (= ~) and log (wapa,) by substracting and adding log ap, respectively, 


in effect i is equivalent to shifting the experimental curves through different 
extents in the direction of a line with slope of —1. Logarithmic plots have 
been used because the relative and not the absolute precision of the data 
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is roughly constant over the ranges covered. A line with a Ben of —lis 
drawn through some convenient reference point on the experimental curve for 
atmospheric pressure (the actual lines drawn are shown in the figure). The 
experimental curve at atmospheric pressure is drawn on a transparent paper 
and the paper is slid along the 45° line until the curve at the next higher pres. 
sure overlaps the atmospheric pressure curve and the second curve is traced, 
This process is repeated to include all the different pressures and the resul. 
tant single curve is shown in fig. 2. The total distance moved along with the 


OF 
0 © m © 0 & 


Fig. 2. log = vslogwa,.a, for polyisobutylene for pressures upto 21.000 psi at 4 me. 
horizontal or the vertical direction gives log a,. Values of a, for polyisobutylene 
were taken from Ferry!) et al’s paper. Values of a, obtained in this way 
are given in the table 1 column 4. 

Comparing the results of a, obtained by the two methods, we find that the 
values obtained by using Doolittle equation are consistently less than those 
obtained by the method of reduced variables. Here it may be recalled that in 
the calculation of a, from Doolittle equation, the value of the constantB 
was taken to be equal to 0.91. This value of B was obtained by Williams 
from the experimental data by the application of linear least square pro- 
gramme. On the basis of certain theoretical considerations one is led to conclude 
that B should be equal to one!’). Taking this as the correct value of B, we 
get the values of a, given in column three of the table. These values are in 
very good agreement with those obtained by the method of reduced variables. 
Thus we see that, if we use the experimental value of B which is more under- 
standable, we get approximate agreement between the theoretical results 
predicted by Doolittle equation, and the experimental results obtained by 
the method of reduced variables. On the other hand if we regard B equal 
to one, as suggested by theoretical considerations, we get perfect aggreement 
between theory and experiment. 


The author wishes to thank Dr. M. Pancholy for useful suggestions and 
Mr. P. D. Gupta for carrying out some of the calculations. 


2) J.D. Ferry, E. R. Fitzgerald, M.F. Johnson and L. D. Grandive, J: 
Appl. 22, 717 (1951). 
" a“ 13) M.L. Williams, J. Physic. Chem. 59, 95 (1955). 
rs er. Chandigarh-3 (India), Department of Physics, Panjab University. 
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Über den Ferrimagnetismus nichtstöchiometrischer mit i 


An synthetisch hergestellten Eisensulfiden der Zusammensetzung FeS, +, 
(0< x < 0,15) wurden magnetische Untersuchungen durchgeführt. Dabei 
interessierte besonders der ferrimagnetische Bereich (0,10 << n < 0,15). 
Die experimentell erhaltenen Ergebnisse lassen sich mit einer Ordnungs—Un- 
ordnungs-Theorie (Bragg-Williams) verstehen, wobei angenommen wird, 
daß bei 7’ = 0 sowohl die Fe+**-Ionen als auch die Leerstellen nur Plätze in 
dem einen der beiden Untergitter besetzen. 

Versuche, in dem ferrimagnetischen Sulfid FeS,,,; das Eisen teilweise oder 
ganz durch Kobalt, Nickel oder Kupfer zu ersetzen, hatten nur beim Kobalt 
Erfolg. Die Mischsulfide sind bis Fe, ‚Co, 5S, 1; ferrimagnetisch, bei höherem 
Kobaltgehalt antiferromagnetisch. 

Vorsichtige Oxydation der Sulfide führt dazu, daß zunächst etwas Sauer- 
stoff in das Schwefelteilgitter eingebaut werden kann, bevor das Gitter zu- 
sammenbricht und Eisenoxyd entsteht. 


Inhaltsübersicht 


1. Einleitung 


Das stöchiometrisch zusammengesetzte Eisensulfid FeS ist in der Lage, 
Schwefel im Überschuß aufzunehmen. Dabei treten Änderungen in den kri- 
stallographischen und magnetischen Eigenschaften auf, die von zahlreichen 
Autoren untersucht worden sind. Eine Zusammenfassung dieser Ergebnisse 
ergibt folgendes Bild. 

Nach Hägg und Sucksdorff!) und Haraldsen?) besitzen die Eisen- 
sulfide FeS,;, im Bereich 0,04 << n < 0,15 NiAs-Struktur. Oberhalb von 
n= 0,15 liegt ein Zweiphasengebiet, in dem neben FeS,,, Pyrit (FeS,) vor- 
handen ist. Unterhalb von n = 0,04 tritt eine ebenfalls hexagonale Über- 
struktur auf, die jedoch bei 120 °C in die NiAs-Struktur übergeht («-Umwand- 
lung). Aus der Abnahme der Dichte und des Zellvolumens mit steigendem 
Schwefelüberschuß schloß Haraldsen, daß bei n > 0 Leerstellen im Eisen- 
teilgitter vorhanden sein müssen (Mischkristalle vom Subtraktionstyp). 
Bertaut?) fand auf der Laue-Aufnahme eines natürlichen Pyrrhotin- 


1, G. Hasen. J. Sucksdorff, Z. physik. Chem. 22 B, 446 (1933). 
*) H. Haraldsen, Z. anorg. allg. Chem. 246, 169 (1941). 
3) F. Bertaut, Acta Cryst. 6, 557 (1953). 
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einkristalles der ungefähren Zusammensetzung Fe,S, eine Anzahl schwacher 
Überstrukturreflexe, die er einer geordneten Verteilung der Leerstellen 
zuschreibt. Aus der Lage der Reflexe schließt er, daß die Eisenatome so in 
übereinanderliegenden Ebenen angeordnet sind, daß nur jede zweite Ebene 
voll besetzt ist, während die dazwischenliegenden Ebenen die Leerstellen in 
geordneter Verteilung enthalten. 


Die magnetischen Eigenschaften der Eisensulfide sind ebenfalls recht 
verwickelt. Auf Grund der Untersuchungen von Haraldsen?)‘) an synthe- 
tisch hergestellten Proben kann man drei Bereiche unterscheiden: 


0<n< 0,07: Die Eisensulfide sind antiferromagnetisch mit einem 
Néel-Punkt von 325 °C, der unabhängig vom SchwefeliiberschuB ist. Unter- 
halb von n = 0,04 zeigt die Suszeptibilität einen Sprung bei 120 °C, der von 
der «-Umwandlung herrührt. Benoit®) fand, daß die Suszeptibilität des 
FeS oberhalb von 480 °C ein Curie-Weisssches Gesetz mit hohem, negativem 
6 befolgt. 

0,10 <n< 0,14: In diesem Bereich tritt eine spontane Magnetisierung 
auf. Der Curiepunkt steigt mit zunehmendem Schwefelgehalt von 300 °C 
auf 325 °C an. Gegenüber den bekannten Ferromagnetika zeigen sich aller- 
dings einige Besonderheiten: 


1. ist die Sättigungsmagnetisierung sehr niedrig, 

2. ist die Kristallanisotropie nach Messungen von Weiss®) an einem natür- 
lichen Einkristall so hoch, daß man in Richtung der hexagonalen Achse bis 
etwa 20 kOe nur einen schwachen linearen Anstieg der Magnetisierung mit 
dem Feld findet, und 

3. fanden Weiss und Forrer’) — ebenfalls an einem natürlichen Ein- 
kristall — ein flaches Maximum der spontanen Magnetisierung bei etwa 
— 160 °C. 

0,07 <n< 0,10: In dem Zwischengebiet fand Haraldsen bei etwa 
200 °C eine Spitze in der Suszeptibilität, die sich mit zunehmendem Schwefel- 
gehalt nach höheren Temperaturen hin verschiebt und dabei immer ausgeprägter 
wird. Bei FeS, „ erstreckt sie sich bereits über ein Temperaturintervall von 
etwa 70°. In diesem Bereich ist die Suszeptibilität stark feldstärkeabhängig. 
Haraldsen nimmt an, daß sich zunächst in einem schmalen Temperatur- 
intervall Ferromagnetismus ausbildet, der sich dann oberhalb von n = 0,10 
über das ganze Temperaturgebiet bis zum Curiepunkt erstreckt. Die Spitze 
ist auch von anderen Autoren beobachtet worden®)?). 

Erklärungen für das magnetische Verhalten der nichtstöchiometrischen 
Eisensulfide sind von Yosida!P) und Néel") gegeben worden. Beide be- 
trachten die Eisensulfide mit n > 0,10 als Ferrimagnetika. Yosida nimmt 
an, daß die beim FeS vorhandene Kompensation der beiden Untergitter durch 
den Einbau von überschüssigem Schwefel aufgehoben wird, indem die er- 


4) H. Haraldsen, Z. anorg. allg. Chem. 281, 78 (1937). 

5) R. Benoit, C. R. Acad. Sci. Paris 234, 2175 (1952). 

*) P. Weiss, J. Phys. (4) 4, 469 u. 829 (1905). 

7) P. Weiss u. R. Forrer, Ann. Chim. Physique (10) 12, 279 (1929). Ir 
8) R. Benoit, J. Chim. physique 52, 119 (1955). 

) T. Hirone, S. Maeda u. N. Tsuya, J. Phys. Soc. Japan 9, 736 (1954). 

10) K. Yosida, Progr. Theor. or 6, 356 ae 
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zeugten Fet++-Ionen nur Plätze in dem einen Untergitter besetzen. Néel 
folgert dagegen aus der Strukturbestimmung von Bertaut, daß die geordnete 
Verteilung der Leerstellen für die Unsymmetrie der beiden Untergitter ver- 
antwortlich ist. Wenn die antiferromagnetische Wechselwirkung zwischen den 
Eisenionen benachbarter Schichten größer ist als zwischen den Eisenionen 
innerhalb einer Schicht, so werden sich die Magnetisierungsvektoren auf- 
einanderfolgender Schichten antiparallel stellen. Dieser Sachverhalt ist, wie 
Lotgering!?) näher begründet, auch vom Standpunkt des Superaustauschs 
verständlich. Daß der Ferrimagnetismus nicht schon bei sehr kleinem Schwe- 
felüberschuß auftritt, rührt davon her, daß erst eine bestimmte Anzahl von 
Fet++-Ionen bzw. Leerstellen vorhanden sein müssen, damit eine geordnete 
Verteilung entstehen kann. 

Benoit®) fand, daß die Suszeptibilität von FeS,,, bei etwa 560 °C steil 
abfällt und oberhalb dieser Temperatur ein Curie-Weisssches Gesetz mit 
der gleichen Curiekonstanten befolgt wie das FeS. Er schreibt diese Ano- 
malie dem Übergang von der geordneten Verteilung der Leerstellen zur regel- 
losen Verteilung zu. Lotgering stellte jedoch fest, daß der von Benoit 
beobachtete Abfall der Suszeptibilität von geringen Verunreinigungen an 
Fe,0, herrührt. Die Übergang Ordnung—Unordnung erfolgt schon wenig 
oberhalb des Curiepunktes, da dort der von der Überstruktur der Leerstellen 
herrührende (0 0 1)-Reflex verschwindet. 

In der vorliegenden Arbeit sollten noch einmal die magnetischen Eigen- 
schaften der ferrimagnetischen Eisensulfide studiert und an Hand der ge- 
wonnenen Ergebnisse eine Deutung — insbesondere des Überganges Anti- 
ferromagnetismus—Ferrimagnetismus beim FeS,, — versucht werden. 
Außerdem interessierte das Verhalten der Eisensulfide gegenüber Sauerstoff 
und der Ersatz des Eisens durch andere Übergangsmetalle wie z.B. Cu, Ni 
und Co. 

Die Eisensulfide wurden durch 24stündiges Erhitzen der je nach Zusammen- 
setzung erforderlichen Mengen an Eisen und Schwefel in evakuierten Quarz- 
röhrchen bei 1000 °C hergestellt. Als Ausgangsmaterial dienten entkohltes 
Carbonyleisenpulver, das zunächst einige Stunden im Vakuum bei 400 °C 
entgast und anschließend bei der gleichen Temperatur in reinstem Wasser- 
stoff nachreduziert worden war, und Schwefel, der durch Lösen in Schwefel- 
kohlenstoff umkristallisiert und dann sublimiert worden war. 

Die Darstellung der Fe—Co-, Fe—Ni- und Fe—Cu-Sulfide erfolgte in der 
gleichen Weise, nur mußten die Mischsulfide nach dem ersten Glühen gut 
pulverisiert und dann nochmals 24 Stunden auf 1000 °C erhitzt werden, 
damit die Präparate homogen waren. Die erforderlichen Metallpulver wurden 
durch Reduktion ihrer Oxalate in reinstem H, erhalten. 


3. Meßergebnisse an den Eisensulfiden 


Es wurden Proben folgender Zusammensetzung untersucht: FeS, 99, 
FeS, 95; 99, usw. bis FeS,,; und FeS, 4. 


18) F, K. Lotgering, Philips Res. Rep. 11, 190 (1956). ES hon 
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Von jeder Probe wurden die Gitterkonstanten aus Debye-Scherrer- 
Aufnahmen bestimmt. Die Ergebnisse sind in Abb. 1 zusammen mit den 


f Werten von Haraldsen wiedergegeben. Dabei stellt die Länge eines Striches 
den durch die ungenaue 
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3.1 Gitterkonstanten 


4 Ausmessung der Linien- 
0% hler 
344 dar; systematische Feh- 
ler sind darin nicht ent- 
Br ° 342 halten. Wie man an den 

9° Ya cja-W erten erkennt, ist 
die relative Uberein- 

5.75 255 168 stimmung sehr gut. Die 

Tle Absolutwerte weichen 

755 dagegen um etwa 3%, 
° voneinander ab. Die 
164 Ursache dafür dürfte 


100 105 ı0 15 100 105 10 115 


wohl sein, daß in beiden 
Fällen mit einer 57,3- 
mm-Kamera gearbeitet 
wurde, deren Genauig- 
keit nicht größer ist. In 
Übereinstimmung mit 
den Angaben von Haraldsen treten beim FeS, ,, bereits Pyritlinien auf, 
und die Gitterkonstanten der hexagonalen Phase sind die gleichen wie beim 
FeS, Zwei weitere Probenreihen mit denselben Zusammensetzungen er- 


gaben innerhalb der Fehlergrenzen die gleichen Gitterkonstanten. 


th 
Die magnetischen Messungen wurden im wesentlichen mit einer Pendel- 
 waage ausgeführt, die im Prinzip der von Foéx und Forrer!?) beschriebenen 
gleicht. Die maximal erreichbaren Feldstärken betrugen für Messungen bei 
Zimmertemperatur etwa 18kOe und für Messungen bei hohen und tiefen 
Temperaturen etwa 7,5 kOe. 
Auf eingehendere Untersuchungen im antiferromagnetischen 
wurde verzichtet, 
berichtet haben. 


Schwefelgehalt 


_ Abb. 1. Gitterkonstanten der Eisensulfide FeS,;„. Auf der 
Abszisse ist 1 + n abgetragen. Striche: Eigene Messungen ; 
Kreise: Messungen von Haraldsen 


= ‘4 


3.2 Thermomagnetische Messungen Ko: 


Gebiet 
da hierüber bereits Haraldsen und Benoit ausführlich 
9, 19,6 2- 


Retox 


Beim Fe weichen die erheblic th ab, 
; da diese Zusammensetzung an der Grenze zwischen Antiferromagnetismus 
und Ferrimagnetismus liegt, so daB schon geringste Schwankungen in der 
ee Zusammensetzung die Werte völlig verändern. Die eigenen Messungen an 
zwei unabhängig voneinander hergestellten Proben ergaben ebenfalls ver- 
schiedene Suszeptibilitäten; bei der einen war sie feldunabhängig und betrug 


13) G. Foex u. R. Forrer, J. Physique Radium 7, 180 (1926). 
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4. — 19,2 - 10% cm? g!, bei der anderen war sie feldabhängig und betrug bei 
17 kOe x, = 50 - 10% cm? g 

Bei den ferrimagnetischen Proben tritt die schon erwähnte Schwierigkeit 
auf, daß man infolge der hohen Kristallanisotropie an polykristallinem Material 
bei den technisch erreichbaren Feldstärken keine Sättigung der Magnetisierung 
findet, sondern nur einen schwachen 


linearen Anstieg, der vom Verhalten in ” 

c-Richtung herrührt. Dies zeigt die ER 
Abb. 2 am FeS,,, bei Zimmertempe- 

ratur. Die spontane Magnetisierung 

kann man aus dieser Kurve auf fol- 

gende Weise abschätzen: 70 


Angenommen, das Material bestehe 
aus sehr vielen kleinen Einkristallen 
ohne gegenseitige Wechselwirkung. Bei 
einem Kristall soll die Normale auf 
die Basisebene mit dem Feld den Winkel 
90°-m bilden. Hat man bis zur Sätti- 0 5 70 75 20 
gung aufmagnetisiert und schaltet dann Feldstärke in k0e 
das Feld ab, so wird der Magneti- Abb. 2. Abhängigkeit der Magnetisierung 


sierungsvektor in die energetisch gün- von der Feldstärke bei FeS, 1, 
stige Basisebene zu liegen kommen, 


und zwar so, daß er mit dem vorher angelegten Feld den kleinstmöglichen 
Winkel — nämlich  — bildet. Die Remanenz o, beträgt dann o, - cosg. 


Spez. Magnetisierung in emE pro 9 


Mittelt man über alle Winkel p, so erhält man op = = 


7 Is 0,8 o,. Dieser 
Wert entspricht in Abb. 2 dem Schnittpunkt A 


der verlängerten Geraden mit der Achse. o, be- 20 
rechnet sich daraus zu 19,5emEg-!. Aus der 
Steigung der Geraden kann man die Suszep- Pr 


tibilität in c-Richtung bestimmen und daraus 
die Kristallanisotropie abschätzen. Man erhält: 


%e = 5-97,5-10-% em? g und daraus 


05° 6 3 
K = & 2- 10 erg 


“Xe 


if 
Uf 
| 


Zum Vergleich wurde K aus den Messungen 
von Pauthenet™), Weiss®) und Benoit) ab- 
geschätzt. Es ergeben sich für X bei Zimmer- 


Spez. Magnetisierung in emE pro g 


temperatur Werte zwischen 2 und 5- 10% erg em”®. . 10 172 1% 116 
Benoit schätzt ebenfalls o, und y, ab. Er Spk 
geht davon aus, daß bei völlig regelloser Ver- KR 


teilung ?/, der Kristallite mit der Basisebene Abb. 3. o,7500,385 in Abhängig- 
und !/, mit der c-Achse in in Feldrichtung liegen. a vers ow 
D: it hält — 93 4 E iat nd areı unabhangig voneinanc er 
ami = a er “0,4 emis g une hergestellte Probenreihen 

= 2,35 - 10% cm? 

Die höchste bisher an Einkristallen gemessene spontane Magnetisierung 
bei Zimmertemperatur beträgt 19,4 emE g—”). Dieser Wert stimmt mit der 
oben gemachten Abschätzung besser überein als mit der von Benoit. 

4) R. Pauthenet, C. R. Acad. Sci. Paris 284, 2261 (1952). 
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Abb. 3 zeigt die Abhangigkeit der Magnetisierung bei 17 kOe vom Schwe- 
felgehalt. Es wurden drei unabhängig voneinander hergestellte Probenreihen 
untersucht. Die Übereinstimmung der Meßergebnisse an den verschiedenen 
Proben gleicher Zusammensetzung ist recht gut. Die Magnetisierung steigt 
zwischen FeS, 4 und FeS,,, sehr steil an und hat ein Maximum bei etwa 
FeS, 33- 

Die Temperaturang yingigkeit der Magnetisierung in einem äußeren Feld 
von 7,5 kOe ist in Abb. 4 für die Zusammensetzungen FeS, 45 bis FeS, ,, dar- 
gestellt. Alle Messungen bei höheren Temperaturen wurden unter Vakuum 
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Abb. 4. Temperaturabhängigkeit der Magnetisierung der ferrimagnetischen Eisensulfide 
bei H = 7,5 kOe 


vorgenommen, um eine Oxydation der Proben weitgehendst zu vermeiden. 
Trotzdem waren die beim Aufheizen und Abkühlen erhaltenen Kurven auch 
bei langsamer Temperaturänderung nur auf etwa 5%, reproduzierbar. Wieder- 
gegeben sind in Abb. 4 nur die beim Aufheizen erhaltenen Kurven. Einzelne 
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Meßpunkte sind nicht eingetragen, da etwa alle 5° gemessen wurde, woraus 
sich praktisch eine zusammenhängende Kurve ergibt. FeS, |, zeigt die gleiche 
eigentümliche Temperaturabhängigkeit, die schon Haraldsen beobachtet 
hatte. Zwischen 475 und 550 °K tritt ein schmaler ferrimagnetischer Bereich 
als Übergang vom Antiferromagnetismus zum Ferrimagnetismus auf. Bei 
den schwefelreicheren Proben erstreckt sich dann der Ferrimagnetismus vom 
absoluten Nullpunkt bis zum Curiepunkt, der sich nur wenig mit der Zu- 
sammensetzung ändert. Er beträgt bei FeS,,, 590 °K, bei FeS,,, 600 °K 
und bleibt von FeS, „, an konstant bei 610 °K. Sämtliche Kurven der ferri- 
magnetischen Sulfide lassen ein flaches Maximum der Magnetisierung bei 
etwa 150 °K erkennen. Leider stand kein flüssiger Wasserstoff zur Verfügung, 
um den Temperaturgang nach noch tieferen Temperaturen hin zu verfolgen. 
Die Beobachtungen sind jedoch in bester Übereinstimmung mit denen von 
Weiss und Forrer an einem natiirlicchen Einkristall und mit der von Benoit 
an FeS, j95 gemessenen Kurve. 

Die Tatsache, daß der Curiepunkt von FeS, ,; an konstant bleibt und 
die Magnetisierung bei der gleichen Zusammensetzung ein Maximum hat, 
läßt vermuten, daß schon hier das Zweiphasengebiet Pyrrhotin—Pyrit be- 
ginnt. Mit Hilfe von Debye-Scherrer-Aufnahmen ist nicht genau zu ent- 
scheiden, bei welcher Zusammensetzung das Zweiphasengebiet beginnt, da 
einmal die Linien des Pyrit erst von einer gewissen Konzentration an auf- 
treten, und zum andern nicht genau zu erkennen ist, von wo ab die Dimen- 
sionen der Einheitszelle sich nicht mehr ändern, da die Gitterkonstanten mit 
einem gewissen Fehler behaftet sind. 


4. Deutung der magnetischen Eigenschaften der Eisensulfide 


Im folgenden soll versucht werden, die recht komplizierten Eigenschaften 
der nichtstöchiometrischen Eisensulfide durch ein einfaches Modell qualitativ 
zu erklären, wobei der Übergang vom Antiferromagnetismus zum Ferrima- 
gnetismus und die anomale Temperaturabhängigkeit der spontanen Magneti- 
sierung von besonderem Interesse sind. Wir werden dabei ausgehen von der 
Verbindung Fe,S, und später die Abhängigkeit der spontanen Magnetisierung 
vom Schwefelüberschuß und der Temperatur betrachten. he 

4.1 Die Sättigungsmagnetisierung der Eisensulfide wily 


Sowohl die eigenen Messungen als auch die von Benoit haben ergeben, 
daß die größte Magnetisierung bei der Zusammensetzung FeS, , liegt. Dieser 
Wert kommt der von Néel auf Grund der Strukturbestimmung von Bertaut 
vorgeschlagenen Formel Fe,S, sehr nahe. Die Sättigungsmagnetisierung 
dieser Verbindung läßt sich experimentell nur näherungsweise bestimmen, 
da zwei Schwierigkeiten auftreten: 

1. ist die spontane Magnetisierung an polykristallinem Material nicht meB- 
bar, da die Magnetisierung wegen der hohen Kristallanisotropie sogar zwischen 
10 und 20 kOe noch linear mit dem Feld ansteigt, und 

2. ist die übliche Extrapolation von o,, „ auf T = 0 nicht durchführbar, 
da die Magnetisierung bei tiefen Temperaturen ein Maximum überschreitet. 

Die erste Schwierigkeit hatten wir durch eine Abschätzung umgangen, 
die auf der statistischen Verteilung der einzelnen Kristallite beruhte und die 
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für Goo, 993 den Wert 19,8 emE g™ ergab. Für den Coens VON Ooo, a93 auf 
00,0 kommt uns zugute, daß das Maximum bei etwa 150 °K ziemlich flach 
ist, so daß man 0%, 9 Goo, 993 20 emE setzen kann. Wir werden daher 
den folgenden theoretischen Betrachtungen für das Fe,S, einen experimentellen 
Wert von etwa 20 emE g~ gegenüberstellen. 

Die Berechnung der Sättigungsmagnetisierung von Fe,S, geschieht in der 
üblichen Weise, indem das Kationenteilgitter in zwei Untergitter aufgeteilt 
wird, die in entgegengesetzter Richtung magnetisiert sind. Ausgehend von 
der Formel Fe}+Fez++LS;— (L bedeutet eine Leerstelle) kommt man zu 
einer Reihe von Grenzfällen, je nachdem wie die Fett, die Fet++ und die L 
auf die einzelnen Untergitter verteilt sind. Dabei soll ein Fe++-Ion ein magneti- 
sches Moment von 4 u, und ein Fe+++-Ion ein magnetisches Moment von 
5a, haben. Fünf solcher Grenzfälle, sowie die dafür berechneten Sätti- 
gungsmagnetisierungen sind in der Tab. 1 wiedergegeben. 


Tabelle 1 


Modell | Untergitter 1 | Untergitter 2 Tog 
1 2,5 Fett 1¥Fet++ 0,5 L 2,5 Fet+ 1¥Fet++ 0,5 L 0 
2 1,5 Fett 2Fett+ 0,5 L 3,5 Fett 0,5 L 17,5 emkg-! 
hing 2 Fett 1Fet+ 1 L 3 Fett 1Fet++ | 35 emEg! 
ase 1 Fett 2Fet+ 1 L 4 Fet+ | 17,5 emEg 
van | 2 Fett 2¥Fett+ 3 Fett 1 L | 52,5 emEg* 


Von den theoretisch berechneten Sättigungsmagnetisierungen stimmen nur 
die der Modelle 2 (Yosida) und 4 mit dem experimentell gefundenen Wert 
überein. Eine Entscheidung zwischen diesen beiden Möglichkeiten ist 
auf Grund der folgenden Betrachtungen möglich: 

1. Aus der Tatsache, daß im ferrimagnetischen Gebiet eine Überstruktur 
der Leerstellen auftritt, ist zu entnehmen, daß die Leerstellen nur auf einem 
Untergitter sein können. Diese Annahme steckt aber nur im Modell 4. Daß 
die Fe*++*-Ionen in dem gleichen Untergitter sitzen, erscheint verständlich, 
da sie sich bevorzugt in der Nähe der Leerstellen aufhalten werden, um La- 
dungsansammlungen zu vermeiden. 

2. Betrachten wir das Verhalten oberhalb von T = 0: Die Temperatur- 
bewegung wird bestrebt sein, einen bei 7 = 0 vorhandenen Ordnungszustand 
zu zerstören, d. h. Gleichverteilung herzustellen. Bei Modell 2 würden Elek- 
tronenplatzwechsel zwischen den Fe++- und den Fet++-Ionen stattfinden, 
die schließlich eine Gleichverteilung der Fe+++-Ionen zur Folge haben werden. 
Dabei nimmt die spontane Magnetisierung mit wachsender Temperatur ab. 
Bei Modell 4 ist der Vorgang etwas komplizierter; da hier sowohl Elektronen- 
platzwechsel als auch Ionenplatzwechsel stattfinden. Von den möglichen 
Platzwechselvorgängen werden die Elektronenplatzwechsel wegen ihrer 
kleineren Aktivierungsenergie bei niedrigeren Temperaturen einsetzen als 
die Ionenplatzwechsel, so daß zunächst eine Gleichverteilung der Fe++-Ionen 
angestrebt wird, bevor mit Hilfe der Ionenplatzwechsel die völlige Gleich- 
verteilung hergestellt wird. Das würde bedeuten, daß wir uns bei tiefen Tem- 
peraturen zunächst in Richtung auf das Modell 3 bewegen, bevor der Übergang 
zu Modell 1 einsetzt. Die spontane Magnetisierung müßte also bei tiefen Tem- 
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peraturen zunächst ansteigen und erst bei höheren Temperaturen wieder ab- 
fallen. Das ist aber genau das Verhalten, das Weiss und Forrer und Benoit 
gefunden hatten, und das wir bei den o (T)-Kurven der Abb. 4 bestätigen 
konnten. 


Da nach den bisherigen Betrachtungen das magnetische Verhalten von 
Fe,S, am besten durch ein Modell beschrieben werden kann, in dem sowohl die 
Fe++*-Ionen als auch die Leerstellen nur Plätze in dem einen Untergitter 
besetzen, soll diese Anordnung auch den folgenden Überlegungen zugrunde 
gelegt werden. 


Um die Abhängigkeit der Sättigungsmagnetisierung vom Leerstellen- 
gehalt zu berechnen, ersetzen wir, da es sich ja um Mischkristalle vom Sub- 
traktionstyp handelt, die bisher gebräuchliche Formel FeS,;, durch Fe,_,S, 
wobei « = n/1-+ n ist. Berücksichtigen wir noch die Zahl der Leerstellen 
und der Fe***-Ionen, so lautet die genaue Zusammensetzung Fej+, „Fejt+L,S. 
Die Sättigungsmagnetisierung berechnet sich unter Zugrundelegung des obigen 
Modells zu: 


fd, 


Dabei bedeuten N, die ‘Dicipasaiail Zahl und M das Molekular- 
gewicht. Da sich das een, icht nur wenig mit der Zusammensetzung 
ändert, müßte die Sättigungsmagnetisierung nach Gl. (1) nahezu linear mit 
zansteigen. Eine Überprüfung dieses Sachverhalts ist leider nicht möglich, da 
wir das magnetische Verhalten der nichtstöchiometrischen Eisensulfide nur 
oberhalb des Siedepunktes der flüssigen Luft kennen. Bei höheren Tem- 
peraturen tritt der Ferrimagnetismus jedenfalls erst von x ~ 0,09 ab auf. 


4.2 Berechnung von oo,, mit Hilfe der Theorie von Bragg-Williams 


Die Te ‚mperaturabhängigkeit der spontanen Magnetisierung der Eisen- 
sulfide wird durch zwei Dinge hervorgerufen: Erstens dadurch, daß die Tem- 
peraturbewegung die völlige Ausrichtung der Elementarmagnete zerstört, 
und zweitens dadurch, daß die geordnete Verteilung der Leerstellen und der 
Fe*++-Ionen zunichte gemacht wird. Der erste Vorgang wird allgemein 
durch die Theorie von Langevin oder deren quantenmechanische Erweiterung 
beschrieben. Den zweiten Prozeß kann man durch die Theorie von Bragg 
und Williams!) beschreiben. Diese Theorie wurde zur Beschreibung des 
Überganges Ordnung —Unordnung bei Legierungen mit zwei Komponenten 
entwickelt. Sie stellt zwar nur eine sehr grobe Näherung dar, ist jedoch in 
der Lage, die wesentlichsten Dinge qualitativ recht gut zu beschreiben. Um 
diese Theorie auf die Eisensulfide anwenden zu können, ist es notwendig, sie 
auf ein System mit drei Komponenten zu erweitern, wobei der Umfang an 
Rechnung so stark zunimmt, daß auf eine ausführliche Ableitung aller Glei- 
chungen verzichtet werden mußte. 


®) W.L. Bragg u. E. J. Williams, Proc. Roy. Soc. London (A) 145, 699 (1934); 
(A) 151, 540 (1935). 
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Der Gang der Rechnung ist folgender: 

Zunächst wird für eine beliebige Verteilung der beiden Ionenarten und der 
Leerstellen auf die beiden Untergitter die freie Energie F, bestehend aus der 
Anordnungsenergie E und dem Entropieterm —T'S, für eine beliebige Tem. 
peratur berechnet. Durch Aufsuchen des Minimums von F in bezug auf die 
Konzentration der Fe***.Ionen und der Leerstellen erhält man den Gleich- 
gewichtszustand bei dieser Temperatur T. Auf diese Weise läßt sich für jede 
Temperatur die Verteilung der Teilchensorten auf die beiden Untergitter be. 
rechnen. Ist diese bekannt, so kann man daraus die Abhängigkeit der spon- 
tanen Magnetisierung von der Temperatur ableiten. Die zweiten Ableitungen 
von F liefern dann Aussagen über den A-Punkt (Umwandlungspunkt Ord- 
nung—Unordnung). 


4.2.1 Der Übergang von der geordneten zur ungeordneten Ver. 
teilung 

Unser Modell liefert für die Zusammensetzung Fej+,,Fejj+L,S— bei 

T = 0 folgende Verteilung der verschiedenen Teilchensorten auf die beiden 


Untergitter: 
Tabelle 2 


| Untergitter 1 | Untergitter 2 | Summe 
L Nu =Nx No = 0 N, Nx 
Fett N,, = N (4. — 32) Nos = N/2 N, = N (1 —32) 
Fett: | Ns = N. =0 N; =2Ner 
Summe | WN/2 N/2 N 


Bei der Temperatur 7 soll bereits ein Teil der Fe++*-Ionen und der Leer- 
stellen in das Untergitter 2 übergegangen sein, so daß auf dem Untergitter 1 
nur noch Nxr, Leerstellen und 2 Nar, Fet++-Ionen übriggeblieben sind. Damit 
variieren r, und r, zwischen 1 und !/,; und zwar bedeuten ry = rz = 1 völlige 
Ordnung und ry = r; = !/, völlige Unordnung. Die Temperaturabhängigkeit 
von r, und r, wird natürlich verschieden sein. Die Verteilung bei der Tem- 
peratur 7 sieht also folgendermaßen aus: 


Tabelle 3 

| Untergitter 1 | Untergitter 2 | Summe 

Fett Ny Ny +27 | N, = N (1 —32) 
+22 rs) | 

'ettt | N= 2Ne8, | N,=2Nx 
L Ny = Nev, | Nog = N x (1 — 1%) 
IN | N 
Summe | N 


Die Anordnungsenergie ist nach Bragg und Williams gegeben durch 
den Ausdruck: 


A N N N 
= = —No%—N2%2— N3 93 — (7- Ns) top -(5 Ny) Wog -(5 No ) ss 
j 2 r r r r r 
N Nog + No Na — Ng, N35) wog + (Nor Noo — No Noe + Nyy N32) 
+ Noy Nog + Nar Noe + Nor N) 
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Die Bedeutung der N, und der N,, ist aus Tab. 3 ersichtlich. Die, - 
die w 4B haben folgende Bedeutung: er 
2y,/2 (<= Zahl der nächsten Nachbarn) ist die Wechselwirkungs- 
zweier benachbarter A-Atome und — (9, + 9, + w, „)/2 ist die Wechsel- 
wirkungsenergie eines A-Teilchens mit einem benachbarten B-Teilchen. 

Diesem elektrostatischen Anteil der Anordnungsenergie ist noch ein zweiter, 
magnetischer Anteil hinzuzufügen, der berücksichtigt, daß beim Übergang 
anes Teilchens mit einem magnetischen Moment in das andere Untergitter 
auch dessen Spinrichtung umgekehrt werden muß. Aus der Molekularfeld- 


theorie erhält man für diesen Term: RE 
Ey = — «I,J, =a 


wobei x die Molekularfeldkonstante ist, und J, und J, die REIN 
der Untergitter 1 und 2 bedeuten. Drückt man J, und J, durch die in Tab. 3 
angegebenen Konzentrationen der Ionen aus und setzt für das magnetische 
Moment des Fe*++-Ions 4 u, und für das Fett+ 5 u,, so erhält man aus 
Gl. (3): 


Ey = (T) + 5 ny (TY) (7) 


32° (T)). (3a) 
y, (1) ist die relative spontane Magnetisierung von Teilchen mit dem Ge- 
samtdrehimpuls J, die sich durch Auflösen der Brillouin-Funktion ergibt. ue re: 
Vernachlässigt man bei 7, (7) den Unterschied zwischen J = 4/2 und 5/2, 
so ergibt sich aus Gl. (3a): 


er: 
Eu = (16 Noy Nog + 20 Noy + 20 Nog + 25 Nay (3b) 


Die Entropie der in Tab. 3 angegebenen Verteilung beträgt: Re 5 

Nor! Nos! N a'Ng! 1 Na! N52! 

Für die Rechnung wird in Gl. 4 In x! durch «In x — x ersetzt. Die freie 

Energie F erhält man aus den Gln. (2), (3b) und (4), wobei noch die N, und die 

N, durch die in Tab. 3 angegebenen Konzentrationen ausgedrückt werden 

müssen. Setzt man die beiden Ableitungen von F nach r, und rz Null, so er- _ R 

hält man die folgenden beiden Bestimmungsgleichungen für r, und r;: = ; 


(ryt 1) + 40 (2 — 1) + 


= —kIn W =~ kin( 


05-zn -2ın 
Nuz 
— ER) (2r——1) +22 +8n—5) — 22 (dr, — 275 


I-n 05-ın-2er 


Das Minimum liegt für 7 = 0 bei ry = rg = Lund wandert mit zunehmender 
Temperatur gegen = 13 1/9. Von einer bestimmten Temperatur 7, 
ab bleibt es dann bei ry = 13 = 1/,. Zwischen 7’, einerseits und den wy, 
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und « eniieetealiia erhält man aus den zweiten Ableitungen von F nach 7 
und r, zwei Beziehungen: 


kT, 


Bei bekanntem 7’, sind zwei der vier Unbekannten (Wgg, Wo3, Wo3, &) bestimm- 
bar, nämlich wo, und etwa w,. « kann aus dem Curie-Weissschen Gesetz 
bei hohen Temperaturen abgeschätzt werden, während man über weg nur aus- 
sagen kann, daß es etwa von der gleichen Größe sein wird wie Wyp. 


4.2.2 Der Verlauf der spontanen Magnetisierung 


Sind ry (T) und r, (T) bekannt, so läßt sich die spontane Magnetisierung aus 
nachstehender Gleichung berechnen: 


Oo, T (5 —axr,(T)—22 ry (7) 4 ny (7) + 2x7,(T)- 5 up (T) 


sey Ubereinstimmung mit Gl. (1) erhalten wir fiir T = 0: 


Durch Einsetzen von GL. (8) in Gl. (7b) ergibt sich fiir die relative spontane 
Magnetisierung: 


-5 4.2.3 Auswertung der Gln. (5) und (9) 


Die Berechnung von r,(T) und r,(T) aus dem Gleichungssystem (8) 
geschah numerisch. Aus den so erhaltenen Werten wurde dann die relative 
spontane Magnetisierung mit Hilfe von Gl. (9) in Abhängigkeit von der Tem- 
peratur bestimmt. Zunächst soll der Fall x = 1/, (Fe,S,) betrachtet werden. 
Zur Diskussion der Ergebnisse wurden einige charakteristische Fälle heraus- 
gegriffen, die aber zu einem allgemeinen Überblick führen. 


Die berechneten Kurven sind in Abb. 5a bis 5f dargestellt, und zwar für 
die beiden Fälle: Magnetische Wechselwirkung Null (a, b, c) und starke anti- 
ferromagnetische Wechselwirkung (d, e, f). 

Wenn keine magnetische Wechselwirkung vorhanden ist, wird 1 (7) 
durch die Größe von wg so gut wie gar nicht beeinflußt, wogegen der starke 
Abfall von r, (7) mit abnehmendem w,, nach tieferen Temperaturen hin ver 
schoben wird. Dies hat zur Folge, daß bei kleinerem wo, die relative spontane 
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Magnetisierung im mittlerem Temperaturbereich ein höheres Maximum 
is bei größerem 

Eine starke antiferromagnetische Wechselwirkung beschleunigt den Ab- — 
fall von ry (7) und verlangsamt den Abfall von r, (7) unterhalb des Curie- _ 


T 


Abb. 5. Berechnete Kurven für TA > 0. Linke Reihe (a, b, c): Ohne antiferromagnetische ‚Re 
Wechselwirkung. Die an den Kurven angegebenen Parameter sind das Verhältnis wgg/Wg- 
a: Temperaturabhängigkeit von r,, 6: Temperaturabhängigkeit von r;, undc:Temperatur- 
abhängigkeit der relativen spontanen Magretisierung. Rechte Reihe (d,e, f): Starke | 
antiferromagnetische Wechselwirkung. Die Bedeutung der Kurven ist die gleiche wie bei er 

(a, b, c) 


sierung in dem mittleren Temperaturbereich erheblich verringert. Nimmt die 
antiferromagnetische Wechselwirkung noch weiter zu, so verringert ich 
schließlich auch die Sättigungsmagnetisierung 0%, , und im Grenzfall unend- _ 
lich hoher Wechselwirkung ist nur noch Antiferromagnetismus vorhanden. 
2. 
Die berechneten Kurven zeigt die Abb. 6. Der wesentliche Unterschied Be 
gegenüber dem Fall 7, > 6 ist, daß mit zunehmender antiferromagnetischer RE 
Wechselwirkung das Verschwinden des Ferrimagnetismus nach tieferen Tem- th 
Peraturen hin verschoben wird. Die Ordnungsenergie ist in diesem Falle nur 
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bei tiefen Temperaturen groß genug, um den ferrimagnetischen Zustand auf. 
rechtzuerhalten; bei höheren Temperaturen überwiegt dagegen die anti. 
ferromagnetische Wechselwirkungsenergie, so daß Antiferromagnetismus 
auftritt, der dann bei 7’ = 6 in Paramagnetismus übergeht. 


Wird 7, kleiner als 0, so nimmt auch die Ordnungsenergie weiter ab, 
Bei Vorhandensein einer. antiferromagnetischen Wechselwirkung bedeutet 
das aber, daß der oben beschriebene Übergang Ferrimagnetismus— Anti. 
ferromagnetismus nach noch tieferen Temperaturen hin verschoben wird. 
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Abb. 6. Berechnete Kurven für 74 — 0. Die Bedeutung der Abbildungen a bis f ist di 


gleiche wie in Abb. 5 


Wird x kleiner als !/,, d.h. nimmt die Zahl der Leerstellen ab, so 
auch die Ordnungsenergie kleiner. Das bedeutet aber, daß der Einfluß de 
magnetischen Wechselwirkung stärker wird und die spontane Magnetisierun 
verringert, eventuell sogar zu Null macht. Die Änderung der thermomagneti 
schen Kurven mit abnehmendem x vollzieht sich in ähnlicher Weise wie der 
im vorigen Absatz diskutierte Übergang von Tı>® zu Ta< 6. En 
besonderer Fall, der eintreten kann, ist in Abb. 7 dargestellt. Die antiferre 
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{| magnetische Wechselwirkung ist so stark, daß die spontane Magnetisierung 
if bereits bei tiefen Temperaturen auf Null abnimmt. Da aber T, > ist, 
sf überwiegt kurz vor dem Curiepunkt die Ordnungsenergie wieder, so daß die 
spontane Magnetisierung noch ein kleines 
Maximum zeigt. Eine Andeutung dieses 

Verhaltens ist schon in Abb. Sfenthalten, 
: wo die Ordnungsenergie noch groB genug 
ist, um die spontane Magnetisierung nicht 
ganz verschwinden zu lassen. Bei sehr Sg 
kleinem x ist dann die Ordnungsenergieso I 
gering, daß man im wesentlichen Antiferro- 05 
magnetismus findet. 


4.2.4 Diskussion der Ergebnisse ee 0 - 
Die Kurve in Abb. 7 gibt bei höheren . te 
Temperaturen annähernd das Verhalten Te ‘Fen 
des FeS, wieder. Der ferrimagnetische 


Bereich bei sehr tiefen Temperaturen ist Abb. 7. Temperaturabhängigkeit der 


: 2 2 relativen spontanen Magnetisierung 
bisher nicht beobachtet worden. Die für den Fall geringer Ordnungsener- 


Curietemperatur beträgt beim FeS, 4) gie und starker antiferromagnetischer 
etwa 275 °C. T, müßte dann theoretisch Wechselwirkung bei Tı > 0 
etwas oberhalb von 275 °C liegen. Beim 

Fe,S, (= FeS, 34) ist 6 = 330 °C und 7, = 320—340 °C?2). Das bedeutet, daß 
sich 7’, mit x nur wenig ändern kann. Schätzt man aus dem Curie-Weiss- 
schen Gesetz die Molekularfeldkonstante x ab, so erhält man für die 
antiferromagnetische Wechselwirkung einen Wert, der zwischen den beiden 
hier diskutierten Fällen liegt. Ein Vergleich zwischen dem Curie-Weiss- 
schen Gesetz für FeS,o, und Fe,S, zeigt, daß sich auch x nur wenig 
mit x ändert). Über w,, kann man nichts aussagen; es ist aber anzunehmen, 
daß es ungefähr gleich dem wp, sein wird. Für den Verlauf der spontanen 
Magnetisierung erwartet man daher auf Grund der Abb. 5, 6 und 7, daß 
der wesentliche Abfall der spontanen Magnetisierung erst nahe am Curie- 
punkt stattfindet, während sich in dem restlichen Temperaturgebiet die 
spontane Magnetisierung entweder nur wenig ändert oder sogar ein Maxi- 
mum überschreitet. Dieses qualitative Ergebnis stimmt auch mit dem wirk- 
lichen Verhalten einigermaßen überein. Mit abnehmendem x wird dann der 
Einfluß der antiferromagnetischen Wechselwirkung stärker, wodurch die 
spontane Magnetisierung im mittleren Temperaturbereich abfällt. Als Bei- 
‚di spiel dafür sei die Kurve für w3/wo. = 0,7 in Abb. 5f genannt, die im Aussehen 
mit der von Haraldsen am FeS,,, gemessenen Kurve vergleichbar ist. 
Nimmt x noch weiter ab, so tritt der oben genannte Fall ein, daß in einem mitt- 
leren Temperaturbereich Antiferromagnetismus vorhanden ist und kurz vor 
dem Curie punkt ein schmaler ferrimagnetischer Bereich auftritt (s. a. Abb.4). 
Bei noch weiterer Abnahme von x tritt dann im wesentlichen nur noch Anti- 
ferromagnetismus auf. 


Man sieht damit, daß die hier benutzte Theorie die magnetischen Eigen- 
schaften der nichtstöchiometrischen Eisensulfide qualitativ recht gut wieder- 
gibt. Eine quantitative Übereinstimmung ist nicht zu erwarten, da einmal die 
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Bragg-Williamssche Theorie nur eine erste Näherung darstellt, und zum 
anderen der hier gemachte Ansatz noch Vereinfachungen enthält (z.B, 
wurde der magnetische Anteil der Entropie nicht berücksichtigt). 

Die Berücksichtigung von Platzwechselvorgängen in der vorliegenden 
Theorie wird auch durch die Tatsache gerechtfertigt, daß die Geschwindig- 
keit von Temperaturänderungen einen Einfluß auf die Magnetisierung bzw. 
die Suszeptibilität hat. Hirone und Mitarbeiter?) haben an FeS, 4 x (7) bei 
drei verschiedenen Heiz- und Abkühlgeschwindigkeiten gemessen. Sie fanden, 
daß bei schnellem Durchlaufen der Temperaturkurve (300°/h und 100°/h) das 
hohe Maximum der Suszeptibilität beim Abkühlen bei tieferen Temperaturen lag 
als beim Aufheizen. Bei 30°/h war dagegen die y (T)-Kurve völlig reprodu- 
zierbar. Bjelow!®) hat ebenfalls am FeS, 1) die zeitliche Änderung der Sus. 
zeptibilität direkt gemessen, indem er die Probe von hohen Temperaturen her 
schnell auf eine Temperatur im Gebiet des Maximums der Suszeptibilität 
abgekühlt hat und dann bei konstanter Temperatur den zeitlichen Anstieg 
der Suszeptibilität verfolgte. Aus beiden Untersuchungen geht hervor, daß 
die Platzwechselvorgänge ziemlich schnell ablaufen. Bei den eigenen Mes- 
sungen stimmten die Kurven‘beim Aufheizen und Abkühlen ebenfalls um so 
besser überein, je langsamer sie durchlaufen wurden; allerdings waren sie auch 
bei sehr langsamen Temperaturänderungen nicht völlig reproduzierbar. Wegen 
des schnellen Ablaufs der Platzwechselvorgänge gelang es nicht, durch Ab- 
schrecken auf Zimmertemperatur den Zustand bei hohen Temperaturen ,,ein- 
zufrieren‘‘. Bei diesen Versuchen ergab sich außerdem, daß sich beim Ab- 
schrecken von Temperaturen oberhalb 400 °C ein Teil des Schwefels an der 
Wand des Quarzröhrchens, in dem die Probe eingeschmolzen war, nieder- 
schlug. Vermutlich ist der Partialdruck des Schwefels bei diesen Temperaturen 
schon erheblich. 
5. Leitfähigkeitsmessungen 

Die im vorigen Abschnitt diskutierten Platzwechselvorgänge in den Eisen- 
sulfiden dürften sich nicht nur in den magnetischen Eigenschaften bemerkbar 

40 machen, sondern müßten auch in den 
elektrischen Eigenschaften eine Rolle 
spielen. Daher wurde die Leitfähigkeit als 
Funktion des Schwefelüberschusses und 
der Temperatur untersucht. 


Fiir die Messungen standen nur Pulver 
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Spez Leitfähigkeit 


20 zur Verfügung. Die feinpulverisierten 
Proben wurden in Glasröhrchen gefüllt 

> und mit Hilfe zweier, genau passender 

10 Kupferstempel, die zugleich als Kontakte 
? dienten, zusammengepreBt. Den Druck 


erzeugte ein auf dem verlängerten oberen 
0 Kupferstempel lastendes Gewicht. Die 
100 105 170 15 erreichte Dichte betrug etwa ?/, der rönt- 


Schwefelgehalt genographischen Dichte. Die gemessenen 
Abb. 8. Spezifische Leitfähigkeit dr 
Eisensulfide FeS,,, in Abhängigkeit 16) K. P. Bjelow, A. W. Saljesski u 


vom Schwefelgehalt. Auf der Abszisse A.S. Lwowa, Dokl. Akad. Nauk SSSR Krist. 
ist 1 + n abgetragen 1, 198 (1956), 
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Leitfähigkeiten waren auf etwa 25%, reproduzierbar. Die Stromspannungs- 
kurven erfüllten über drei Zehnerpotenzen das Ohmsche Gesetz innerhalb 


der Meßgenauigkeit von etwa 1%. 


Die Ergebnisse der Leitfähigkeitsmessungen sind in Abb. 8 und 9 dar- 
gestellt. Abb. 8 zeigt die Leitfähigkeit in Abhängigkeit vom Schwefelüber- 
schuß bei Zimmertemperatur. Sie liegt bei allen Proben ziemlich hoch — 
nur etwa drei Zehnerpotenzen unter der der Metalle — und nimmt mit wach- 
sendem Schwefelüberschuß zu. Der beobachtete Anstieg der Leitfähigkeit 
mit dem Schwefelüberschuß ist durchaus verträglich mit der Annahme von 


Platzwechselvorgängen zwischen den Fet+-, 
Fe+++-Ionen und den Leerstellen, da ja mit 
zunehmender Abweichung von der stöchio- 
metrischen Zusammensetzung die Zahl der 
Fett*-Ionen und der Leerstellen ansteigt, 
und sich damit die Möglichkeit für Platz- 
wechsel erhöhen muß. Von den möglichen 
Platzwechseln dürften bei Zimmertempe- 
ratur die Elektronenplatzwechsel den größten 
Anteil liefern, da die Aktivierungsenergie 
für Ionenplatzwechsel wesentlich größer ist 
als die für Elektronenplatzwechsel (s. a. die 
Kurven für r,(T) und r,(T) in Abb. 5 
und 6). Beim stöchiometrischen Sulfid ent- 
fallen die oben diskutierten Vorgänge. Daß 
dessen Leitfähigkeit trotzdem noch erheblich 
ist, mag einmal daran liegen, daß es unmög- 
lich ist, eine Probe der exakten stöchio- 
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Abb. 9. Spezifische Leitfähigkeit 
der Eisensulfide in Abhängigkeit 
von der Temperatur 


metrischen Zusammensetzung herzustellen, ganz abgesehen von den Ver- 
unreinigungen, die in den Ausgangsmaterialien enthalten sind. Zum anderen 
spricht der metallische Charakter aller Verbindungen mit NiAs-Struktur 
dafür, daß die stöchiometrischen Verbindungen noch eine eigene Leitfähigkeit 


besitzen. 


Die Temperaturabhängigkeit der Leitfähigkeit einiger Eisensulfide ist 
in Abb. 9 dargestellt. In dem untersuchten Temperaturbereich von 100 bis 
300 °K fällt die Leitfähigkeit monoton mit abnehmender Temperatur. Ein 
Minimum bei 160 °K, wie es Benoit’) gefunden hatte, konnte nicht fest- 


gestellt werden. 


Eine genaue Analyse des elektrischen Verhaltens der Eisensulfide erscheint 
gegenwärtig wenig sinnvoll. Dazu müßten eingehendere Untersuchungen an 
größeren Kristallen oder besser Einkristallen gemacht werden, wozu sicher 
auch andere Messungen wie z.B. der Hall-Effekt herangezogen werden 
müßten. Es sollte hier nur geprüft werden, ob sich aus dem elektrischen Ver- 
halten ernsthafte Einwände gegen die im vorigen Kapitel gemachte Annahme 
von Platzwechseln zwischen den Fet+-, Fet++-Ionen und Leerstellen er- 


geben. Das scheint aber nicht der Fall zu sein. 


6. Einfluß von Sauerstoff auf die Eisensulfide 


Lotgering!?) und Bjelow!®) haben gefunden, daß die Eisensulfide unter 
der Einwirkung von Sauerstoff teilweise in Magnetit umgewandelt werden 
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und re so dem von Benoit?) bei 550 °C gefundenen steilen. Abfall der 
Suszeptibiltität als Curiepunkt des gebildeten Magnetits. Die eigenen Unter- 
suchungen deuten darauf hin, daß bei vorsichtiger Oxydation der Sauerstoff 
zunächst in das Schwefelteilgitter eingebaut wird. Erst bei weiterer Oxydation 
bilden sich dann Eisenoxyd. - 
Zu dieser Annahme führte die folgende Beobachtung: 
Eine Probe der Zusammensetzung FeS, j9, die bei Zimmertemperatur eine 
feldstärkenunabhängige Suszeptibilität von 19,2 - 10% cm?/g zeigte, wurde 
an Luft mehrmals auf etwa 400 °C aufgeheizt und abgekühlt und dabei der 
Verlauf der Suszeptibilität verfolgt 


1 (Abb. 10). Nach zweimaligem Aufheizen 
50- fa und Abkühlen ist die Suszeptibilität bei 
(d) Zimmertemperatur sehr stark ange- 


stiegen, und die Temperaturabhängig- 
keit ähnelt der eines Ferromagnetikums, 
wobei in der Gegend von 500 °K noch 
4 ein Maximum überlagert ist. Auffällig 
‘ad ist, daß sich sowohl der Curie punkt 
SQ (c) als auch die Suszeptibilität oberhalb 

; des Curie punktes nicht geändert haben. 
(b Daraus folgt, daß der Luftsauerstoff 

U keine ferromagnetischen Oxyde gebildet 
haben kann, sondern in das Schwefel- 

teilgitter eingebaut worden sein muß, 


Magnetisierung in willk. Einheiten 


10 da er sich so bemerkbar macht, als hatte 

a der Schwefelgehalt zugenommen. Die 
(a) Gitterkonstante dieser Verbindung hatte 
nach dem Aufheizen an Luft etwas ab- 


300 400 500 600 


Temperatur in °K genommen, wie es ebenfalls einer Er- 


Abb. 10. Anderung der Magnetisierung hos ng des Schwefelgehaltes entsprechen 
von FeS, 39 bei mehrmaligem Aufheizen 

an Luft. a und b: Erstes Aufheizen und Zur genaueren Überprüfung dieses 
Abkühlen, ce und d: Zweites Aufheizen Ergebnisses wurden noch zwei weitere 
By und Abkühlen Untersuchungen vorgenommen: 


A Rest des FeS, ‚, wurde in vier gleiche Teile geteilt. Zwei davon wurden 
im Vakuum, die anderen beiden Teile an Luft auf 400°C erhitzt. Die eine 
von den im Vakuum erhitzten Proben wurde dann sehr schnell, die andere 
langsam auf Zimmertemperatur abgekühlt. Mit den beiden an Luft erhitzten 
Proben wurde in gleicher Weise verfahren. Die Ergebnisse der vor und nach 
dieser Behandlung ausgeführten magnetischen Messungen sind in Tab. 4 
aufgeführt. 


Aus der Tabelle geht deutlich der Unterschied zwischen der Behandlung 
an Luft und im Vakuum hervor. Die Magnetisierung der beiden an Luft 
erhitzten Proben ist etwa die gleiche wie bei FeS, ,3, während sich die Suszep- 
tibilität der im Vakuum erhitzten Proben gegenüber dem Ausgangswert 
kaum geändert hat. Die kleinen Diskrepanzen zwischen den langsam und 
schnell abgekühlten Proben rührt möglicherweise von den in 4.2.4 IUMEROFNENEE 
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Tabelle 4 


Suszeptibilität des | 


Suszeptibilität bzw. Magneti- 
Ausgangsmaterials 


Art der Behandlung sierung nach der Behandlung 


A Re | 1 him Vakuum bei 400 °C, lang- | feldunabhängig 


sam abgekühlt Y; = 19,0 10% 
1 him Vakuum bei 400 °C, | leicht feldabhängig, nach Extra- 
| schnell abgekiihlt polationaufH=co 
1 han Luft bei 400 °C, langsam : 

1h an Luft bei 400 °C, schnell 17590,293 = 17,5 emEg™ KE 
abgekühlt 


Die zweite Untersuchung sollte darüber Aufschluß geben, wieviel Schwefel 
durch Sauerstoff ersetzt werden kann. Zu dem Zwecke wurden Proben der 
Zusammensetzung FeS, 99, FeS, 9; und FeS, 4) je 1 bzw. 2 Stunden lang auf 
400 °C an Luft erhitzt und dann der genaue Gehalt an Sauerstoff aus der Ge- 
wichtszunahme, die Gitterkonstanten, sowie die Änderung der u 
hilität bzw. der Magnetisierung bestimmt. Die Ergebnisse enthält die Tab. ; 


Ausgangssubstanzen 
were cla | x,bei293 °K 
FeS; oo 1,704) 12,9 - 104g! 
es, 1 1,667| 54 -10* 
bei 17,5 kOe 
1 hauf 400 °Can Luft 2h auf 400 °C an Luft 
Zusammen- | Zusammen- | 
setzung a 17800,298 setzung 
FeS, ‚0005.13 1,675 9,0 emEg™ Fe 00,17 | 1,655) 14,2 emEg 
1,668, 9,5 | 1,655 12,9 emEg™ 
FeS; 1090, 1,657 14°8 emEg- FeS, 10% 06 1,657 


13,1 emEg™ 


Von den röntgenographischen Ergebnissen sind der Übersichtlichkeit halber 
nur die Achsenverhältnisse angegeben, die jedoch zum Vergleich mit den in 
Abb. 1 enthaltenen Daten der reinen Eisensulfide ausreichen. Auf den De- 
bye-Scherrer-Aufnahmen sämtlicher Proben waren noch keine Linien 
eines Eisenoxydes zu sehen, was jedoch einen geringen Oxydgehalt nicht aus- 
schließt. 

Die c/a-Werte nehmen beim Erhitzen an Luft deutlich ab. Die erreichten 
Endwerte stimmen fast genau mit dem c/a-Wert von FeS,,; überein. Die 
Ergebnisse der magnetischen Messungen deuten darauf hin, daß die Magneti- 
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sierung mit wachsendem Sauerstoffgehalt etwas abnimmt. 
punkte liegen alle bei etwa 570 °K, d.h. etwas niedriger als bei den reinen 
Eisensulfiden. Nach dem zweistündigen Erhitzen scheinen die ersten beiden 
Proben schon etwas ferromagnetisches Oxyd zu enthalten, da die Suszeptibili- 
täten oberhalb des Curiepunktes höher liegen als vor und nach dem ein- 
stündigem Erhitzen. 


Alle aufgeführten Ergebnisse sprechen dafür, daß bei vorsichtiger Oxyda- 
tion ein Teil des Schwefels in den nichtstöchiometrischen Eisensulfiden durch 
Sauerstoff ersetzt werden kann, wobei die Zusammensetzung Fe,S, ,0, 
angestrebt wird. Erst bei stärkerer Oxydation wird Eisenoxyd gebildet. 
Es mag zunächst verwunderlich erscheinen, daß in den Eisensulfiden ein 
Teil des Schwefels durch Sauerstoff ersetzt werden kann, da der umgekehrte 
Fall — Ersatz des Sauerstoffs durch Schwefel in Magnetit — nicht möglich 
ist. Es dürfte aber leichter sein, das relativ große Schwefelion durch das 
kleinere Sauerstoffion zu ersetzen als umgekehrt. 


7. Ersatz des Eisens durch Kupfer, Nickel und Kobalt 


Vom FeS,,; ausgehend wurde versucht, das Eisen durch andere Uber- 
gangsmetalle wie Cu, Ni oder Co zu ersetzen. 


Beim Kupfer und Nickel gelang dieser Ersatz nicht. Schon bei geringen 
_ Gehalten an Cu oder Ni bestanden die Mischsulfide jedesmal aus zwei Phasen, 
wobei im Falle des Kupfers neben Pyrrhotin noch Chalkopyrit (CuFeS,) 


enthalten war, wahrend 
a-Achsel 3, im Falle des Ni als 
} : zweite Phase Pentlandit 
s 342 fein kubisches Eisen— 
5,60. 340  Nickel-Sulfid der Zu- 
| sammensetzung (Fe, 
I ° T 171338 Ni))S, mit 40 bis 75% 

540 - 166 Fe] auftrat. 
Ya CoS,;5 und FeS,5 
: ; 162 sind dagegen im festen 
f Zustand in jedem Ver- 
# hältnis mischbar. Zur 
14754 Untersuchung der Git- 


0 20 40 60 80 1000 20 40 60 80 100 terkonstanten und der 
magnetischen Eigen- 


Co- Gehalt in Atomprozent schaften wurden Proben 
Abb.11. Gitterkonstanten der Fe—Co-Mischsulfide der Zusammensetzung 


Fe,_,Co,8,,; in Abhängigkeit vom Co-Gehalt Fe,_,Co,,S) 15 herge- 
stellt, wobei n in Stufen von 0,1 zwischen 0,0 und 1,0 variiert 
wurde. 

Die ermittelten Gitterkonstanten sind in Abb. 11 als Funktion des Kobalt- 
gehaltes aufgetragen. Bis etwa 50 At.-%, Co nehmen sowohl a und ¢ als auch 
das Achsenverhältnis c/a ziemlich stark ab. Oberhalb von 50 At.-%, Co ändern 
sich c und c/a nicht mehr so stark, während a sogar praktisch konstant bleibt. 


Die für das CoS, ıs gefundenen Gitterkonstanten stimmen recht gut überein 
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mit Werten von Lundquist und Westgren!”) (a = 3,361 A und ce = 5,160) 
und von Vogel und Au!) (a = 3,385 A und c/a = 1,54). Die sehr große 
Fehlergrenze bei dem Sulfid mit 20 At.-%, Co rührt davon her, daß auf der 
Debye-Scherrer-Aufnahme nur sehr verwaschene Linien zu sehen waren, 
die keine genauere Auswertung zuließen. Auch an anderen, unabhängig 
voneinander hergestellten Proben der gleichen Zusammensetzung wurden 
keine besseren Aufnahmen erhalten. 


Die Abb. 12 bis 14 enthalten die Ergebnisse der magnetischen Messungen. 
In Abb. 12 sind die bei Zimmertemperatur und H = 17,5 kOe gemessenen 
Magnetisierungen eingetragen. Der gen verschwindet bei 
Zimmertemperatur schon bei etwa ” At.-% Co 
Die im Bereich von 50 bis 90 At.-°4 Co gemes- 
senen Suszeptibilitäten liegen sieh im para- 
magnetischen Bereich (22 bis 13 - 10% em3/g bei 
17,5 kOe); sie sind jedoch noch schwach feld- 
abhängig. Nur die Suszeptibilität des CoS, 1s 
war feldunabhängig; sie betrug 4,35 - 10% em3/g. 

Die Temperaturabhängigkeit der Magneti- 
sierung der ferrimagnetischen Proben zeigt die 
Abb. 13. Das angelegte Magnetfeld betrug wieder 
7,5kOe. Ob sich das bei den ersten beiden 
Proben noch vorhandene Maximum mit zu- 
nehmendem Co-Gehalt nach tieferen Tempe- 
raturen hin verschiebt oder ganz verschwindet, 
konnte nicht geklärt werden, da keine Messungen 
unterhalb von 100 °K ausgeführt werden konnten. 2 
Auffällig ist die bei 20 und 30 At.-%, Co auftretende ee Kohängickait 
thermische Hysterese am Curiepunkt (Abb. 13c vom Co-Gehalt 
und 13d). Sie war bei mehrmaligem Durchlaufen 
der Temperaturkurve völlig reproduzierbar. Die aus den o—T-Kurven extra- 
polierten Curietemperaturen enthält die Abb. 14. 
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Spez. Magnetisierung in emE pro 9. 
© 


0 20 40 60 80 100 


Co-Gehalt in Atomprozent 


Eine ausführliche Diskussion der magnetischen Eigenschaften des Systems 
FeS, ;;—CoS, ,, stößt auf erhebliche Schwierigkeiten, da sich die Zahl der 
beteiligten Ionenarten gegenüber den Eisensulfiden erhöht hat. Wir wollen 
uns daher mit einigen Bemerkungen begnügen, die den Übergang Ferrimagne- 
tismus—Antiferromagnetismus bei etwa 50 At.-%, Co betreffen. 


Der Antiferromagnetismus des CoS, „, spricht dafür, daß hier keine geord- 
nete Verteilung der Ionen und Leerstellen mehr vorliegt. Für die Leer- 
stellen wird dies durch die Angabe von Lundquist und Westgren!”) be- 
stätigt, daß sie im Homogenitätsbereich des CoS keine Überstrukturreflexe 
beobachten konnten. In den Mischsulfiden muß daher mit zunehmendem 
Co-Gehalt die antiferromagnetische Wechselwirkungsenergie im Verhältnis 
zur Ordnungsenergie zunehmen, und zwar so, daß von etwa 50 At.-% Co ab 
der antiferromagnetische Zustand stabiler ist. Betrachtet man die Gitter- 
konstanten, so fällt auf, daß c im Bereich von 0 bis 50 At.-%, Co sehr stark 


1) D. Lundquist u. A. Westgren, Z. anorg. allg. Chem. 239, 85 (1938). a 
8) R. Vogelu. R. Au, Z. Metallkde. 40, 290 (1949). 
“ 


§ 
1 
1 
D 
x 
| 
4 


294 20 & Annalen der Physik. 7. Folge. Band 5. 1960 


nd mer 


2 


N 


> 


E 
So 


€ 
; 8 8 
Q 
7 
\ 
N 
a 
1 0 
12 
6 


200 400 6000 20 400 600 
Temperatur in °K 


Abb. 13. Temperaturabhängigkeit der Magnetisierung der ferrimagnetischen Fe—Co- 
Sulfide bei H = 7,5 kOe 


600 abnimmt, d. h. die entgegengesetzt magneti- 
U sierten Schichten rücken näher zusammen. 
Damit wäre aber die Zunahme der magnetischen 
400 Wechselwirkung in diesem Konzentrationsgebiet 
\ zu erklären. Möglicherweise werden durch die 


Abnahme der Gitterkonstanten die Platz- 
wechselvorgänge zwischen den einzelnen Ionen 
\ erleichtert, so daß außerdem noch die Ord- 
nungsenergie erniedrigt wird. 


200 


Curietemperatur in °K 


0 : Besonderer Dank gebührt Herrn Prof. Dr. 

9 20 40 60 80 0 W. Holzmüller vom Physikalischen Institut 
Co-Genalt in Atomprozent der Karl-Marx-Universität Leipzig für wertvolle 
26 Anregungen und die Unterstützung, die er der 
‘e—Co-Mischsulfide in Abhän-. Arbeit zuteil werden ließ, sowie Herrn Dr. 
gigkeit vom Co-Gehalt E. Röttig für wertvolle Hinweise. 
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Gleichzeitig möchte ich auch Herrn W. Keilig danken, der mir bei der 
Durchführung der Messungen und deren Auswertung behilflich war. 


Jena, Institut für Magnetische Werkstoffe der Deutschen Akademie der 
Wissenschaften zu Berlin. 


Ei Bei der Redaktion eingegangen am 6. Juni 1959. 
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On the Statistical Theory ot 


By V. B. Magalinskij and Ja. P. Terletskij 


Summary 


A general method is worked out on the basis of principles of Gibbs’ 
statistical mechanics which allows to find stationary probabilities and transi- 
tion probabilities for physical quantities provided either the behaviour of 
their mean values or the general form of corresponding equations of 
motion (Langevin equations) is known. 


vil 


The proposed method is free from ordinary restrictions of the theory of 
fluctuations and Brownian motion, such as, modest fluctuations, linearity of 
phenomenological equations, Markov’s character of the random process and 
is applicable also to postacting systems. 


1) _ The method is applied to some physical systems. 


For the system, described by a non-linear equation of motion the general 


coordinate-velocity equation of motion of the transition probability density 
is obtained. This equation coincides with the well-known Einstein-Fokker- 
Plank coordinate-velocity equation in the case of linear friction law. The 
diffusion equation is obtained for the electrical circuit consisting of capacity 
and nonlinear resistance. Its general solution is found. 


It is well that, using methods of statistical we 
can derive exact relations which enable us to reduce a calculation of 
fluctuations and correlations of various quantities (time correlations also) 
to the determination of average values of these quantities in the presence of 
(or after turning on) additional constant forces!*), 


It was shown®) that all principal momenta determined in fluctuation 
theory, as well as in the theory of Brownian motion, can be computed by this 
method. 


1) J. W. Gibbs, Principles of Statistical Mechanics (Collected Works, 1), New York, 
Longmans Green and Co., 1928. 

2) V. V. Viadimirskij, J. Exptl. Theoret. Phys. (USSR) 12, 199 (1942). 

3) V. V. Vladimirskij and Ja. P. Terletskij, J. Exptl. Theoret. Phys. (USSR) 
15, 258 (1945). 

4) Ja. P. Terletskij, mHamuyeckne u 3aKOHbI (Dyna- 
mical and Statistical Laws of Physics, Moskow 1949). 

5) Ja. P. Terletskij, Nuovo Cimento 7, 308 (1958). 


M a 


This 
density 
tained it 

The : 
which n 
as well 
density 
systems. 

In tl 
the form 
(Lange 
the rand 


Let 1 
gate X ¢ 
to be so! 

Then 
F,(X) = 


mechani 


where W 
For 
exp 


6-functi 


one can 
Wi 


where a 


Z 
Zis 


On 
distribu 


| 
| 
V. 
729 (1958 
7) V. 
V. 
a 86, 1731 


Magalinskij and Terletskij: The Statistical Theory of Nonequilibrium Processes 297 


This method was worked out for the calculation of transition probalility 
density in®), and general diffusion equations for nonlinear systems were ob- 
tained in °)8). 

The systematical account of the general statistical method is given below 
which may be used to calculate directly the transition probability density 
as well as to obtain general equations of motion for transition probability 
density in coordinate-momentum (velocity) space for arbitrary nonlinear 
systems. 

In the latter case it is sufficient to make only one hypothesis concerning 
the form of an averaged phenomenological equation of motion of such a system 
(Langevin equation), and it is unnecessary to use any hypothesis about 
the random fore es correlation or about the character of stochastic processes. 


Let us consider a system, microstates of which are described by the ies 
gate X of all its canonical variables. Let us suppose F, (X), F, (X),..., F, (X) 
to be some functions of canonical variables of the system. 

Then the probability for the system to have given values F, (X) = F. 


= F,,...,; F, (X) = F,, according to general principles of statistical 
mechanics, is 
8 
W(F...... F,) = I 6 {F,, — F,, (X)} w (X) dX, (1) 


where W (X) is the probability density of the microstate X. 

For the system, according to Gibbs’ canonical distribution w (X) = 
exp [Y — H (X)]/6. In this case, making use of an integral representation of 
ö-funetions entering into (1): 


6{F, —F, (X)} = f exp {ié, [F, —F, (X)}} dé, ( 


bo 


) 
one can write (1) in the form 
W(F,...., F)=(2n)* I dé, exp Fy); 
where a, = iQ &,; | 
Z (a,,...,@,) (exp 5) [exp | + Fy (x) ax. (4) 
(x 


Z is a characteristical function for W, as 


On the other hand, one can consider the integrand in (4) as a Gibbs’ 
distribution for a system with the Hamiltonian 


(X) =H (X)+ (X). (6) 


hy V.B. Magalinskij and Ja. P. Te a J. Exptl. Theoret. Phys. (USSR) 34 
129 (1958); Soviet Phys. J.E.T.P., 7, 501 (1958). 

?) V. B. Magalinskij, J. foxptl. Theoret. Phys. (USSR), 36, 1423 (1959). 

*) V.B. Magalinskij a. Ja. P. Terletskij, J. Exptl. Theoret. Phys. (USSR) 
86, 1731 (1959). 
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Differentiating (4) with respect to a, we obtain Lint ing 


ny Ns 6) Z 2} 
(5') 
Z day? daze 


The index “a” at the averaging line means that the averaging is carried out 
over an equilibrium distribution with the Hamiltonian function (6) (over 
the state ‘‘a’’). 

Particularly, if F,—Q, are generalized coordinates, the state “a” 
is such a state, if an additional constant force — a, acts in the direction 
of each generalized coordinate Q,. 

Let us assume that we know (either from experience or from theoretical 
considerations) the macroscopic equations of state, to which the quantities F, 
obey. Differently speaking some connections between various moments of 
these quantities are known, which are determined by relations of the form 


O(F,,.. 57, Fi,...,.)=0. (7) 


Then, writing these equations for the state “a” and making use of relations 
(5’), we obtain differential equations for the statistical integral Z: 


(-7& = 9 82 (—6)? #2 )=0 
Zu Din Zur ba, =0. (7) 


One can obtain Z from previous expressions and then, according to (3), 

As a simple example, we compute the distribution of coordinate q of a 
particle in an external field of force K m = — dU/dq = = im 


equilibrium condition in the state “a” is: 
(3) 


Aye: 


K (q) = =K = 


or, according to (5’), 


(—6 


Going over from Z to W by means of Fourier-transformation, we find 
K (q) W = —@" dW/dq, where W (q)= W,exp [—@" U (q)], i.e. the 
Boltzmann distribution. 

The problem becomes easier, when the equation of state connects only 
first-order moments, i.e. mean values of quantities F,. In this case equa- 
tions (7’) are integrated directly, and we get 


= exp AY (a) =f 


0,...,0 
+00 = 
W (F,. F,) (2 f eve [exp iz &. F, AY (a, a,)/0| IT 


(10) 


Equation (10) enables us, particularly, to determine the probability distri- 
bution of generalized coordinates from empirical data which give the depen- 
dence of its mean values on the external force in an isothermal N: Exam- 
ples of the application of this equation see®). Sa Br 
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Let W (F, t; Fo, to) be the probability density of the transition of the con- cae 
sidered system m the state ; 
F (X) = F, (i. e. F, (X) = Fy; Er F, (X) = F,o) = > 

at the initial instant t, to the state F (X) = F (F,(X)=F,...,F,(X)=F,. 
at the time ¢, and let W,(F) be the probability density of equilibrium 
values F over the Gibbs’ ensemble [see (1), (10)]. a 
Then, according to*), we have 


} 6 {F, — F (XP) dX°®, 
(12) 
where X and X° are canonical variables of the system, connected by means 
of relations x(t; X°,t,), which represents solutions of dynamical 
equations with the Hamiltonian function H (X). According to (12), when _ 
[=t), obviously, W (F, ty; Fo, to) = 6 (F — Fy). 
If at the initial instant ¢ = ¢, F isn’t given exactly, but with some proba- 
bility density o(F),then the probability density at an arbitrary time t is defined 


W(F,t; Ft) = = [dir -F (x9) 


as 
(13) 
0 AN ay 
F — F (X49} exp — H (X°) aX". 
The equation (12) will be obtained from (13) if o(F) = 6(F — F,). 
Similarly to*), one can also represent the expression (13) i in the form a 
where @ (£, t) is the characteristical function: PETE a 
(é.t) = R {F (X exp {— i + [W — H (XP axXe, (1) 
R (F) = 0 (F)/W, (F). (16) 
Moments of the distribution (14) are expressed in terms of p as 
Let us introduce the function be su 
Z (a, t) = R {F (X°)} exp {— [H (X°) + aF (X96}dX°. (18) 
(19) 
On the right hand of expression (18) one can also substitute H (X*) inste ad 


of H (X°) according to the energy conservation and dX instead of dX° accor- 
ding to Liouville’s theorem, and then, taking into account Xt = x (t; X", a MM 
and X° — x (tj), X”, rt), one can consider the expression 


w (X41) =Z REF (X)} (A(X) +aF 
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as a microscopical probability density over some non-equilibrium process 
which obeys the equation of motion of the cone ensemble 
) 


Obviewty. moments over the ensemble (2c 0) are defined by means of 
relations, similar to (5’): 
Z dal! ... dare 
“= 0, 2) coincides with (17). 
Moreover, Pia EN Z with respect to ¢ we obtain the relation 
Provided moments (22), (23) computed by means of the non-equilibrium 
ensemble (20), are known, one can to compute Z and, consequently, W. 
All the information about this moments, which may be obtained only from 
the macroscopical experience, one can divide in two categories: 
a) we know the behavior of mean values of quantities F,,..., F, or the 
—a 
behavior of corresponding mean velocities F,, ... , F, in the non-equilibrium 


process corresponding to (20). 

b) We know the relation between various moments of quantities F,,..-F, 
which arises as a result of averaging of well known phenomenological equations 
for this quantities (Langevin equations) over the ensemble (20). 

(It is well known, that for linear systems cases (a) and (b) are equal.) 

According to this fact the problem arises how to calculate the transition 


probability density either on the basis of given mean values (F* er F, 


—a 
or F,,..., F,), or given Langevin equations, to which the quantities F obeys 

In the first case the transition probability density is computed by direct 
integration of equations (22), (23), in the second one the solution of the pro- 
blem is reduced to the construction of an equation of motion for the transition 
probability density. 

In the following sections we give the solution of the first problem as well 
ae as of the second for particular physical systems. 

r ie IV. Calculation of the Transition Probability Density for a Generalised 

2 


Coordinate according given Mean Velocity of its Change 


Let F, =@Q be some generalised coordinate (from now on we confine 
ourselves to the one-dimensional case for simplicity). 
Then the expression (20) for the microscopical probability density at the 
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One can consider this expression as a non-equilibrium probability 
density for the phase ensemble, which is formed from the equilibrium one 
with the Hamiltonian function H (X!) + aQ (X*) by means of definite fixation 
of canonical variables X‘ at the time ¢ being given by the function R {x (ty; X, t)}. 

Obviously, one can consider the equilibrium system with the Hamil- 
tonian 


H+aQ (25) 


as a system being under the action of an additional constant force — a. If 
o= Wo, i.e. R= 1, the ensemble (24) with the Hamiltonian (25) is an 
equilibrium one, and the mean value of velocity Q is zero. But if o + Wo, 
—a 
the equilibrium is disturbed and the mean value Q, obviously, is different 
from zero. Let us suppose ®) one can obtain this mean value from some ma- 
eroscopical experiment with an additional real force — a. 
Let us suppose we know the mean velocity of change of the coordinate in 
the transition, under the action of force —a. 
4 9° 
The relation (23) gives for Z 
Za =aQ. (26) 


AY=—-HAlnZ 


According to the definition of Z 


(26) 


+00 
W (Q, t; Qo; to) = (? exp {i EQ — AY (a, t)/6} =i 66. (27) 
From the normalizing condition it follows 


AP (0, t) = 0. (28) 


If W obeys the initial condition W (Q, ty; Qo. to) = 6 (Q — Qo). we have also 
AP (a, th) = aQy.- (29) 


Integrating (26) and taking into account (28), (29), we obtain 


So formulas (27) and (30) enables us to determine the transition probability 
density for an arbitrary generalized coordinate, provided the mean velocity 
of change of latter in the presence of an additional constant force is known 

phenomenological equations or from empirical data. 


— 

®) It is possible that one can prove this hypotesis or some equivalent one if one can 
e more distinctly the ,,macrocopical quantities“ and variables entering the Langevin 
équation. In present paper we put this hypotesis forward as some rather obvious one and 
justifying oneself by obtaining results. 
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ern " The ‘Equation of Motion for Transition Probability Density in 
Configuration Space 


Let us consider the non-inertial brownian particle being under the action 
of an external potential force K (q) and a friction force — y (q) q, where y (q) is 
the coefficient of friction, depending in general from the coordinate q of the 
particle. Therefore, q obeys the equation: 


vy = K (q). (31) 


a Let us suppose a motion av eraged over the ensemble (24) to be described 
by the averaged equation (31), i. e. 

The following procedure is typical for the applied method. We evaluate 

the right and left sides of the equation (31) into power series. After this we 

average each term of this evaluation, pickung out in the explicite form the 

connection between different moments of coordinate q: 


(31! 


(33) 


(33) 


which we get by successive differentiation of Z/a with respect to ¢ and a, we 
obtain an equation for Z: 


Finally, using the connection between Z and W, and after applying a 
Fourier-transformation to this equation, we obtain an equation for the 
transition probability density: 


or 
ow 


So, using only general laws of Statistical mechanics we obtain the well 
known Einstein-Fokker-Planck equation for a system, described by 
the equation (31). Here it is interesting to note that the introduction of @ 
constant force —a formally introduces a diffusion flow on the right hand of 
the equation (36). 
; The consideration given above was restricted by a non-inertiality of the 
ye system. In the following part we will get rid of this restriction, carrying out 
the consideration in a phase space. 
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VI. The Equation of Motion for Transition Probability Density in the ps a 


Phase Space 


4 


Let us introduce in addition to the generalized coordinate FR, = @Q the 
corresponding generalized momentum F, = P and let us consider the transi- 
tion probability density in the phase space. In this case the microscopical pro- 


bability density (20) at the time t = t has the form 
w (X¢, t) R {Q (X°), P(X} exp|— 5 TH (X) + aQ(X*) + b P( 


(a, =a, a, = 5). 


)} 
(37) 


If we consider the momentum P to enter into the Hamiltonian H only as 


aterm P?/2 M, then 


(P+Mb% MB (38) 
= H,+ + aQ, 
and, consequently, for Hamiltonian (38) the canonical variables P and @ 
satisfy the equations 
P=- a. (39) 


Further the macroscopical equations of motion of particle (equations averaged 
over molecular parameters, or Langevin equations) may have the form 


P=F(Q,Q), Q= P/M, 


(40) 


where F is a sum of the external potential force and the dissipative force 


depending arbitrary from Q as well as from Q. 
We suppose, analogously to § 4, that the motion, averaged over 
semble (37), obeys averaged equations (40), i. e. according to (39): 


—a,b — a,b —a,b 


the en- 


u, —a,b 
P =F(@Q) =P |M+b. (41) 
On the ER relations (22) and (23) give in our case for a 
da” ob" 
43 
Z ot (43) 


Now, combining (41), (42) and (43), one can construct the in of 

motion for the transition probability density in an phase space similarly as 

it was carried out in the previous part of this paper. a EEE A to 
Substituting (41) into (43), we obtain: EST 0 


a P 0 82 
z + = 
Evaluating F into a power series, we have LE, 
P P\m 
+ SF b) = - 5 or (zz) bi-m 4! 
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where om ae I!/m! ( 1—m)!. Substituting (45) and using (42), we otain the 
equation for Z: 


Since, according to (14) and (19), 

W (Q, P; t) = (27) exp (¥/0) [[ Z t) exp [(aQ +6 P)/O]dadb, (42) 


then, applying to (46) Fourier-transformation (47), we obtain an equation 
for W: 


ow ¥ oW _ x FC cr “Qi "(pm W). (48) 
~ OPEL m kl M” apı-m 

With (J—m)! Cj" Pm = d!-m Pl/dP!-”, according to (45), one can re. 

present the right-hand side of equation (48) in the form: 


oP (!— m)! ap!” op!" 


(49) 


If we A Kai a ity V = P/M instead of momenta P, and take into 
account (49), we can write the ae (48) in the form: 


So we have obtained the most nn equation of diffusion in the coor- 
dinate-velocity space for a system with an arbitrary macroscopical dissipative 
force and a potential force. 

Obviously one can obtain an expression for the transition probability 
density W (Q, P, t; Q5. Po; fo) as a solution of the equation (50) if the initial 
condition is 


W (Q, P, ty; Qo, Po to) = 6 (Q — Qo) 6 (P — Po). 


(51) 


the equation (50) becomes the usual coordinate-velocity Einstein-Fokker- 
Planck equation, which has the form 


ow py ow y @ ew 
(see, for example ®)). 

Unlike all well-known methods of derivation of equations for the proba- 
bility density, equations (36), (50) and (52) are obtained by using only the 
general principles of statistical mechanics and the natural hypothesis (31’), (41) 
about averaged macroscopical equations of motion. = =——> 
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Herein we did not make any assumptions about the Marcovian character 
of random processes or about the character of correlation of random forces. 
Consequently, equations (36), (50), (52) are obtained from most general 
assumptions. 


VII. The exact Solution of Equation (50) for one Special Case 


When an external force and a friction force does not depend on coordinates, 
i.e. F(Q,V) = K —G(V), one can integrate equation (50) exactly. At 
first we write this equation in new variables: 


E ad ew „er 4K Kew ey og 


where g = Q/0,u = V = (M6), B(u) = B, w=G (u Y6/M). 


The expression for Y one can transform into the tenn 


j-2k 
Y = Y, exp (—u?/2), Y, = - (54) 


and for W, = W (exp (u?/2 — Kq)) we have ae ask i a 
+u ‘4 K = uY;. (55) 
OT öq eu eu Daw 
After substituting ML 


—2 oq+su a6 
= (221) W, (a,8;t)dods, (56) 
we obtain for W, 
OW, 


Ks W,= 8 WwW, 4 W, Ar (57) 
where 
Y,(s) = ; H, (sfi V 
(H,, (x) is the n-th Hermite polynomial). 
As a result of the substitution 


we obtain an equation 


the solution of which, satisfying the initial condition 
W, (a, 8; 0) = Way (0, 8) 

Ann. . Folge, Bd. 
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(0, 852) = Way fo, 
where p is defined by 
J Y,(® +0 4 ray (63) 
When W (q, u; 0) = 6 (¢ — 9) 6 (u — up), Wa is 
Wo = exp {ug/2 —8U, + + K) qQ} (64) 


by means of (59). 

Returning from W, to W and from the variables g, u,t to variables Q, V ,t, 
we obtain the expression for Green’s function for equation (53), which may 
be written as following: 


W(Q— Qo, ts = 
= | — (V — V,)* +2 (Q- — Qe) do 


Here Q, (a, &;t) and V, (a, &;¢t) are coordinate and velocity, computed 


from the equation 
MV,+4G,(V,)=K-—o (66) 


with the initial conditions Q,|;~9 = 0, V ole= o = € and 


M(V-—V 
G,(V) - Vix G (V ) exp |— (67) 


It is easy to see, that if 6 > 0 
W (Q—Qp, V.t; Vo) > 6 [Q — Q Vo; 0 [LV — V (Vo, (68) 


where Q (Qo, Vo; t) and V (Vo, t) are solutions of the non-averaged equations 
of motion (Langevin equations) which determines the motion of a Brownian 
particle, i.e. equations 
P-+6(Q)= K, Q= P/M (40) 
with the initial conditions 
hr, Qlı-0o = V leno = Ve 


_ Hence it is easy to obtain the asymptotic expression for W if 9 is small, 
evaluating Q, and V, in power-series of § — V, and o. 


VIII. The Application of the Method to an Electric Circuit, Consisting of 
Capacity and Nonlinear Resistance 


Following the given method, it is easy to obtain equations, similar to the 
equation (50) also for non-mechanical physical systems choosing properly in 
each special case generalized coordinates, velocities and forces. 

Let us consider as a simple example the electric circuit, consisting of a 
capacity and a resistance with a non-linear current-voltage characteristic 


J =F (V)= SF, Vi. 
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The charge Q concentrated in the capacity, is a generalised coordinate in 
our case, the voltage V is a corresponding generalized force, and consequently 
the averaged equation of motion for@ in presence of an additional constant 
force — a is 


a 

Q+F(Q/C+a) =0. (69) 
By means of the procedure described above we find the differential equation 
for Z (a, t): 


22 ct (— 0/0): (70) 
et 6 af ca 
and then for W: 
er ög kl eq dq 


» 


/ yt 7 ark 7 
= F , 


This equation is similar to the equation (53), and it can be integrated exactly. 
The Green’s function is given by the expression 
+#00 
1 
Qo) = CO? dE, (72) 
-i00 
where Q, satisfies the equation 
with an initial condition Q,|,-0 = &:; 


Similarly one can treat other, more complicated systems. 

Appendix. It follows from a more thorough analysis of equation (50). 
that if the friction force F(Q, V) is an odd function of velocity, the obtained 
equation (50) satisfies general requirements of statistical theory for none- 
quilibrium processes of the Brownian movement type, namely it has as a 
steady-state solution the Maxwell-Boltzman distribution which is also 
the limiting distribution for t + oo. ‘a 


#0) S. Chandrasekhar, Stochastic Problems in Physics and Astronomy. 


Bei der Redaktion eingegangen am 19. Juni 1959. 
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mee. or Beiträge zur Klärung des Mechanismus 
Von Eugen Badareu, lovitzu Popescu und Iancu lIova 


In der vorliegenden Arbeit werden eine Reihe von theoretischen und 
experimentellen Ergebnissen gegeben bezüglich der Vorgänge, welche bei 
Glimmentladungen in Argon an doppelter und einfacher Kathode stattfinden. 
Es wird die Wichtigkeit der Rolle unterstrichen, welche der Resonanz- 
umladungsprozeß in der Kinetik der positiven Ionen in eigentlichen Gasen 
und, im besonderen, in der Theorie des Kathodenfalls spielt. Es wird auch 
gezeigt, wenigstens in den Bedingungen anomaler Entladungen, daß ein 
Großteil der Anregungen und Ionisierungen vor der Kathode in negativem 
Licht stattfinden, wobei diese Formation ein praktisch äquipotentiales Plasma 
darstellt, das von Primärelektronen durchzogen und unterhalten wird. 


F 


Ja der Doppelkathcdeneffekt bekannt ist, kann er in Kiirze in folgender 
Weise beschrieben werden: Legt man zwischen eine Anode und eine aus zwei 
gegenüberliegenden parallelen Scheiben gebildete Kathode eine Spannung an, 
die in dem betreffenden Gas bei genügend vermindertem Druck eine Glimm- 
entladung hervorruft, dann wächst unterhalb eines gewissen Wertes des 
Produktes p-a (p ist der reduzierte Druck und a der Abstand zwischen den 
Scheiben) die Stromdichte rasch an und erreicht für genügend kleine p-a 
Werte, die um 2 bis 3 Größenordnungen höher liegen als jene bei der gewöhn- 
lichen Entladung (mit einfacher Kathede) unter den gleichen Bedingungen. 
Gleichzeitig wächst auch die Lichtstärke jenes charakteristischen Lichtge- 
bildes stark an, die aus der Verschmelzung der beiden negativen Glimm- 
lichter entsteht. 

Der keschrietene Effekt spielt sich in ähnlicher Weise auch an einer 
zylirdrischen Kathcde ab, und im allgemeinen bei jeder Hohlkathode (für 
vellstindige Literatur siehe 2)). 

Was das Verständnis der Vorgänge betrifft, die sich an der Kathode einer 
Glimmentladung abspielen, ist man stufenweise zu einem vollständigen und 
qualitativ befriedigenden Bild gelangt; von den unzähligen Beiträgen wollen 


1) Vv orgetragen auf der Tagung der Deutschen Physikalischen Gesellschaft in Leipzig 
am 17. April 1959. 


2) P. F. Little und A. v. Engel, Proc. Roy. Soc. London (A) 224, 209 (1954). 
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wir nur 2) hervorheben. 7) enthält eine Aufzählung der bis jetzt vorge- 
schlagenen Theorien des Kathodenfalles. 

In der vorliegenden Arbeit versuchen wir einen Schritt vorwärts in Richtung 
einer quantitativen Erklärung des Mechanismus der Vorgänge bei der Doppel- 
kathode zu machen; wie wir sehen werden betrifft dies auch die Vorgänge an 
der einfachen Kathode einer anomalen Glimmentladung, und zwar werden 
wir eine Reihe von Ergebnissen besprechen, die wir in Argon erhalten haben. 
Zu diesem Zweck wollen wir im folgenden Abschnitt einige Beziehungen an- 
führen, die die Grundlage unserer Abhandlung bilden werden. 


2. Theoretischer Teil 


Auf Grund jüngster Forschungen*~!) ist man zu einer Beziehung gelangt, 


die die Fortschreitungs-Geschwindigkeit v der positiven Ionen im eigenen Gas, 
als Funktion von E/p (E ist die elektrische Feldstärke und p der reduzierte 
Druck) und vom Umladungs- Querschnitt Q, richtig ausdrückt 


=) (2 (1) 


(m und e sind die Masse bzw. die Ladung des Ions). Diese Beziehung ist fiir 
E/p größer als 100 Volt/cm mm Hg gültig; sie ist auf Grund des Modells der 
Bewegung des Ions als ,,anisotrope Stafette‘‘ abgeleitet. In dieser Hinsicht 
werden wir die Tatsache unterstreichen, daß die Kinetik der positiven Ionen 
im eigenen Gase praktisch nur vom Umladungsprozeß bestimmt ist. Diese 
Hypothese, die vor verhältnismäßig kurzer Zeit aufgestellt worden ist!!) und 
sodann bestätigt wurde §) 9) 12-15), ist eine Folge der für den Stoß zwischen Ion und 
Atom gleicher Art charakteristischen Winkel- und Energieverteilung!®). 

Die quantenmechanische Berechnung des Umladungs- Querschnittes führt 
zu einem Ausdruck folgender Form: 


) Q, = A — Blog v. (2) 
(4 und B sind Konstanten, die vom ersten Ionisierungspotential abhängen); 
dieser Ausdruck bewährt sich in einem breiten Geschwindigkeitsgebiet 
(10°—108 cm/s) für eine ganze Serie von Dämpfen und Gasen, unter denen 
sich auch das Argon befindet !7). 


3) A. Güntherschulze, Z. techn. Physik 11, 49 (1930). rei a 
4) A. Lompe, R. Seeliger u. E. Wolter, Ann. Physik 36, 9 (1939). h Eee 


5) W. Weizel, R. Rompe u. M. Schön, Z. Physik 112, 339 (1939). 

6) M. J. Druyvesteyn u. F. M. Penning, Rev. mod. Physics 12, 88 (1940). 

?) G. Francis, Handbuch der Physik, B41. XXII, S. 53 (1956). 

8) Yu. M. Kagan u. V. I. Perel, Dokl. Akad. Nauk. 98, 575 (1954). 

®) Yu. M. Kagan u. V. I. Perel, J. Exp. Theor. Phys. 29, 884 (1955). 

0) T. Popescu, Rev. de Phys., Acad. R.P.R. 4, 211 (1959). 

1) L. A. Sena, J. Exp. Theor. Phys. 16, 734 (1946). - : aire 

2) B.M. Paliuh und L. A. Sena, J. Exp. Theor. Phys. 20, 481 (1950). cae 

13) E. Badareu u. M. Hagiescu, Bul St., Mat., Fiz., Chim. 603 (1950). ¢ 

4) E. Badareu u. M. Hagiescu, Stud. Cere. de Fiz., Acad. R.P.R. 3, 429 (1959). a 

15) I. Popescu, Rev. de Phys., Acad. R.P.R. 4, 199 (1959). 

16) H.S. W. Massey u. E.H. S. Burhop, Electronic and Ionic Impact Phenomena, 
Clarendon Press, 1956. 

Fi Popescu u. N. J. Ionescu, J. Techn. Phys. 29, 866 (1959); Proc. Phys. Soc., 
zur Veröffentlichung abgesandt (1959). 
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Das Auftreten der Raumladung und des elektrischen Feldes vor der 


Kathode ist eine Folge der durch (1) beschriebenen Fortschreitungsbewegung 
i der positiven Ionen sowie der Poissonschen Gleichung, die man fiir die ebene 
Er: Kathode folgendermaßen schreiben kann: 
u” (3) 


dabei geht die x-Achse vom Rand des negativen Lichtes in Richtung der 
Kathode aus; 7 ist die elektrische Stromdichte. Diese Art die Poissonsche 
4 Gleichung auszuschreiben ist durch folgende Tatsachen gerechtfertigt: 
EN a) die Konzentration der Elektronen ist infolge ihrer großen Geschwindig- 
keit im Dunkelraum unbedeutend gegenüber der Konzentration der positiven 
Ionen; 
Se b) der elektrische Strom wird im Dunkelraum praktisch nur von den im 
negativen Glimmlicht gebildeten positiven Ionen getragen; dies ist eine Folge 
der Tatsache, daß die von den Elektronen hervorgerufene Anzahl von Ionisie- 
oy rungen im Dunkelraum (wie wir weiter unten zeigen werden) sehr klein ist 
und die Elektronenemission aus der Kathode sehr schwach ist. 

ü Wenn wir den Ausdruck (1) von v in (2) einführen, dann erhalten wir nach 
2 der klassischen Methode?°) eine einfache Gleichung, die durch Integration mit 
der experimentell gut begründeten Randbedingung 2—=0, E=0 uns zu 
folgendem Ausdruck des elektrischen Feldes im Dunkelraum führt: 

6 


Die Vereinfachungen, die man um den Ausdruck (4) zu erhalten machen 
muß, sind folgende: 
me a) Q, ist im betrachteten Geschwindigkeitsgebiet unabhängig von v¥), 
in. eine Vereinfachung, die durch die Form der Beziehung (2) gerechtfertigt ist. 
Deshalb werden wir in allen folgenden Beziehungen unter Q, einen Durch- 
Br. schnittswert für den Dunkelraum verstehen. 
Br: b) Der Ausdruck (1) für v ist nur dann richtig, wenn die Ionen sich im 
Gleichgewicht mit dem Feld?) bewegen, was in sehr großer Nähe des Randes 
des negativen Glimmlichtes nicht mehr gilt; eigentlich strebt in diesem Ab- 
sehnitt zu schwacher Feldstärke die Bewegung der Ionen zur Isotropierung 
und » wird allmählich proportional E!p 
> Die Integ (d = Dicke des Dunkel- 
s raumes), ergibt den Kathodenfall aT" auf diese Weise erhält man den gewöhn- 
Be lichen Ausdruck zwischen den Ähnlichkeitsparametern j/p?, pd und V.: 


= { —-- (5) 
3 P | 18 ( m Q: (p pay" 


hut 
grundlegenden Abänderungen an dem folgenden bedingt. 
1%) Das Anwendungskriterium von (1) ist in erster Näherung gegeben durch 


<8E?®), 

a 2, Yu. M. Kagan u. V. I. Perel, Dokl. Akad. Nauk. 108, 222 (1956). aii 


18) Es kann gquecigt werden, daß die Inbetrachtziehung der Abhängigkeit von v keine 
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Der Ausdruck (5) ist dem der in 2) abgeleiteten ungefähr ähnlich, der 
Zahlenfaktor auf der rechten Seite ist jedoch hier ~ 3,5mal kleiner. Unter 
Zuhilfenahme von (5) kann man den Ausdruck der Feldstärke (4) auch folgen- 
dermaßen schreiben: 5 Ve 


3d! (*) 

zum Vergleich erinnern wir an den Näherungsausdruck, der manchmal in der 
Theorie des Kathodenfalls angewendet wird?) 


Wenn man (4) von 0 bis x integriert, erhält man die a ai Potentials 
V(x) im 


amQ, 


eine Beziehung, die mit Hilfe von (5) auch in der Form td 


v=v.[1-(5) ] (6) 

geschrieben werden kann. 
Beziehung (4) kann noch zusammen mit der Poissonschen Gleichung (3) 


zu einem Ausdruck für die Konzentration der positiven Ionen n = j/ve im 
Dunkelraum führen: 


amQ, p / 

Natürlich hört der Ausdruck (7), in sehr großer Nähe des Randes des 
negativen Glimmlichtes, auf, gültig zu sein, da in diesem Gebiet die Bedingung 
«= 0, E=0, die zur Integrierung von (3) angesetzt wurde, gemäß der 
Beziehung (1) zu v = 0, also n = oo führt; eigentlich ist die minimale Ge- 
schwindigkeit der Ionen die Geschwindigkeit vmin # 0, so daß die maximale 
Schätzung von n(0) durch j/e vmin gegeben ist. 

Geht man zur Bon hnung von n von den Beziehungen (8 )) und (7) aus, dann 


erhält man 
18 ned": 


Dieser PRATER ist grundsätzlich verschieden von jenem, idee aus ik linearen 
Näherung des Feldes (4’’) folgt und der zu einer konstanten Konzentration 


der Ionen im Dunkelraum führt: 


J 

Wir bemerken noch, daß der aus (7’) folgende Mittelwert von x im Dunkel- 
raum, in der Nähe des Wertes (7’’) liegt: 

— 5 Vo ZZ 

Was die Geschwindigkeit der positiven Ionen im Dunkelraum betrifft, 

erhilt man diese sofort: v = j/ne. Daraus folgt unter Benutzung von (7): 
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(8) und (5) führen sodann zu KR Aue Pr 
3mm Q, 


R Mit Hilfe von (8) kann man die Energie der Ionen berechnen: 


woraus zu ersehen ist, daß die Energie der Ionen an der Kathode 


5 e Vo 
W u 3x Q, pd (8 ) 


ist; sie ist somit ungefähr der Hälfte der auf dem letzten, freien, mittleren 
Umladungsweg gewonnenen Energie gleich. 

Die Anzahl N, der von einem den Dunkelraum durchlaufenden Ion 
erlittenen Umladungen kann man einfach 


N,=Q,pd (9) 
schreiben; Q, ist, wie oben gesagt der Mittelwert des Umladungsquerschnittes 
im Dunkelraum. 

Die Anzahl N, der Ionisierungen, die von einem von der Kathode aus- 
gehenden Elektron ausgelöst werden, kann man mit Hilfe der Beziehung: 
d 


N, = pf Q; (x) dx (10) 


berechnen, wobei man Q; (x) durch Eliminieren von V, zwischen der durch das 
Experiment gegebenen Funktion Q,(V) und dem durch (6) oder (6’) gegebenen 
V(x) erhält. 

Auch ist die von einem primären Elektron zusammen mit den sekundären, 
tertiären usw. Elektronen im Dunkelraum ausgelöste Zahl N, der Ionisierungen 
sehr leicht angenähert durch den Ausdruck: 


gegeben, wobei Q, der Mittelwert des Ionisierungsquerschnittes im Dunkel- 
raum ist: 


Die Energie 7, die von einem primären Elektron und seinen im Dunkel- 
raum gebildeten Nachfolgern an das negative Glimmlicht abgegeben wird, 
kann man unter Vernachlässigung der Verluste durch elastische Stöße und 
=e Anregungen im Dunkelraum unter die Form 


(12) 


V (2792 — 


Er 


bringer 


ist. 

Im 
zweite 
Beitrag 
durch |] 

Unt 
negativ 


ausdrüc 

Die 
ein 
gelöst ' 


wo € = 
fiir ein 
völlig \ 
N, und 
hältnis 
Strom 


= d 

Der 
seinen | 
tiven ( 


Bez 
maßen 


wo die 
der Ph 
Atome, 
der Ke 
beuten 
Nachfo 
usw. ; f, 
auffällt 
Im 

die Ion 
rufen, 


= 
f 
E 
| 
@ 
4 


(12) 


Im Ausdruck (12) bedeutet e V, den Beitrag der Primärelektronen, das 
zweite Glied den Beitrag der Sekundärelektronen, das dritte Glied den 
Beitrag der anderen Elektronen aus der Lavine und das vierte die Verluste 

durch Ionisierung. 
Unter Benutzung von (12) kann man die mittlere Energie V 


ist. 


m der in das 


negative Glimmlicht eintretenden Elektronen durch % 
T 


ausdriicken. 

Die Anzahl N, der Ionisierungen, die in dem negativen Glimmlicht durch 
ein primäres Elektron und seine im Dunkelraum gebildeten Nachfolger aus- 
gelöst werden, ist 


N, =(1+N,)— —, (14) 


Vv. 
woe = e(V „) die von einem Elektron mit einer Anfangsenergie e V „im Mittel 
für eine lonisierung verbrauchte Energie darstellt (diese Energie wird 
völlig verbraucht). Unter Benutzung der oben berechneten Ausdrücke von 
N, und N, kann man am kathodischen Rand des negativen Lichtes das Ver- 
hältnis zwischen dem von den Elektronen getragenen Strom und dem totalen 
Strom finden 
(f; = der Bruchteil von N,, der auf die Kathode auffällt). ar 
Der Bruchteil (in Prozenten) der Energie 7’, die von einem Elektron und 
seinen im Dunkelraum gebildeten Nachfolgern für die Jonisierungen im nega- 
tiven Glimmlicht verbraucht wird, ist 


Ne V, Se 
T ion = 100 (16) 


Bezüglich der Stationaritätsgleichung kann diese im allgemeinen folgender- 
maßen geschrieben werden‘): 


(N + hi N;) Yi + thm N m Ym NaYn = |, (17) 
wo die vier Glieder der linken Seite der Reihe nach den Beitrag der Ionen, 
der Photonen, der metastabilen Atome und der beschleunigten neutralen 
Atome, die aus dem Umladungsproze8 stammen, zur Elektronenemission aus 
der Kathode bedeuten; y;, ;, 7. Y„ bezeichnen die entsprechenden Aus- 
beuten und N,, N,, N, die in Glimmlicht von einem Primärelektron und seinen 
Nachfolgern aus dem Dunkelraum gebildete Anzahl von Ionen, Photonen 
usw.; f;,f, und f,, bedeuten den Bruchteil von N,, N,, N „, der auf die Kathode 
auffällt. 

Im Falle des Doppelkathodeneffektes und unter der Annahme, daß nur 
die Ionen und die Photonen die Elektronenemission an der Kathode hervor- 
rufen, vayelenont sich die Beziehung (17) zu: 


(Ny +N)y, + Nyy 
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3. Die Versuchsergebnisse und ihre Diskussion 


Im folgenden werden die Ergebnisse der Versuche, die an einer Glimm- 
entladung mit Doppelkathode in Argon erhalten wurden, diskutiert werden. 

Tab. 1 enthält in den ersten vier Reihen die Werte der vier Größen 
V.. j, p und d; die Ablesungen geschahen gleichzeitig, d. h. im gleichen Zu- 
stand der Oberfläche der Kathode und in kurzen Zeitabständen (Sekunden- 
bruchteile), um die thermischen Effekte auszuschalten. j/j, bedeutet das Ver- 
hältnis zwischen der Stromdichte an der doppelten Kathode und jener an der 
einfachen Kathode unter gleichen Bedingungen, d.h. gleichem Druck und 
gleichem Kathodenfall. Die Versuche wurden in einem kugelförmigen Rohr 
aus Pyrex mit einem Durchmesser von 10 cm durchgeführt, die im Zentrum 
gelegenen und auf der hinteren Seite mit Glas abgedeckten Kathoden waren 
Eisenscheiben mit einem Durchmesser von 3 em und einer Dicke von 0,15 cm. 
Der Abstand a = 1 cm zwischen den beiden Scheiben, die parallel gegenüber 
lagen, wurde während der Versuche konstant gehalten. Die Dicke des Dunkel. 
raumes, die man durch ein ebenes Fenster aus Tempax, das seitlich an der 
Röhre .angebracht war, beobachten konnte, wurde mit Hilfe eines Katheto- 
Während der Messungen war die Entladungsröhre mit 
einer in flüssigen Stickstoff getauchten Falle in Verbindung. 
wurde so eingestellt, daß die an das Rohr angelegte Spannung praktisch dem 
Kathodenfall gleich war. 

Unter Benutzung der Beziehung (5) und der Versuchsergebnisse für j/p*, 
pd und V. wurde der mittlere Umladungsquerschnitt im Dunkelraum be- 
rechnet, eine Größe, deren Werte man in Reihe 6 der Tab. 1 verfolgen kann. 
Diese steht mit der mit Hilfe des Ausdrucks (8) berechneten und in Reihe 7 
befindlichen Mittelgeschwindigkeit v,, der positiven Ionen im Dunkelraum 
in Beziehung. Man ersieht daraus, wie wir auch weiter unten feststellen werden, 
daß der in obiger Weise gefundene Umladungsquerschnitt in guter Überein- 
stimmung mit den durch (2) berechneten Werten ist (für Argon A = 14,4: 
103 em, B= 1,6 - 10-8 em). 


Tabelle 1 
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17 | 0,18 
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Die nach (7) berechnete und in Reihe 8 wiedergegebene Konzentration der 
positiven Ionen an der Kathode n, ist von der Größenordnung 10%— 1019 cm-. 

In Reihe 9 sind die mit Hilfe von (9) berechneten Werte von N, eingetragen; 
sie liegen, wie man sieht, zwischen einigen Einheiten und einigen Zehnern. 

Was N, betrifft (Reihe 10, Ausdruck (10)), ist dieses im allgemeinen klein, 
genauso wie N, (Reihe 11 Ausdruck (11)). Diese Feststellung ist im Einklang 
mit der Annahme auf Grund deren die Poissonsche Gleichung angeschrieben 
wurde (3) und zwar, daß der elektrische Strom im Dunkelraum praktisch 
nur von den positiven Ionen, die in dem negativen Glimmlicht gebildet wurden, 
getragen wird. Wir werden sehen, daß diese Behauptung in guter Näherung 
für die Kathode der anormalen Glimmentladung im allgemeinen gültig ist. 
Die Werte des Querschnittes Q, sind aus ?!) entnommen worden. In Reihe 12 
ist die nach (12) berechnete von einem primären Elektron und seinen im 
Dunkelraum gebildeten Nachfolgern an das negative Glimmlicht übertragene 
Energie 7 eingetragen. Wie man sieht, liegt 7 unter den der Tab. 1 ent- 
sprechenden Bedingungen, ungefähr zwischen e V, und 2e V,. 

Die mittlere Energie V,, der Elektronen, die aus dem Dunkelraum ins 
negative Glimmlicht eintreten, ist nach (13) berechnet und in Reihe 13 ein- 
getragen worden. In Reihe 14 findet man die aus den in ®) angegebenen 
interpolierten Werte von 2e=e(V „). Mit Hilfe dieser Angaben und des Aus- 
druckes (14) wurde die Zahl N, der Ionisierungen (per Primärelektron) im 
negativen Glimmlicht berechnet (Reihe 15). 

Unter Benutzung des Ausdrucks (15) konnte man das Verhältnis (j,/)), 
berechnen. Dieses stellt die relative Größe des Elektronen-Stromes an der 
Grenze zwischen dem Dunkelraum und dem negativen Glimmlicht dar 
(Reihe 16). Wie ersichtlich, sind die Werte im allgemeinen klein. 

Um einen Anhalt über die relative Anzahl der Ionisierungen im Dunkel- 
raum und jenen im negativen Glimmlicht zu gewinnen, ist in Reihe 17 das 
Verhältnis N,/N, eingetragen, dessen Wert, wie man sieht, zwischen einigen 
Einheiten und einigen Zehnern liegt. 

Reihe 18 enthält den Bruchteil 7;,, von 7, der für die Auslösung der 
N,-Ionisierungen verbraucht wird; der Rest der Energie wird wahrscheinlich 
praktisch nur für die Anregungen im ee Glimmlicht verbraucht. 

In Reihe 19 sind die Werte der Größe (N, + N,) y; wiedergegeben, die den 
Beitrag des Jonenaufpralls zur Elektronenemission aus der Kathode angeben. 
Natürlich hängt diese Schätzung stark von der Genauigkeit der Bestimmung 
von y, ab. Wir schreiben, den jüngsten Messungen gemäß, y, den Wert 
2.10?) zu, ein Wert, der beim Aufprall von Argonionen auf eine mit adsor- 
bierter Argonschicht bedeckten Platinoberfläche erhalten wurde. y, bleibt 
innerhalb der Grenzen von ungefähr 5 bis 60 eV [Beziehung (8’’)] der Energien 
der positiven Ionen an der Kathode beinahe konstant; dieses Gebiet entspricht 
den Bedingungen von Tab. 1. Desgleichen gründet sich obige Schätzung auf 
die spezielle Geometrie der Kathode, die den Verlust von Ionen durch Re- 


2) M. Knoll, F. Ollendorff u. R. Rompe, Gasentladungs-Tabellen. Verlag 
J. Springer, Berlin 1935. 
2) A.v. Engel u. M. Steenbeck, Elektrische Gasentladungen. LI. Verlag J. Springer, 
sittle "hd. yvsik, Bd. XXI, S. 574 (1956). 
) P. F. Little, Handbuch d. Physik, Bd. XXI 574 (1956) ae n weh. 
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kombination an den Wanden verhindert. (Die Rekombination im Raume 
kann man vernachlässigen.) 

Nehmen wir an, daß die Elektronenemission aus der Kathode praktisch 
nur durch den Aufprall der Ionen und durch Photoeffekt ausgelöst wird, 
dann erlaubt die Beziehung (17’) die ungefähre Berechnung der Zahl der 
aktiven Photonen, die per Primärelektron gebildet werden, ein Wert, der in 
Tab. 1, Reihe 20 eingetragen ist; auch hier hängt die Genauigkeit der Ab- 
schätzung von N, von der Genauigkeit der Bestimmung von y, ab; dieses 
wurde hier (nach jüngsten Ergebnissen) gleich 10° angenommen); 
dieser Wert von y, entspricht der photoelektrischen Ausbeute an einer Platin- 
oberfläche beim Auffallen von Photonen, deren Energie der Resonanzlinie des 
Argons (A = 1066 A) entspricht. In großen Zügen folgt somit, daß N, von 
gleicher Größenordnung wie N, + N,, aber etwas größer ist. 

Aus dem über Tab. 1 Gesagten, kann man folgern, daß die Mehrzahl der 
Ionisierungen und Anregungen, die die Entladung unterhalten, sich in dem 
durch die Verschmelzung der beiden negativen Glimmlichter entstandenen 


Lichtgebilde abspielen. Dieses Gebilde, das dem Doppelkathodeneffekt 


charakteristisch ist, stellt ein äquipoten- 
550 Volt tielles Plasma dar, was man nach der 
Nullstrom-Methode mit Hilfe einer Lang- 
muir-Sonde feststellen kann). In Abb. | 

500 geben wir als Beispiel den von uns in 

Argon erhaltenen Verlauf des Potentials 

200 V,, dem man eine zylindrische Sonde auf- 

4 ladet, in den Raum zwischen den beiden 

100 } 3 ; Kathoden, wenn der Sondenstrom gleich 

Null ist; in Abb. 2 ist das Ergebnis einer 


> 


rad unter den gleichen Bedingungen bei ein- 


messung wiedergegeben. Wir konnten 


facher Kathode durchgeführten Sonden- 

goer feststellen, daB die Energieverteilung der 

} H Elektronen im gemeinsamen negativen 
200 Glimmlicht der beiden Kathoden eine 

| ! maxwellsche ist. So fanden wir z. B. 

100} für V, = 550 Volt, p=0,2mm Hg und 

a = 1cm, mit Hilfe einer ebenen Sonde 

a in der Mitte des Plasmas eine Elektronen- 


2 a temperatur 7’, = 8eV; unter gleichen Be- 
Abb. 2 dingungen ist die mit der üblichen Be- 

Abb. 1 u. 2. Das Potential Vo der ziehung errechnete Elektronenkonzen- 
Sonde, gegenüber der Anode, als tration, also auch die der positiven Ionen 
Funktion des Abstandes von der . N Mitte d Plassshs — 39. 
Kathode, wenn der Strom an der der 4 it wf es 
zylindrischen Langmuir-Sonde 10! cm, während die mit Hilfe des 
gleich null ist. Der Rand des nega- Ausdrucks (7’) errechnete Konzentration 
tiven Lichtes ist gestrichelt angegeben der positiven Ionen an der Kathode 

ungefähr dreimal kleiner ist. 

Ein anderer von uns unternommener Versuch, aus dem man viele Schlüsse 
ziehen kann, besteht im Studium der Ausbildung der Entladung mit Doppel- 


2) E. Badareu u. L. Popescu, J. Electr. a Control 4, 503 (1958). 
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kathode C,—C,, wenn man anfangs an die Kathode C, eine der Anoden- 
spannung V, gleiche Spannung V¢, anlegt und sie sodann allmählich bis zur 
an die Kathode C, angelegte Spannung V«, verkleinert. Abb. 3 zeigt den 
Verlauf der Stromdichten in Argon und der Dunkelraumdicken d, und d, 
vor den Kathoden C, bzw. C,, für V, = 450 Volt, Vo, = 0, p = 0,35 mm Hg 
und a = 1 cm; einen ähnlichen Ver- 


lauf konnten wir auch in Luft er- 03 11 
zielen Wie ersichtlich, beobachtet 
man für Vc,=V,, an C, einen am 


negativen Strom, was normal ist, 
da ja C, in diesem Falle eine zweite 
Anode ist. Bei Verringerung von 02 N 
Vo, geht der Strom 7, durch den 
Wert null (bei Vo, = 440 bis 
435 Volt) und wird dann positiv. 
Gleichzeitig fangt sich vor C, ein 


#50 


angrenzender Dunkelraum zu bilden ar 

an. Bei weiterer Verringerung von B 

Vc, nähert sich der Strom j, all- 14 
mahlich dem Werte von j,, bis jı 44 
und jg praktisch gleich werden; 
gleichzeitig verringert sich der ove 12 


Dunkelraum vor C,, während jener 
vor C, unverändert bleibt. Selbst- 450 400 30 20 100 0 
verständlich haben wir für Vo, = Vc, -— volt 

h fy Abb. 3. Die Stromdichten j, und 7, und die 
dieser Entwicklung von einer nor- Längen der Dunkelräume d, und d, an der 
malen Kathode zu einer Doppel- Kathode C, bzw. C,, als Funktion des Poten- 
kathode kann man folgendes sagen: tials von C, gegenüber der Anode 
sobald der Kathodenfall V,— Vo, 
den Wert von ungefähr 10—15 Volt erreicht, wird C, mit Unterstützung 
der Entladung von C, zur Kathode; in den Kathodenfall bzw. in den 
vor C, sich ausbildendem Dunkelraum dringen die aus dem negativen 
Glimmlicht kommenden positiven Ionen ein. Wegen des Fehlens eines 
elektrischen Feldes im äquipotentiellen Plasma gelangen diese Ionen 
bis zum Rande des negativen Glimmlichts nur durch Diffusion. Die Energie 
der zur Kathode C, gelangenden Ionen ist innerhalb eines großen Spannungs- 
bereichs (350 > Vo, > Vc, = 0) der Energie der zur Kathode C, gelangenden 
Ienen gleich, weil die Energie der Ionen an der Kathode [Beziehung (8’”’)] 
praktisch bei C, und bei C, für ein gegebenes Vc, den gleichen Wert hat; 
infolgedessen sind auch die Ausbeuten y,; gleich; sobald Vg, den Wert 
erreicht bei dem j, =j, ist, kann man somit behaupten, daß die ganze 
den positiven Ionen zuzuschreibende Elektronenemission an beiden Kathoden 
dieselbe ist; das gleiche kann man auch bezüglich der photoelektrischen Emis- 
sion behaupten, und zwar auf Grund der Gleichheit der mittleren Raum- 
winkel, die die beiden Kathoden, vom Lichtgebilde aus gesehen, einnehmen. 
Die Tatsache, daß sich der an C, liegende Dunkelraum bei Vc, kleiner als 
10-15 Volt ausbildet, zeigt, daß dies eine Folge der Abstoßung der Mehrzahl 
der Elektronen aus dem Plasma des vor C, liegenden negativen Glimmlichtes 
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ist; je mehr Vg, abnimmt, um so mehr werden auch die schnellen Elektronen 
abgestoßen, die aus dem Dunkelraum vor C, stammen und diese sind derart 
gezwungen, die Energie über die sie verfügen, für die Vergrößerung der Zahl 
der Anregungen und Ionisierungen in den gemeinsamen negativen Glimmlicht 
zu verausgaben. Die von C, ausgesandten Elektronen spielen natürlich zu. 
gleich eine ähnliche Rolle, aber solange Vc, = Vo,, werden die beiden Elek- 
tronenflüsse, die aus entgegengesetzten Richtungen von beiden Kathoden 
injiziert werden, asymmetrisch zu den Anregungen und Ionisierungen im ge- 
meinsamen negativen Glimmlicht beitragen. Dies erklärt sich dadurch, daß 
für 200 > Vc, > Vo, die Elektronen in den beiden Kathodenfällen verschiedene 
Energien aufnehmen; für Vo, > 200 kommt noch dazu, daß die Intensität 
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Abb. 4 und 5. Die Verteilungen der relativen Intensität in dem Raum zwischen C, und C;, 
fiir zwei Linien des Argons. Der Rand des negativen Lichtes ist gestrichelt angegeben 


der beiden Elektronenflüsse verschieden ist (j,< j,). Diese Asymmetrie 
kann man in Abb. 4 und 5 verfolgen, in denen, in logarithmischen Einheiten, 
die Verteilung der relativen Intensitäten im Raum zwischen C, und (, für 
zwei Linien des Argons, wiedergegeben ist: (Ionenlinie A = 4764,89 A und 
Atomlinie 2 = 6752,83 ie 
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Wie im weiter oben erörterten 
Falle wurden an C, verschiedene 
Spannungen angelegt, die zwischen 
der Anoden- und der Kathodenspan- 
nung lagen. Man sieht also, wie bei 
Verringerung von Vc, bis Vo, die 
Verteilung der Anregungen und 
wahrscheinlich auch jene der Ioni- 
sierungen zwischen den beiden Rän- 
dern des negativen Lichtes allmäh- 
lich zu einem Ausgleich strebt. 

Um das Bild der Vorgänge in 
den Kathodenteilen einer Glimm- 
entladung zu vervollständigen, wer- 
den wir noch bei der Entladung 
mit einer einzigen Kathode ver- 
weilen. Es wird daraus folgen, daß 
das über die Doppelkathode Ge- 
sagte, sowie auch die in Abschnitt 2 
erhaltenen Beziehungen in guter 
Näherung auch für die gewöhnliche 
jedoch anormale Entladung gelten. 
Um dies zu zeigen, haben wir wie 
im Falle der Doppelkathode die 
vier Größen: V., j, p und d in 
Argon gemessen und mit Hilfe der 
Beziehungen aus Abschnitt 2 die 
Tab. 2, ähnlich mit Tab. 1 aufge- 
stellt. Um die Größe der Abwei- 
chung von den Ähnlichkeitsgesetzen 
ersichtlich zu machen, sind in Abb. 6 
und 7 j/p® bzw. pd beide als Funk- 
tion von V.. wiedergegeben; was 
(v,,) betrifft, sind diese 
Werte zusammen mit jenen, die bei 
der Doppelkathode erhalten wurden, 
in Abb. 8 durch Punkte dargestellt; 
an gleicher Stelle sind durch Kreuze 
alle durch die Methode des Ab- 
bremsens des monokinetischen 
lonenbiindels erhaltenen Meßergeb- 
nisse eingezeichnet 25-30), sowie auch 


2 J.B. Hasted, Proc. Roy. Soc. 
London (A) 205, 421 (1951). 

%) B. Ziegler, Z. Physik 136, 108 
(1953). 

27) J.A. Dillon, W. F. Sheridan, 
H. D. Edwards u. S. N. Ghosh, J. 
chem. Physics 28, 776 (1955). 

2-30) s, S. 320. 
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jene, die pos Messung der Driftgeschwindigkeit als von /p?)') 
(durch kleine Kreise eingezeichnet) erhalten worden sind. Wie man aus Abb, 8 
ersehen kann, ist die Beziehung (2) für ein großes Geschwindigkeitsintervall 
der Ionen gültig. Als Abschluß der Besprechung über Tab. 2 machen wir 
noch darauf aufmerksam, daß die Größe (N, + N,)y, bei der einfachen Ka. 
thode (Reihe 17) den Maximalwert des Beitrages des Ionenaufpralls zur 
Elektronenemission aus der Kathode darstellt. Eine genauere Besprechung 


2 
4 pd mm 
pa mm Ng 
2) mm 
q A-Fe 
03 
‘| 089mmttg 
02 \ + 
0 
500 Wvot 1000 0 500 1000 Volt 
5 Er Abb.6. Die reduzierte Dichte des Abb. 7. Der reduzierte Dunkelraum 
ae Ba Stromes der Glimmentladung in an der einfachen eisernen Kathode 
ot Argon an der einfachen eisernen der Glimmentladung in Argon als 
jon 3 Kathode als Funktion des Ka- Funktion des Kathodenfalls. Die älteren 
: thodenfalls. Die älteren Daten Daten von A. Güntherschulze 


von A. Giintherschulze(siehe’)) 


(siehe 7)) sind gestrichelt angegeben 


> 
waa 


sind gestrichelt angegeben 


der Vorgänge an der einfachen Kathode kann nur an Hand von Beziehung (17) 
durchgeführt werden; dies ist jedoch vorläufig mit großen Schwierigkeiten 
verbunden, um so mehr als die ‚geometrischen‘ Faktoren f,. f,, f,,, nicht genü- 
gend genau bekannt sind. 


Wir gehen jetzt zur Beschreibung und Auslegung der spektroskopischen 
Beobachtungen über, die von uns im sichtbaren Gebiet in Argon ausgeführt 
wurden. 


Die Aufnahme des Spektrums wurde auf Kodakplatten mit einem Spek- 
trographen kleiner Dispersion für das sichtbare Gebiet (Zeiß, Modell Löwe) 
gemacht. Der Zweck der Aufnahmen war, den Verlauf der Intensität der 
Spektrallinien in dem kathodischen Teilen der anomalen Glimmentladung 


2) N. V. Fedorenko, V. V. Afrosimov u. D.M. Kaminker, J. techn. Physik % 
1929 (1956). 

2) H. B. Gilbody u. J. B. Hasted, Proc. Roy. Soc. London (A) 288, 334 (1957). 

%) J. Flacs u. E. Soloviev, J. techn. Physik 28, 599 (1958). 
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mit einfacher und doppelter Kathode zu verfolgen. Die Entladung zusammen 
mit den beiden Kathoden (im Falle der Doppelkathode) wurde auf den Spalt 
des Spektrographen so projiziert, daß der Spalt in die Achse des Bildes der 
Kathoden zu liegen kam und die Größe des Bildes nicht die Größe des Spaltes 
überschritt. Bei dem relativ niedrigen Arbeitsdruck (0,35 Torr) ist Selbst- 
absorption unbedeutend. 

Unter diesen Bedingungen VQ, 10%m 
haben wir die Intensität der % T 
Linien mit der Intensität 
des Untergrundes verglichen. 
In den Abbildungen haben | 
wir relative Intensitäten m Argon. 
in logarithmischen Einheiten > 
wiedergegeben. Die Inten- 
sititen wurden mit einem ‘ 
ZeiB II, gemessen. Die ; x 
großen Unterschiede zwi- in 
schen den Schwärzungen bei 

verschiedenen Versuchsbe- 
dingungen erhaltene Schwär- 0, 5 r 7 
zungen haben eine sorg- log (vom/s) 
fältige Wahl der Expositions- Abb.8. Der Umladungsquerschnitt in Argon als 
zeiten gefordert, um noch Funktion der lonengeschwindigkeit 
genügend genau meßbare 

Schwärzungen zu erhalten. In dem von uns untersuchten Spektralgebiet, 
das zwischen 4131,73 A und 6965,43 A lag, konnte man folgende Verteilung 
der Spektrallinie des Argons feststellen: 30 Atomlinien mit einer mittleren 
Anregungsspannung von 14,5 Volt und 21 Linien des einmal inonisierten- 
Argons (darunter 16 Linien mit einer Anregungsspannung von etwa 19 Volt, 
4Linien mit einer solchen von etwa 21 Volt und 1 Linie, die etwa 23 Volt 
entspricht). Wir bemerken, daß ein großer Teil der Linien der Ionen im nega- 
tiven Glimmlicht ungefähr gleich oder stärker erscheinen als die Atomlinien, 
in Übereinstimmung mit älteren Bestimmungen des Anregungsquerschnittes 
der Argonlinien*)%?). Ionenlinien, die größeren Anregungsspannungen (etwa 
37 Volt) entsprechen, wurden in einer Entladung mit zylindrischer Hohl- 
kathode®®) in Neon beobachtet. Diese Beobachtungen geben uns einige An- 
halte über die Elektronenenergie im negativen Glimmlicht. 


Die Intensitäten der Argonlinien, von denen weiter oben die Rede war, 
deren Verlauf besonders charakteristisch ist, sind in Abb. 4 und 5 dargestellt. 
Auch -die H,-Linie des Wasserstoffs, die zufällig als Verunreinigung im ver- 
wendeten Argon erschien, zeigt eine interessante Intensitätsverteilung, 
die in Abb. 9 wiedergegeben ist. Wie ersichtlich, befindet sich das Maximum 
der Kurve im Dunkelraum und nicht an der Grenze des negativen Glimm- 


1) O. Fischer, Z. Physik 86, 646 (1933). 


8) 1. P. Bogdanova u. Tschen Gu-Thek, Opt. i. Spektr. 2, 681 (1957). 
7. Folge, Bd. 5 
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dadurch erklären, daß die von der Kathode ausgehenden Elektronen i im star. 
ken elektrischen Feld des Dunkelraumes eine allmähliche Beschleunigung er. 
fahren und sodann das Maximum der Anregungsfunktion der H,-Linie über. 
schreiten. Es ist jedoch kein Beitrag zur Anregung dieser Linie von seiten 
der zahlreichen aber langsamen Elektronen, die sich im Plasma bilden, zu 
beobachten; diese Feststellung macht uns glauben, daß ihre Temperatur, die, 

wie wir sahen, etwa 8eV beträgt, nicht 

ausreichend ist um das betreffende Niveau 

uk, mit einer genügend großen Wahrscheinlich. 

en ) keit anzuregen. Der monotone Verlauf 

der Intensität im negativen Glimmlicht 

ist eine Folge der Abbremsung und des 

allmählichen Verschwindens der Primär. 

elektronen. Ein ganz ähnliches Verhalten 

zeigt auch die Verteilung der Intensität 

ar der H;- und H,-Linien®). Unter Be. 

' nutzung dieses Verhaltens sowie des Aus. 

- druckes (6’) und unter Vernachlässigung 

2250 | G der Diffusion der angeregten Atome kann 

1 man Anhalte über den Verlauf der An- 

- regungsfunktion für die betreffende Linie 

} für zwischen 0 und eV, liegende Elek- 

: G tronenenergien erhalten. Selbstverstand- 

- | lich muß man zu diesem Zwecke das 


Spektrogramm einer genügend anomalen 
Entladung aufnehmen, um praktisch nur 
G die Anregung seitens der Primärelektronen 

zu erhalten und nicht auch jene seitens 
V,=430 Volt ihrer Nachfolger aus dem Dunkelraum. 


verlauf, der jenem aus Abb. 4 fiir die 
G Linie A= 4764,89A des Argons für 
r r . . . 
4500! (Ve, = 0 V) ähnlich ist. Die Anregungs- 
2 ' funktion dieser Linie hat wahrscheinlich 
shoes fit keine zu hohen Werte, so daB die wenigen 
Abb. 9. Die Verteilung der relativen Primärelektronen kaum zur Anregung bei- 
Intensität in dem Raum zwischen C, t kö d sich desbeib u 
und C, für die H,-Linie des Wasser- "Tagen Sonnen und sic 
stoffes unter den Bedingungen der Intensitätsmaximum im Dunkelraum b - 
Entladung in Argon. Der Rand des den kann. Andererseits ist es aber wahr- 
negativen Glimmlichtes ist gestrichelt scheinlich, daß das Maximum der Funktion 
e angegeben einem bedeutenden Teil der im negativen 
Glimmlicht gebildeten Elektronen zu- 
_ gänglich ist, die dank ihrer großen Anzahl einen großen Beitrag zu den An- 
regungen im Plasma liefern. 


; Andere Linien zeigen einen Intensitäts- 


%) Von diesem Standpunkt aus sind auch einige jüngst erhaltenen Ergebnisse*) 
über die Anregung der Wasserstofflinien in den Kathodenteilen der Glimmentladung be- 
merkenswert. 

%) F. Schneider, Ann. Physik 18, 297 (1956), 
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star. Die allmähliche Abnahme der Intensität der Linie vom negativen Glimm- 
g er.| licht bis zur Kathode kann man nur durch eine Diffusion der angeregten Atome 
iber.| aus dem Plasma erklären; weil es den langsamen Elektronen aus dem 
eiten | Plasma vollkommen unmöglich ist, einen größeren Weg gegen das Feld zu 
1, gy} durchlaufen. Wir können dabei versuchen die mittlere Lebensdauer 7 des 
die, | angeregten Atoms auf dem betreffenden Niveau zu schätzen; dies folgt aus 
icht | der kinetischen Theorie 

veau t= zy" 
lich- 
lauf | wo & die mittlere V Verschiebung des Teilchens nach einer Koordinate und D 
licht | der durch RAR 

miir- 

ten | gegebene Diffusionskoeffizient ist. (C = thermische Mittelgeschwindigkeit 
sität | und A, = mittlere freie Weglänge.) Die auf solche Art geschätzte mittlere 
Be. | Lebensdauer 7 für des oberen Niveaus der betrachteten Argonlinie beträgt 
Aus. | ungefähr 10-7 Sek. Für die anderen sichtbaren Linien des Argons erhielten 
zung wir für 7 Werte gleicher Größenordnung. 

ann Schließlich haben Linien, wie z. B. A = 6752,83 A des Argons bei Vc, = 
An. | 450 Volt, Intensitätsverteilungen wie jene aus Abb. 5. Die Erklärung dieses 
‚inie | Verlaufs fußt auf der Deutung der beiden von uns oben analysierten Fälle: 
‘ek. | das Maximum im Dunkelraum ist den Primärelektronen zuzuschreiben (einen 
ind. | gewissen Beitrag steuern sicherlich auch ihre Nachfolger aus dem Dunkel- 
das | raum bei), während das Maximum im negativen Glimmlicht praktisch nur 
alen | das Ergebnis der Wirkung der zahlreichen und relativ langsamen Elektronen 
nur | aus dem Plasma ist; auch hier kann man die Diffusion der angeregten Atome 
onen | aus dem Plasma in den Dunkelraum beobachten. 


tens 
vum. 4. Weitere Ergebnisse. Einige Einschränkungen 
täts- Obige Betrachtungen bestätigen die Annahme, daß die von der Kathode 


die | ausgesandten Primärelektronen und die von ihnen im Kathodendunkelraum 
für | gebildeten Sekundärelektronen ihre Energie praktisch nur im negativen 
ings- | Glimmlicht verausgaben. 
lich Auf Grund dieser Annahme konnte schon vor seit langem gezeigt werden, 
igen | daß die Länge des negativen (anomalen) Glimmlichtes in verschiedenen Gasen 
'bei- | und bei verschiedenen Kathodenfällen gleich dem Relaxationsweg jener 
kein | Elektronen ist, die eine dem betreffenden Kathodenfall entsprechende Energie 
| bil- | besitzen). Die ersten Berechnungen, die in dieser Hinsicht®) gemacht 
ahr- | wurden, haben diesen Mechanismus bei der anomalen Glimmentladung be- 
‘tion | stätigt; es ist aber interessant zu bemerken, daß dieser Standpunkt bis jetzt 
iven | noch nicht allgemein anerkannt ist, ein entgegengesetzter diesbezüglicher 
zu- | Standpunkt wird in ?) vertreten. Eine andere Grundannahme der obigen Ab- 
An-| schnitte besteht darin, daß der Umladungsprozeß die Kinetik der positiven 
Ionen im eigenen Gas beherrscht; in dieser Hinsicht scheint es keine Un- 
einigkeit mehr zu geben!!),8), ®), 37), 14), 15), 10), Auf dieser Grundlage erhält 
ges) | Man leicht die notwendigen analytischen Ausdrücke der Driftgeschwindigkeit 


%) A.K. Brewer u. J. W. Westhaver, J. Appl. Phys. 8, 779 (1937). sty 2 
#®) G. H. Wannier, Physic. Rev. 96, 831 (1954). ae 
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als Funktion von E/p; diese Ausdriicke wurden ups bestätigt und | (pd ,, 
führten auf Grund der Versuchsergebnisse von Beweglichkeitsmessungen zu | em) 
richtigen Werten der Umladungsquerschnitte®) 17). Wa 
Die Entstehung der Raumladung gegenüber der Kathode ist ein direktes | Gleich 
Ergebnis der Bewegung der Ionen und in letzter Analyse des Umladungs- | ständi; 
prozesses. Für relativ große E/p bekommt man für die Driftgesehwindigkeit 
den Ausdruck (1), der die Form des Feldes E ~ x”): festlegt. Es ist interessant 
zu bemerken, daß auch das Modell des elastischen Stoßes mit isotroper 2 
Streuung nach jedem Stoß zu einem Ausdruck von v derselben Form (v~(E/p)'!) Die 
führt3® 4%), woraus sich die gleiche Form des Feldes vor der Kathode er. | Lang! 
gibt; doch der absolute Wert der experimentell bestimmten Feldstärke konnte | Stöße 
auf diese Weise nicht erhalten werden®!)42). Der richtige Wert folgte nur aus | in Er 
der Inbetrachtziehung des Umladungsprozesses der zu den Ausdrücken (1) | (bzw. | 


bzw. (4) führt; Q, ist der entsprechende Umladungsquerschnitt!P). muß. 
Was die Einschränkung der angegebenen Beziehungen betrifft, so müssen 

wir darauf aufmerksam machen, daß man bei sehr großen pd-Werten, wie sie Au: 

ungefähr der normalen Glimmentladung entsprechen, sowie auch bei sehr | Mecha 

kleinen p d-Werten gewisse Korrektionen anbringen muß. der Ke 
Tatsächlich kann man für große pd die Ionisierung im Dunkelraum nicht a) 


mehr vernachlässigen, in diesem Falle müssen wir die exponentielle Ent- | Glimm 
wicklung der von der Kathode ausgehenden Elektronenlavine berücksichtigen, | relativ 
und der die Feldstärke bestimmende Ausdruck muß, wie in ®) gezeigt wurde, | ist, die 


folgendermaßen vervollständigt werden: Die Se 
darin ist y das Verhältnis zwischen dem Elektronen- und dem Ionenstrom an ry 

der Kathode, und « der mittlere Townsendsche Koeffizient im Dunkel- ant 
raum. Man kann zeigen, daß der Ausdruck (4’’’) die ganze Spanne der ex- ung 
perimentell bestimmten Felder umfaßt. So ergibt z. B. (4”’) für die normale b) 
Entladung in Argon bei pd = 0,33 mm Hg - em, mit «/p = 9 (em-mm Hg) ay 

und y= 4y,=8-10-**) eine praktisch lineare Abhängigkeit der Feld- | „ge | 
stärke E von x. Im genügend anomalen Gebiet, wo der Beitrag des zweiten en: 
Gliedes der quadratischen Klammern vernachlässigbar ist, gilt dies nicht mehr; sdeotin 

in diesem Fall vereinfacht sich (4’’’) zu (4). Auslös 
Mehr als dies; wenn man das Feld (4’’’) integriert, erhält man einen all- en 
gemeinen Ausdruck zwischen V ., j/p?, pd, &/p, Q, und m, der in jedem Gebiet und 2: 

der Glimmentladung giiltig ist (auch ist dieser Ausdruck durch (5) unter den oe 
Bedingungen der anormalen Entladung gut angenähert). So, z. B. unter Be- sleiche 
nutzung des von (4’”) gegebenen Feldes, kann man leicht das gleichzeitige StoBes 
Vertragen unter den Bedingungen der normalen Entladung in Argon fol- Veriine 
gender Werte überprüfen: (V = 165 Volt, (j/p?), = 160- 10-8 A (cm -mm Hel 

38) G. H. Wannier, Physic. Rev. 83, 281 (1951). | positiv 

3) G.H. Wannier, Physic. Rev. 87, 795 (1952). N sobald 

40) G. H. Wannier, Bell. Syst. Tech. Jour. 32, 170 (1953). im Du 


41) R. Warren, Physic. Rev. 98, 1650 (1955). Tor. 1 

42) R. Warren, Physic. Rev. 98, 1658 (1955). ha a 

ou gr Dieser Wert folgt aus der Angabe der ersten senkrechten Reihe der Tab. 2 für He Ka 
(N, + No) Umlad 


FR 


2 für 
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(pd „= 0,33 mm Hg-em, (a/p), = 9 (mm Hg- em)-!, Q,= 170 (mm Hg - 
m), y = 8,4 - 102. 

Was die sehr kleinen Werte von p d betrifft, konnte von der gaskinetischen 
Gleichung ausgehend gezeigt werden?®), daß man die Beziehung (5) voll- 


Diese Beziehung erweist sich bei der Untersuchung des Ionenstromes an 
Langmuir-Sonden mit Inbetrachtziehung der von Umladung begleiteten 
Stöße im Dunkelraum als außerordentlich nützlich. Die Beziehung (5’) ist 
ein Ergebnis der Inbetrachtziehung der starken Inhomogenität des Feldes 
(bzw. der Kleinheit von Q,- pd), so daß der Ausdruck (1) korrigiert werden 
muß. 

5. Schlußbemerkungen 


Aus den oben Erörterten ergeben sich folgende allgemeine Züge des 
Mechanismus der Vorgänge an der Doppelkathode und im allgemeinen an 
der Kathode der anomalen Glimmentladung. 

a) Der größte Teil der Anregungen und Ionisierungen geht im negativen 
Glimmlicht vor sich; dieses Gebilde ist ein äquipotentielles Plasma, das von 
relativ wenigen aber energiereichen Elektronen durchkreuzt und unterhalten 
ist, die von der Kathode oder aus dem Kathoden-Dunkelraum abstammen . 
Die Sekundärelektronen, die im Plasma gebildet werden, bewegen sich chao- 
tisch und sind für die Mehrzahl der im negativen Glimmlicht gebildeten An- 
regungen verantwortlich, sie können leicht durch eine dem Plasma gegen - 
über mit mehr als ungefähr 15 Volt negativierte Elektrode abgestoßen werden; 
zugleich bildet sich an der betreffenden Elektrode ein Dunkelraum; dieser 
spielt die Rolle eines Beschleunigers für die Primärelektronen und ihre im 
Dunkelraum gebildeten Nachfolger, der sie mit der zur Erhaltung der Ent- 
ladung notwendigen Energie versieht. 

b) Bei Benutzung der letzten Angaben über die Elektronenemission unter 
dem Aufprall von positiven Ionen auf metallische Flächen, die in das betref- 
fende Gas getaucht sind, reicht die Zahl der per Primärelektron gebildeten 
Ionen nicht aus, um aus der Kathode ein neues Elektron auszulösen. In Über- 
einstimmung mit jüngsten Forschungen?), müssen wir auch andere bei der 
Auslösung von Elektronen aus der Kathode wirksame Faktoren in Betracht 
ziehen, und zwar in erster Reihe Photonen mit einen Energie zwischen 8 
und 25 eV. 

c) Der Grundprozeß bei der Wechselwirkung zwischen Ion und Atom 
gleicher Art ist die reine Umladung, durch die zwischen den Partnern des 
StoBes ein Ladungstausch vor sich geht, ohne jedoch von einer merklichen 
Veränderung der Geschwindigkeit der Partner begleitet zu sein. 

Im äquipotentiellen Raum des negativen Glimmlichtes bewegen sich die 
positiven Ionen in ,,isotroper Staffette‘‘ und, gemäß der Diffusionsgesetze, 
sobald sie den Rand des negativen Glimmlichtes erreicht haben, unterliegen 
im Dunkelraum sie kräftigen Feldern der Größenordnung 10?—10% Volt/em 
Torr, und ihre Bewegung verwandelt sich in eine ,,anisotrope Stafette“. An 
die Kathode gelangen sie mit einer etwa der Hälfte der dem letzten freien 
Umladungsweg entsprechenden Energie. 
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d) Da bei der anomalen Glimmentladung die Zahl dee seninlabiiage im 
Dunkelraum unbeträchtlich ist, folgt, daß in diesem Gebiet der Strom praktisch 
nur von positiven Ionen getragen wird. Weil aber die Geschwindigkeit der 
Ionen allmählich in Richtung auf die Kathode zu wegen des immer stärker 
werdenden Feldes wächst, so folgt weiter, daß zugleich die Konzentration 
der positiven Ionen in Richtung auf die Kathode zu abnehmen muß. Ein 
Anwachsen der Konzentration in Richtung auf die Kathode zu, wie es in 
neueren Arbeiten angegeben wird") *), ist unter den Bedingungen der gewöhn- 
lichen anomalen Glimmtentladung schwer zu erklären. Im Gebiet der normalen 
Glimmentladung kann wegen der zahlreichen Ionisierungen, die im Dunkel- 
raum stattfinden, die Konzentration der positiven Ionen in größerer Nähe vor 
der Kathode beträchtlich zunehmen, so daß ein quasilinearer Verlauf des 
Feldes oder sogar eine ,,nach oben‘‘ konkave Kurve möglich ist. Dies kann 
man auch aus der Beziehung (4’’’) folgern. 

e) Die im Vergleich zur einfachen Kathode Stromvervielfachung beim 
Eintritt des Doppelkathodeneffektes bei gleichem Druck und Kathodenfall 
ist sowohl eine Folge der Vergrößerung der Ionengeschwindigkeit an der 
Kathcde (in dem eingeengten Dunkelraum erleiden die Ionen weniger Um- 
ladungszusammenstöße) als auch der Zunahme ihrer Konzentration (beinahe 
die gleiche Energie der Elektronen wird jetzt einem viel kleineren Raum 
abgegeben). 

f) Schließlich bemerken wir noch, daß Prozesse wie die lonisierung des 
Gases durch beschleunigte positive Ionen oder die Elektronenemission aus 
der Kathode unter dem Aufprall energiereicher neutraler Atome, die aus 
Umladungen stammen, im allgemeinen nur in sehr großen Kathodenfällen 
und bei genügend kleinen Drucken vorkommen; unter den Bedingungen der 
in vorliegender Arbeit analysierten Entladung in Argon spielen derartige Vor- 
gänge infolge der 60 eV nicht überschreitenden Energie der positiven Ionen 
eine vollkommen vernachlässigbare Rolle. 


44) J. Wilhelm, Ann. Physik 17, 225 (1956). abe 

43) (+. Brederlow, Ann. Physik 1, 359 (1958). gee 
Mer 


ı ay Bukarest, Physikalisches Institut der Akademie der RVR. 


Ss Hbiat rm Bei der Redaktion eingegangen am 22. Juni 1959. 
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Zum Ubergangseffekt der neutronenauslösenden 


Komponente der kosmischen Strahlung in Blei — 


Von I. Hauser 


Inhaltsübersicht 


Mit borbeladenen Agfa K3-Kernemulsionen wurde die Neutronenerzeu- 
gung durch die kosmische Strahlung in Blei bis zu 2,4 cm Dicke gemessen. 
Die mit der verwendeten Absorberanordnung erhaltenen Ergebnisse zeigen 
einen glatten Verlauf der integralen Neutronenerzeugungskurve entsprechend 
einer exponentiellen Abnahme der erzeugenden Strahlung. Ein Hinweis auf 
ein früher von anderer Seite angegebenes Maximum zwischen 1—2 cm Blei 


konnte nicht festgestellt werden. 


Einleitung 


Verschiedene Autoren!-7) haben hinter festen Absorbern (Pb, C, Fe, Sn 
und Al) Intensitätsmessungen der sternauslösenden Komponente der kosmi- 
schen Strahlung in Kernemulsionen in Abhängigkeit von der Absorberdicke 
in Höhen von 3000—4000 m durchgeführt. Sie geben an, daß die Intensität J 
mit zunehmender Absorberdicke x nicht, wie zu erwarten war, nach einem 
Exponentialgesetz 

I=I, (1) 


abfällt (A = Absorptionslänge der Strahlung), sondern 1 oder 2 Maxima durch- 
läuft. 

Die Lage des ersten Maximums hinter Bleiabsorbern schwankt in den 
verschiedenen Arbeiten zwischen 0,5 und 2cm. Die Zunahme der Intensität 
gegenüber dem Wert bei 0 cm Pb liegt zwischen 0,5 und 2cm. Die Zunahme 
der Intensität gegenüber dem Wert bei 0 cm Pb liegt zwischen 5 und 35% 
Ein zweites Maximum ist bei ungefähr 20 cm Pb angegeben worden. 

Da heute als gesichert angenommen werden kann, daß die Neutronen der 
kosmischen Strahlung vorwiegend aus Kernzertrümmerungen herrühren, 

') G. Bernardini, G. Cortini u. A. Manfredini, Physic. Rev. 76, 1792 (1949). 

*) L.Tomäskovä, Csechosl. Journ. Phys. 2, 102 (1953). 

3) E. Schopper, K. H. Höcker u. G. Kuhn, Physic. Rev. 82, 444 (1951). 

‘) E. Schopper, K. H. Hicker u. E. Rössle, Z. Naturforschg. 6a, 603 (1951). 

5) L. Malaspina, M. Merlin, O. Pierucciu. A. Rostagni, Nuovo Cimento 7, 145 


6. Belliboni, L. Fabbrichesi, L.de Marco u.M. Merlin, Nuovo Cimento 8, 374 


%. Rössle u. E. Schopper, Z. Naturforschg. 9a, 836 (1954). ar 
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während der Beitrag der anderen Komponenten nur gering ist, haben es 
mehrere Autoren®-1!) unternommen, die Übergangseffekte durch Messungen 
der Neutronenintensität in festen Absorbern (Pb, Al und C) weiter zu unter- 
suchen. Als Nachweisinstrumente bei den Messungen in Blei wurden Zähl- 
rohranordnungen benutzt. 

Die gewonnenen Meßergebnisse stimmen nicht miteinander überein. 
Einige Resultate®)!%) weisen auf das Vorhandensein des ersten Maximums 
der Übergangskurve hin; bei anderen Messungen 9) !1)12) konnte dieses Maximum 
nicht festgestellt werden. Eine befriedigende Deutung für das Auftreten der 
Maxima in den Übergangskurven der stern- und neutronenauslösenden Kom- 
ponenten konnte bisher nicht gegeben werden. 

In der vorliegenden Arbeit wurde der Ubergangseffekt der neutronenaus- 
lösenden Strahlung in Blei mit borbeladenen Agfa K3-Emulsionen mit einem 
Borgehalt von 0,063 g/cm über die Reaktion B!® (n, x) Li’ untersucht. Es 
wurde eine integrale Intensitätsmessung durchgeführt. 


Experimentelle Anordnung 


Der benutzte Absorber (Abb. 1) bestand aus drei konzentrischen Blei- 
zylindern mit einer Höhe von 50 cm. Die Zwischenräume und der Innenraum 
waren mit Paraffin ausgefüllt. Der innere Paraffinzylinder besaß in Richtung 
der Achse eine Bohrung, in der sich der Platten- 


ee behälter, bestehend aus Messingblech mit einer 
ge: 3 5 Wandstärke von 1 mm, befand. Die Emulsionen 
vor = = (3cm x 6em x 100m) wurden vertikal in die 
% = S|: = Plattenbehälter eingelegt. 
ale Die Messungen wurden mit Bleizylindern, 
= F335 8G deren Wandstärken im Mittel 5,9, 12,1 und 
3: Sim 24,4 mm betrugen, durchgeführt. 

sy Die Dauer der Exponierung betrug 38 Tage 


in einer Höhe von 2630m über dem Meeres- 
spiegel (Lomnicky Stit). Über den Absorbern 
befand sich ein Holzdach mit einer Wandstärke 
von ungefähr 0,8 g/em?, das mit einem 1 mm 
dicken Aluminiumblech bedeckt war. 

Die Messung des Nulleffektes, d. h. der Zahl 
der Neutronen aus der Umgebung, wurde mit 
einem Aufbau vorgenommen der kein Blei enthielt. 

Vor Beginn jeder Exponierung wurde eine 
Kontrollplatte entwickelt, die sich bis zum Be- 
Abb.1. Absorberaufbau, im ginn der Versuche unter gleichen Bedingungen 
Schnitt gezeichnet befand wie die Emulsionen in den Absorbern. 


Auswertung und Korrekturen 
Die Auswertung der Emulsionen erfolgte mit einem Lumipan-Mikroskop 
von Zeiß, Jena, mit einer 450fachen Vergrößerung. Im Okular befand sich 
i 8) S.B. Treiman u. W. Fonger, Physic. Rev. 85, 364 (1952). 4 
ae 9) K. Hogrebe, Z. Naturforschg. 7a, 772 (1952). = 
2) H.M. Weiss, Z. Naturforschg. 10a, 21 (1955). 
2) H. Reich, Z. Naturforschg. 10a, 914 (1955). 
12) V. Cocconi-Tongiorgi, Physic. Rev. 76, 517 (1949). 
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I. Hauser, Kosmische Strahlung in Blei 


ein | Okularmikrometer, das bei der verwendeten Vergrößerung eine Kreis- 
fläche von 1,09 - 107 2 mm? der Emulsion begrenzte. 

In den exponierten Platten wurden alle Spuren gezählt, die folgenden Be- 
dingungen genügten: 

1. Ihre Reichweiten mußten innerhalb eines begrenzten Intervalls liegen 
(siehe unten). 

2. Sie mußten im Innern des Zylinders liegen, der durch den Kreis im 
Okular aus der Emulsion ausgeschnitten wurde oder konnten die Mantel- 
fläche dieses Zylinders schneiden (eine Berührung von außen war ausge- 
schlossen). 

3. Die Spuren mußten bezüglich Basis und Oberfläche vollständig in der 
Emulsion verlaufen. 

Schrag zur Emulsionsebene liegende Spuren wurden getrennt registriert. 
Die Länge R dieser Spuren wurde nach der Formel 

R=/VP2+ 
aus ihrer Längenprojektion / und der Höhenprojektion h berechnet. Der 
Schrumpfungscharakter S der K3-Emulsion wurde zu 2,64 + 0,04 bestimmt. 
Dieser Wert gilt fiir eine relative Luftfeuchtig- 
keit von 64°%,. } Anzahl ger Spuren 

Zur Bestimmung der Reichweiteverteilung a-u? u 
der Spuren aus der Reaktion B!® (n,a) Li? 
wurden in einer mit einer Neutronenquelle 
bestrahlten K3-Emulsion die Längen von 
660 x-Li?-Spuren gemessen, die parallel zur zs9L 
Emulsionsebene verliefen. Das Ergebnis ist in 
Abb. 2 dargestellt. Auf Grund dieser Verteilung 
wurde das Reichweiteintervall zu 4—10 u fest- — go} 
gelegt. 

Weiter wurde, um eine mögliche Fehler- 
quelle auszuschließen, die Homogenität des 50} 
Borgehaltes in den Emulsionen geprüft. Zu 
diesem Zweck wurden K3-Emulsionen mit 
einer Neutronenquelle bestrahlt. In diesen 
Emulsionen wurden an 10 verschiedenen, be- 
liebig verteilten Stellen die Häufigkeiten der i 
x-Li?-Spuren in gleichen Volumen bestimmt. 2. 

der Spuren aus der Reaktion 
Jedem Meßwert entsprachen etwa 200 Spuren. BY (n, x) Li? 
Ein Vergleich der Werte zeigte, daß die Ab- 
weichungen der Spurdichten innerhalb der erlaubten Fehlergrenzen lagen. 

Zur Berechnung der Spurdichten ist eine möglichst genaue Bestimmung des 
ausgewerteten Emulsionsvolumens notwendig. Da nicht nur von einer Platte 
zur anderen, sondern auch innerhalb einer Platte größere Dickenschwan- 
kungen auftreten können, wurden in festen Abständen Dickenmessungen an 
den Stellen der Emulsionen durchgeführt, in denen die Bestimmung der 
Spurzahlen erfolgte. Diese Messungen wurden mit einem ‚„Zeiß-Orthotest‘‘. 
Jena (Teilungswert 1 u), vorgenommen. 

Da sich die Dicke mit ihrem Feuchtigkeitsgehalt ändert, befanden sich die 
Emulsionen während dieser Messungen unter kontrollierten Feuchtigkeits- 
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7. Folge. Band 5. 1960 


a Mit Hilfe einer Eichkurve (Abb. 3), die die prozentuale Dicken- 
änderung in Abhängigkeit von der relativen Luftfeuchtigkeit angibt, wurden 


Prozentuole Abweichung 
50 
+30 
10 a 
Relative Lufffeuchte?% 


Abb. 3. Prozentuale Dickenänderung 
der Emulsion in Abhängigkeit von 
der relativen Luftfeuchte 


alle Dickenwerte auf eine Feuchtigkeit 
von 45%, umgerechnet. 

Bei der Auswertung der exponierten 
Platten wurden nur Spuren gezählt, die 
einen Neigungswinkel mit der Emulsions- 
ebene einschlossen, der nicht größer als 
70° war. Es mußte also an den Meß. 
werten eine Korrektur vorgenommen 
werden, die die nicht endenden Spuren 
berücksichtigt. 


line einfache geometrische Betrachtung?) ergibt für den Bruchteil dieser 


nicht gemessenen Spuren: 
2d 
9 
N,= ,(53- - Rsin y : (2 
N, = Zahl der hinausgehenden Spuren (unbekannt), = 
N, = Zahl der endenden (registrierten) Spuren, _ engen 
d = Dicke der Emulsion, are ai 
R = Reichweite der Spuren, 
y = größter beobachteter Winkel der Spuren. ie 
Diese Methode stellt eine Annäherung 
=--sofort nach Bestra- «dar, da es wegen des diskontinuierlichen 
Spurbildes nicht immer möglich ist zu 
a-u? Besratlung u. Ent- entscheiden, ob die Spuren in der Emul- 
L ZA wicklung 28 Tage sion enden oder nicht. 


\ - 


er 


J 4 5 6 ? 6 9 on 
Länge (Skt) 
Abb.4. Reichweiteverteilung der 
Spuren o-Li? in Abhängigkeit von 
der Lagerungszeit zwischen Bestrah- 
lung und Entwicklung 


Sämtliche Meßwerte wurden unter 
Verwendung der Formel (2) korrigiert. 
Für die Länge R wurde der Wert Tu 
(Mittelwert des Längenintervalls) einge- 
setzt. 

Da die Exponierungsdauer der Emul- 
sionen 38 Tage betrug, können nach Ab- 
lauf dieser Zeit durch allmählichen 
Schwund des latenten Bildes (Fading) 
Spuren, deren ursprüngliche Länge 4—10 u 
betrug, nicht mehr in dieses Intervall 
hineinfallen bzw. Spuren, deren ursprüng- 
liche Länge größer war, hinzukommen. 

Um eine Aussage über die Größe 
dieses Effektes machen zu können, wurden 
zwei K3-Emulsionen im Abstand von 
28 Tagen mit einer Neutronenquelle be- 
strahlt und dann gemeinsam entwickelt. 


In jeder Emulsion wurden dann die Reichweiteverteilungen von 760 will- 
kürlich ausgewählten «-Li’-Spuren gemessen. In die erhaltenen Diagramme 
sind die entsprechenden Normalverteilungen eingezeichnet worden (Abb. 4). 


13) E.M. Binggeli, Helv. phys. Acta 28, 3 (1955). - 
gg phy ( ) 
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mee I. Hauser: Kosmische Strahlung in Blei i 


en- 7 Ein Vergleich der beiden Kurven zeigt, daß ein Fadingetfekt vorhanden ist. 
den | pie Zahl der Spuren, die aus dem vorgeschriebenen Intervall herausfällt, 
keit | peträgt jedoch nicht mehr als 1,5%. 

Berücksichtigt man nun, daß das Fading nur auf ungefähr die Hälfte der 
ten | während der Gesamtbestrahlungszeit aufgezeichneten Spuren einwirkt, so 
die | kann eine Korrektur bezüglich Fading vernachlässigt werden. 

INS- Ein Teil der in den Bleizylindern erzeugten Neutronen wird beim Durch- 
als queren der Absorber im Blei reflektiert oder absorbiert. Der Anteil dieser 
leB- | Neutronen wird mit zunehmender Bleidicke anwachsen. Um den Einfluß ab- 
nen | shitzen zu können, den dieses Verhalten auf die MeBwerte ausübt, wurden 
ren folgende Messungen durchgeführt: in Radium-Beryllium- Präparat wurde 
an Stelle der Bleizylinder in den Absorber gebracht und in den Zwischen- 
sser Jriumen zwischen den Paraffinzylindern vertikal bewegt. Gleichzeitig wurde 
der Plattenbehälter mit den Emulsionen kontinuierlich gedreht. Die Be- 
(2) strahlungszeiten in den verschiedenen Lagen waren den Radien der Blei- 
wlinder proportional. Bei einer zweiten Messung befanden sich die Blei- 
ıylinder mit einer Wandstärke von 12,1 mm in dem Aufbau und die gleiche 
Messung wiederholte sich. Um das Präparat vertikal bewegen zu können, 
waren die Pb-Zylinder jeweils mit einem senkrechten Spalt von etwa 1 cm 
Breite versehen worden. Die Auswertung der Emulsionen führte zu folgendem 


Ergebnis: 
ung Tabelle 1 
hen | | 
zu Pb-Dicke | . = Ausgewertetes Spurdichte 
[mm] Spurzahl Emulsionsvolumen [mm-? h-1] 
| [mm3] 
ern 0 1550,0 + 39 0,238 + 0,30, 1447,95 + 2,60, 
ger 12,1 1419,2 + 38 0,200 + 0,3%, 1612,54 + 2,6%, 
1ge- Der Anteil der nicht in der Emulsion endenden Spuren ist bereits in den 
Werten der Tab. 1 enthalten. 
nd Die Korrekturwerte fiir 0,6 und 2,4 cm Blei wurden erhalten, indem durch 


Ab. | tie beiden Meßpunkte für 0 und 1,2 cm Pb eine Gerade hindurchgelegt wurde. 
Bei den geringen Bleidicken und der Größe 


“a des MeBfehlers ist diese lineare Näherung 10 Aorrekturfoktor 

0 1 durchaus zulässig. 06 

vall |. In Abb. 5 sind die Werte des Korrektur- 02 

ng. aktors, mit dem die beobachteten Spur- ? 
en, | Uchten multipliziert werden müssen, über —_ 


öße | et Bleidicke aufgetragen. Diese Korrektur Abb. 5. Korrekturfaktor, der die 
den | *tzt voraus, daß die Energie- und Richtungs- 
verteilungen der Neutronen aus einer Ra—Be- ~ = 
be. | Muelle mit den entsprechenden Verteilungen satis 

elt. der Verdampfungsneutronen aus Kernzertrümmerungen übereinstimmen. 
vill- Diese Voraussetzung ist innerhalb der Fehlergrenzen sowohl fiir die Energie 
als auch für die Winkelverteilung erfüllt ™) 1). 


4). 4) W. Heisenberg, Vorträge über kosmische Strahlung, Springer-Verlag, Berlin- 
böttingen-Heidelberg 1953. 


8) U. Schmidt-Rohr, Z. Naturforschg. 8a, 470 (1953). 


| 


Luftdruck von 550, 0 mm wt, —— Zur ae wurde der von 
Cocconi-Tongiorgi!?) angegebene Barometerkoeffizient 


B= — 11,2 + 0,6 %/em Hg 


In der Tab. 2 sind die beobachteten Spurzahlen fiir die verschiedenen Ab- 
sorberdicken angegeben. Der Fehler entspricht der mittleren ea 


Bleidicke beobachtete Versuchsdauer = Spurdichte 
[mm] Spurzahl | [h] [mm=3 h=1] 
[mm*] 
0 1355 +37 | 9141 | 2,035 + 0,1% | 0,718 + 2,8%, 
5,§ 1554 + 39 | 914,1 1,370 +0,1% | 1,241 + 2,6% 
12,1 1897 + 44 914,1 1,352 + 0,1% | 1,535 + 2,3% 
0 1198 + 35 915,5 1,393 + 0,1% | 0,939 + 2,9% 
24,4 2190 + 47 915,5 0,980 + 0,1% 2,441 + 2,1% 


Nach Durchführung der oben angegebenen 


Tabelle 3 
Korrekturen und Umrechnung des Wertes 


Bleidicke |  Spurdichte bei 24,4mm Pb auf eine Spurdichte von 
[mm] | 0,730 Neutronen/mm?- h bei0 mm Pb ergeben 
0,730 + 3,0%, sich die endgültigen Spurdichten der Tab. 3. 
1,196 + 4,5%, In den Fehlern der Spurdichten der Tab. 3 
iit + 440% sind die Volumenfehler sowie die Fehler der 


Es soll jetzt ein Vergleich der Meßwerte der Tab. 3 mit einer exponentiellen 
Abnahme der erzeugenden Strahlung in Blei vorgenommen werden. Die 
Gleichung für die integrale Erzeugungskurve lautet: 


J,=A(l—e*) + B. (3) 
hire Für die Berechnung von (3) wurde für A die von George und Jason”) 
bei großen Bleidicken bestimmte Absorptionslinge von (310 + 20) g/em? 
eingesetzt. Dieser Wert stimmt innerhalb der Fehlergrenzen mit den Re- 
sultaten anderer Autoren überein (siehe z. B.17)). Der Wert von A wurde aus 
den experimentellen Daten nach der Methode der kleinsten Quadrate zu 
15,425 Neutronen/mm® h bestimmt. Die Konstante B entspricht dem Wert 
0,730 Neutronen/mm? h für x = 0 (Tab. 3). io 
16) E.P. George u. A. C. Jason, Proc. physic. Soc. 62, 243 (1949). alk at 
1) K.H. Lindenberger u. P. Meyer, Z. Physik 189, 372 (1954). abet ed 
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Die Werte der theoretischen Verteilung lauten: 


Jpurdignte Tabelle 4 
i | Pb-Dick | 
-Dicke 
| [mm] "theor. mm? h | "neor.| Yo] 
15 a | 
0,70 | 100 
0 5,9 1,060 
“123,1 1,398 115 
10 20 cm Pb 


; Diese Verteilung und die MeBwerte 
Abb. 6. Integrale Neutronenintensität der Tab. 3 sind in Abb. 6 eingetragen. 
hinter Blei. Die ausgezogene Kurve a i 
entspricht der theoretischen Verteilung Wenn mit der benutzten Absorber- 
nach Gl. (3) anordnung ein Maximum zwischen dem 

I. und 2. em Pb registriert worden 

wäre, müßten die Werte bei 1.2 und 2,4cm Pb oberhalb der theoretischen 
Kurve liegen. 
Die Fehler der in Abb. 6 eingezeichneten Werte setzen sich nur aus den 
einfachen statistischen Fehlern Vn der gemessenen Spurzahlen und den 
Fehlern der Korrekturgrößen nach dem Fehlerfortpflanzungsgesetz zusammen 
(iehe Tab. 3). Rechnet man nicht mit den einfachen statistischen Fehlern — 
eine Forderung, die für eine sichere Aussage über die Existenz eines Effektes 
notwendig ist — sondern legt 96%, Sicherheit entsprechend 2 VN fest, so 
stimmen alle experimentellen Werte innerhalb dieser Fehlergrenzen mit der 
theoretischen Kurve überein. 


Dieses Ergebnis befindet sich in Übereinstimmung mit den Messungen von 
Reich") und Hogrebe®) mit Zählrohranordnungen, die ebenfalls kein 
Maximum der Übergangskurve in Blei feststellen konnten. Treiman und 
Fonger®) haben. gleichfalls Übergangsmessungen der neutronenauslösenden 
Strahlung in Pb mit Zählrohren durchgeführt. Ein statistischer Vergleich 
ihrer Meßwerte mit einer exponentiellen Abnahme der Neutronenintensität 
ergibt Übereinstimmung, wenn man mit einem Fehler von 1,3) N rechnet. 

Damit kann zusammenfassend gesagt werden, daß das Maximum der 
Übergangskurve der neutronenauslösenden Komponente hinter kleinen Blei- 
dieken nicht als statistisch gesichert angesehen werden kann. 

Dieses Ergebnis veranlaßte uns, den Übergangseffekt der Sterne hinter 
kleinen‘ Bleidicken unter Beachtung sämtlicher möglicher Fehlerquellen in 
einer weiteren Arbeit!) zu messen. Die Ergebnisse zeigen keinen Hinweis auf 
die Existenz eines Maximums zwischen 1—2 cm Pb. Dieses Resultat befindet 
sich in Übereinstimmung mit den theoretischen Erwartungen. 


Der Slowakischen, der Tschechoslowakischen und der Deutschen Akademie 
der Wissenschaften, dem Hydrometeorologischen Dienst der CSR und dem 
Büro für Technisch-Wissenschaftliche Zusammenarbeit danke ich dafür, daß 
sie die Durchführung der Arbeit ermöglichten und unterstützten. 


_ '®) I. Hauser, P. Landrock, K. Lanius, L. Mitrani u. A. Peeva, Ann. Physik 
(1) 5, 335 (1960). 
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=) Herrn A. Mikes von der Hydrometeorologischen Station auf dem Lom. 
- Fr Stit schulde ich Dank für die Hilfe, die er mir während der Durch. 
führung der Messungen leistete. 
Herrn Dr. L. Lanius danke ich für viele anregende Diskussionen und 
Hinweise. 
Vor allem möchte ich Herrn Engelhardt, Fräulein Kussatz, Fraulein 


Ww agner und Herrn Landrock für ihre große Hilfe bei der Auswertung der Zum 
Platten danken. 
Von I 
mst Pe Zeuthenbei Berlin, Kernphy sikalisches Institut der Deutschen Akademie 
RE der Wissenschaften zu Berlin. 
Bei der Redaktion eingegangen am 11. Juli 1959. 
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Zum Übergangseffekt der sternauslösenden Komponente 
der kosmischen Strahlung in Blei 


Von I. Hauser, P.Landrock, K.Lanius, L. Mitrani und A. Peeva 


Inhaltsübersicht 


Der Übergangseffekt der sternerzeugenden Komponente der kosmischen 
Strahlung wurde hinter Blei in 3000 m Höhe untersucht. Ein Maximum der 
Sternhäufigkeit hinter 1—2 cm Pb wurde nicht beobachtet. Dieses Ergebnis 
befindet sich in Übereinstimmung mit den theoretischen Erwartungen. stars) 


Einleitung 


Mehrere Autoren!~’) untersuchten in Höhen von 2600-4600 m den Häu- 
figkeitsverlauf der in Kernemulsionen registrierten Kernzertrümmerungen 
(Sterne) hinter Bleiabsorbern. Sie beobachteten ein Maximum der Stern- 
häufigkeit bei einer Absorberdicke von 1—2cm Blei. Dieses Maximum soll 
vorwiegend durch neutral ausgelöste Sterne kleinerer Armzahl hervorgerufen 
sein. 

Theoretische Abschätzungen®)?) ergaben, daß keine der bekannten neu- 
tralen sternauslösenden Komponenten der kosmischen Strahlung in Gebirgs- 
höhe ein derartiges Maximum verursachen kann. 

Die Absorption der sternauslösenden Komponente der kosmischen afi 


lung erfolgt gemäß der Beziehung: 
T= 


Die Absorptionslinge in Blei wurde von George und Jason!) zu 4 = 
310 + 20 g/cm? bestimmt, ein Wert, der innerhalb der Fehlergrenzen auch 


Er G. Berner. G.Cortini u. A. Manfredini, Physic. Rev. 76, 1792 (1949). 
2) L. TomaSkova, Csechosl. Journ. Phys. 2, 102 (1953). 

3) E. Schopper, K. H. Höcker u. G. Kuhn, Physic. Rev. 82, 444 (1951). 

4) E.Schopper, K.H.Höcker u. E. Rössle, Z. Naturforschg. 6a, 603 (1951). 
5) L. Malaspina, M. Merlin, O. Pierucci u. A. Rostagni, Nuovo Cimento 7, 145 


6) G. Belliboni, L. Fabbrichesi, L.de Marco u. M. Merlin, Nuovo Cimento 8, 
374 (1951). 

7) T. W. Warsimaschwili u. N. I. Kostanaschwili, Zurn. eksper. teor. Fiz. 33, 
1530 (1957). 

8) E. Réssle u. E. Schopper, Z. Naturforschg. 9a, 836 (1954). 

®) K. H. Hécker, H. Kuhn u. M. Ritzi, Z. Naturforschg. 10a, 386 (1955). 

1%) E.P. George u. A. C. Jason, Proc. physic. Soc. 62, 243 (1949). 
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= den Resultaten anderer Autoren übereinstimmt (s. z. B. Lindenhouill 
und Meyer!?). 

Wir führten zunächst einen statistischen Vergleich zwischen der theore. 
tisch zu erwartenden Verteilung nach (1) und den experimentellen Ergebnissen 
der verschiedenen Autoren!~’) durch. Dazu benötigten wir die tatsächlich 
bei den verschiedenen Bleidicken gemessenen Sternhäufigkeiten. Diese Werte 
wurden nur von Tomäskovä?) angegeben. Für alle anderen Arbeiten mußten 
wir die gemessenen Häufigkeiten aus den Fehlern der Relativwerte abschätzen. 
Mit der Methode der kleinsten Quadrate ermittelten wir den /,-Wert, der die 
beste Anpassung der theoretischen Kurve nach Gl. (1) an die jeweiligen Meß. 
ergebnisse ergab. Sodann berechneten wir die Wahrscheinlichkeit P, daß die 
beobachtete Verteilung als zufällige Schwankung der theoretischen Verteilung 
auftreten kann. Mit Ausnahme der Arbeit von Schopper, Höcker und 
Rössle®) (für den statistischen Vergleich ihrer Daten wurden nur die zwischen 
0 und 11cm Pb gelegenen Meßwerte benutzt) ergab sich in allen anderen 
Experimenten P>1%. Der zugehörige J,-Wert der besten Anpassung 
stimmte in allen Fällen innerhalb des dreifachen statistischen Fehlers mit dem 
jeweiligen Beobachtungswert überein. Die Übereinstimmung ließe sich in 
allen Fällen verbessern durch Berücksichtigung der Erzeugung der sekundären 
Nukleonen in Blei, die ihrerseits zur Sternbildung in der Emulsion beitragen. 

Auf Grund des Ergebnisses einer Untersuchung des Übergangseffektes 
der neutronenauslösenden Komponente durch Hauser?) und der vorstehend 
skizzierten Überlegungen wiederholten wir die Messungen von Schopper und 
Mitarbeitern®)*) über den Häufigkeitsverlauf der Sterne hinter Blei. Diese 
Autoren zählten in ihren Messungen alle Sterne mit > 2 Armen, während alle 
anderen Autoren nur Sterne mit > 3 Armen berücksichtigten. Die Bestimmung 
der Zahl der Sterne mit 2 Armen kann stark subjektiv beeinflußt sein, da 
eine eindeutige Trennung eines Sternes von einer einfachen Streuung in der 
Mehrzahl aller Ereignisse nicht möglich ist. ei 

Experimentelle Anordnung und Ergebnisse A 

Der verwendete Bleiabsorber are dieselben Dimensionen wie der von 
Se hopper und Mitarbeitern®)*) (s. Abb. 1). Langs der vertikalen Achse be- 
fanden sich in einer (1mm Wandstärke) 4 Ilford-Gö- 
ömulsionen (5 x 9cm?, 600, Schichtdicke) parallel nebeneinander. Die 

stn Emulsionen wurden vor der Exponierung 
etwa drei Tage bei konstanter Luftfeuchtig- 
du | keit aufbewahrt und dann luftdicht in der 


Messingkassette verschlossen. Durch eine 
Gammastrahl-Markierung konnte nachträglich 
die Lage der oberen Kante des Absorbers auf 


use ea Stalin (2960 m) in Bulgarien exponiert. Auch 
Abb. 1. Form des Bleiabsorbers hierin bestand also Übereinstimmung mit den 


11) K.H. Lindenberger u. P. Meyer, Z. Physik 139, 372 nes 
12 97 0 é 
) I. Hauser, Ann. Physik (7) 5 327 (1960). 


> 
S 
; den Emulsionen mit einer Genauigkeit von 
ZA YY 0,1—0,2 mm fixiert werden. Die Photoplatten 
a wurden im Absorber 28,4 Tage auf dem Vrah 
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Experimenten von Schopper und Mitarbeitern®)®), die ihre Messungen auf 
der Zugspitze (2950 m) durchführten. Etwa 75 cm oberhalb des Absorbers 
befand sich ein Holzdach (1,3 g/em?). 

Abgesehen von der 28tägigen Exponierung im Absorber befanden sich die 
Kernemulsionsschichten zwischen dem Tage der Herstellung und dem der 
Entwicklung etwa 18 Tage ohne Absorber in Meereshöhe. Da die während 
dieser Zeit in der Emulsion aufgezeichneten Sterne auf ein eventuelles Ma- 
simum der Sternhäufigkeit ohne Einfluß sind und ihre Anzahl sich auf < 10% 
abschätzen läßt. wurden sie bei der Auswertung der Platten nicht berück- 
sichtigt. 

Für jede der vier Kernemulsionen wurden die mittlere Dicke vor der Ent- 
wicklung und der Schrumpfungsfaktor bestimmt. Die Entwicklung erfolgte 
nach dem 'Temperaturentwicklungsverfahren. 

Die Platten wurden unter dem Mikroskop bei 300facher Vergrößerung 
abgesucht. Zwei Platten wurden in Sofia und zwei in Berlin bearbeitet. Wir 
registrierten alle Sterne mit > 3 Armen. Spuren von Rückstoßkernen mit 
Reichweiten < 6 4 wurden nicht als Arme gezählt. Sterne. bei denen eine 
telativistische Spur vom Ursprung in die obere Halbebene führte, betrachteten 
wir als durch ein geladenes Teil- 


chen ausgelöst. Dabei definierten l 

wir alle Spuren als relativistisch, a | 

die eine Korndichte 9< 1,5: Gain I S Absorberkante 
(min = 22 Körner pro 100 u) be- I | 

saßen. In der Gruppe der drei- 

armigen Sterne vorhandene Proton- § Z 

Proton-Streuungen wurden elimi- r a 

niert. Um radioaktive Zerfalls- 
sterne auszuschließen, wurden vo. 
Sterne, bei denen alle Arme Reich- 14 7 ur. 
weiten <50 u besaßen, nicht regi- 
striert. Jedes registrierte Ereignis 

wurde von mindestens 2 Mit- Sen 


arbeitern identifiziert. 

Abb. 2. Lage der abgesuchten Bereiche 

suchens wurden die in den Bereichen 

I-VI jeder Platte liegenden Sterne zusammengefaßt (s. Abb. 2). In Tab. 1 

ind die in jedem Bereich beobachteten Sternhäufigkeiten angegeben. (Die 

Bereiche V und VI wurde n nur in zwei der vier Platten abgesucht.) 


Tabelle 1 


Bereich Ill IV V | 
Absorbe rdicke in em Pb 0 C—1 1—2 9-8 3—4 
Zahl der Sterne mit 
- 3 Armen 326 329 333 302 290 


Die Werte der Tab. 1 müssen auf gleiche Emulsionsvolumen bezogen 
ierden, da beispielsweise die Emulsionsdicke in einer einzelnen Platte um 
twa 10%, in den Bereichen variiert. Diese Variation der Emulsionsdicke von 


Ann. Physik. 7. Folge, Bd. 5 
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Platte zu | Platte bzw. innerhalb jeder einzelnen Schicht muß bei der Bestimmung 
der Zahl der Ereignisse pro Volumeneinheit berücksichtigt werden, denn die 
Dickenschwankungen um den von der Herstellerfirma angegebenen Nennwert 
können bis zu einem Fehler von ~10% der Häufigkeit pro Volumeneinheit 
führen. Die Bestimmung des Emulsionsvolumen jedes Bereichs wurde mit 
einem Fehler von 2,5%, durchgeführt. 

In Tab. 2 sind die in den Bereichen I—VI ermittelten Sternhäufigkeiten 
pro cm? und Tag mit ihren einfachen statistischen Fehlern angegeben. 


Tabelle 2 
Absorberdicke in cm Pb 0 | 0—1 1—2 | 2—3 | 3—4 | 4—5 


Zahl der Sterne mit 
> 3 Armen pro cm?, d | 11,4+0,6 


11,2+0,6 | 11,6-49,6 | 10,5-+0,6| 10,140,8 , 9,0-+0,8 


Diskussion der Ergebnisse 


In Abb. 3 sind die Werte der Tab. 2 eingezeichnet. Die ausgezogene Kurve 
wurde nach Gl. (1) mit A = 310 g/em? berechnet. Vergleicht man mit Hilfe 
des y?-Testes die experimentelle mit der 
theoretischen Verteilung, so ergibt sich, 
daß die gemessenen Werte i in 66% aller Fälle 
als zufallige Abweichungen der hen 
4 - N 
Kurve auftreten können. Führt man den 
Vergleich nur für die neutral ausgelösten 
J Sterne mit 3—5 Armen durch, so wird die 
_ Übereinstimmung zwischen den MeB- und 
den Erwartungswerten noch besser. 


Loni der Sterne mit 23 Armen 
pro cm’und Tag 


0 1 2,9 4 SAbsorberdiche 


inem Pb Zusammenfassend können wir fest- 
stellen, daß unsere Messungen keinen Hin- 


hinter Blei ur weis auf die Existenz eines Maximums der 


Sternhäufigkeit hinter 1—2cm Blei ent- 
halten. Dieses Ergebnis befindet sich in Übereinstimmung mit den theo- 
retischen Erwartungen. 


Wir danken der Bulgarischen und der Deutschen Akademie der Wissen- 
schaften für die Ermöglichung der Gemeinschaftsarbeit. 

Herrn H. W. Meier und Herrn A. Strigatscheff danken wir für die 
Hilfe, die sie bei dem experimentellen Aufbau des Absorbers leisteten. 


Zeuthen bei Berlin, Kernphysikalisches Institut der Deutschen Akademie 
der Wissenschaften zu Berlin und 
Sofia, Physikalisches Institut der Bulgarischen Akademie der Wissen- 
schaften. 
Bei der Redaktion eingegangen am 11. Juli 1959. 


Die 
tration 
verallg 
gewinr 
tionsal 
zierend 
Schliiss 
mittel 


Lut 


als Sz 
Lumin 
UV-Lic 
bierte § 
iibertre 
in den 
Autore 
Swan! 
deutig 


Das 
durch 


ı) H 

2) M 

3) M 

4) H 
quiums | 
Braunse 
5) J 
J 
7) R 
8) A 
%) H 


. 
Die 1 
| q 
= 
44 
| 
E 
= 


en- 


die 


mie 


. . . . 
Lösungen mit einer gelösten lumineszierenden Substanz _ 


7 
Von H.K. Bothe 
> Mit 12 Abbildungen 


Inhaltsiibersicht 


Die von Kallmann u. Furst?) abgeleitete Funktion für die Konzen- 
trationsabhängigkeit der Lumineszenzintensität flüssiger Lösungen wird so 
verallgemeinert, daß sie auch für die festen Polystyrol-Lösungen Gültigkeit 
gewinnt. Aus dem charakteristischen Unterschied zwischen der Konzentra- 
tionsabhängigkeit der Lumineszenzintensität fester und flüssiger lumines- 
zierender Lösungen bei Anregung mit energiereicher Strahlung werden 
Schlüsse auf die Art des Mechanismus der Energieübertragung vom Lösungs- 
mittel auf die gelösten lumineszierenden Moleküle gezogen. 


Einführung 


Lumineszierende organische Lösungen finden in immer stärkerem Maße 
als Szintillatoren für Szintillationszähler Verwendung. Ihr relativ gutes 
Lumineszenzvermögen bei Anregung mit Kernstrahlung oder kurzwelligem 
UV-Licht beruht in erster Linie darauf, daß die vom Lösungsmittel absor- 
bierte Strahlungsenergie von diesem auf die gelösten lumineszierenden Moleküle 
übertragen wird (Kallmann und Furst!)). Diese Energieübertragung war 
in den letzten Jahren Gegenstand eingehender Untersuchungen verschiedener 
Autoren (Kallmann, Furst!)?)?) und Brown‘), Birks®) und Cameron®), 
Swank und Buck’), Schmillen®), Knau%) und andere), ohne daß ein- 
deutig geklärt werden konnte, welcher Mechanismus ihr zugrunde liegt. 


Das Ziel der vorliegenden Arbeit war, Kriterien ausfindig zu machen, 
durch welche eine Entscheidung zwischen den verschiedenen bisher vorge- 


1) H.Kallmannu.M. Furst, Physic. Rev. 79, Nr. 5, 857 (1950). 

2) M. Furstu. H. Kallmann, Physic. Rev. 85, Nr. 5, 816 (1952). 

3) M. Furstu. H. Kallmann, Physic. Rev. 94, 503 (1954). 

4) H. Kallmann, M. Furst u. F. H. Brown, Vortrag anläßlich des internat. Kollo- 
quiums „Halbleiter und Phosphore““ 1956 in Garmisch-P., ersch. Verlag Friedr. Vieweg, 
Braunschweig 1958, S. 269. 

5) J. B. Birks, Physic. Rev. 94, 1567 (1954). Fund 

6) J. B. Birksu. A. J. W. Cameron, Proc. physi. Soc. 72, Nr. 1, 53 (1958). 

7) R.K. Swank u. W. L. Buck, Physic. Rev. 91, 927 (1953). 

°) A. Schmillen, Z. Physik 150, Nr. 2, 123 (1958), 

®) H. Knau, Z. Naturforsch. 12a, 881 (1957). 
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schlagenen und im folgenden zunächst kurz beschriebenen Energieübertra- 
gungsmechanismen herbeigeführt werden kann. 


Da eine Energieübertragung nur bei solchen Lösungsmitteln erfolgt, die 
in reinem Zustand ein Lumineszenzvermögen besitzen, liegt es nahe anzuneh- 
men, daß die Übertragung der Anregungsenergie vom Lösungsmittel auf die 
gelösten Moleküle durch die Reabsorption der Lösungsmittellumineszenz erfolgt. 
In den meisten Fällen ist jedoch die Quantenausbeute der von den reinen 
Lösungsmitteln emittierten Lumineszenz viel zu gering, um die Energie- 
übertragung allein auf die Reabsorption dieser Lumineszenz zurückzuführen. 
Kallmann, Furst und Brown?) nehmen daher an, daß die vom Lösungs- 
mittel absorbierte Anregungsenergie durch intermolekulare Stöße (collisions) 
zunächst von Lösungsmittelmolekül zu Lösungsmittelmolekül und schließlich 
auf die gelösten lumineszierenden Moleküle übertragen wird. 


Birks stellt dagegen die Hypothese auf, daß auch die höheren Anregungs- 
zustände der organischen Moleküle unter Lumineszenzemission in den Grund. 
zustand zurückkehren, diese Lumineszenz aber infolge der starken Reab- 
sorption nicht aus dem Szintillator austreten kann. Unter dieser Voraus- 
setzung ist es ihm möglich, die Energieübertragung in lumineszierenden organi- 
schen Lösungen ausschließlich auf eine Emission und Reabsorption von 
Photonen zurückzuführen. 


Da selbst bei Szintillatoren von nur wenigen u Dicke die von Birks po- 
stulierte Lumineszenz der höheren Anregungszustände nicht beobachtet 
werden konnte, muß die Reichweite der diesen Anregungszuständen ent- 
sprechenden Photonen kleiner als ein w und somit von der Größenordnung 
ihrer Wellenlänge sein. Es ergibt sich die Frage, ob man in diesem Falle noch 
von einer Emission und Reabsorption von Photonen sprechen kann, oder ob 
bereits eine unmittelbare quantenmechanische Resonanz (Förster!")) zwischen 
den emittierenden und absorbierenden Molekülen vorliegt. 


Schmillen®) zeigte, daß sich die Energieübertragung in Lösungen auch als 
Diffusion angeregter Zustände (Excitonen) auffassen läßt. Er läßt aber 
offen, in welcher Weise sich diese Excitonen in organischen Substanzen 
bilden können. Wie Heller und Markus!!) am Beispiel der anorganischen 
Kristallphosphore nachwiesen, läßt sich das Modell der Excitonenwanderung 
auch auf eine Energieübertragung durch quantenmechanische Resonanz 
zurückführen. Exeitonendiffusion und quantenmechanische Resonanz sind 
daher als zwei verschiedene mathematische Modelle des gleichen Vorganges 
‚aufzufassen, von denen allerdings das der quantenmechanischen Resonanz 
wesentlich konkreter ist. Die vorliegende Arbeit wird sich daher vor allem 
mit den Energieübertragungsmechanismen 


A. Intermolekularer Stoß (collision) (Kallmann, Furst und Brown’). 
B. Emission und Reabsorption von Photonen (Birks®)), 
C. Quantenmechanische Resonanz (Förster!®)) 


befassen. 


1) Th. Förster, Ann. Physik 6, Nr. 2 55 (1948). . ‘| 
W.R. Heller u. A. Marcus, Physic. Rev. S84, 809 (1951). 
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Die Konzentrationsabhängigkeit der Lumineszenzintensität von Lösungen 


Kallmann und Furst und auch Birks führen die Konzentrationsab- 
hängigkeit der Lumineszenzintensität von Lösungen auf die Häufigkeit be- 
stimmter an der Energieübertragung beteiligter Prozesse zurück (Tab. 1 
und 2). Auf Grund der jeweils berücksichtigten Prozesse ergeben sich für 


Tabelle 1 
Prozesse, von denen nach der Theorie von Kallmann und Furst?) das Lumincszenz- 
vermögen einer Lösung der lumineszierenden Substanz A im Lésungsmittel Z abhängt 


Wahrschein- 


Bezeichnung des Prozesses Prozeß 
lichkeit 


Absorption der anregenden Strahlung E durch das 


Lösungsmittel . . ‘ E+L—>L 
Innere Löschung im Lösungsmittel. te L>L 
Ü Zr der Anregungsenergie auf die Mole - 

Innere Lése hung der Mole küle A+—>A 
Konzentrationslöschung . . . ........ AP +A, 
Eumineszenzemission . . . 4+>4A-+hr, 


Tabelle 2 
Bei der Energieübertragung nach Birks?) berücksichtigte Prozesse in eirer Lésurg der 
Substanz 4 im Lösurgsmittel Z 


Wahrscheinlich- 


Bezeichnung des Prozesses Prozeß keitelakton 
Absorption d.anregend. Strahlung .. .... E+L-L: 1 
Innere Löschung im Lösungsmittel. . . . . . L+>L | K,, 
Lumineszenz der Lösungsmittelmoleküle L+—>L+hr, Ky, 
Selbstabsorption der Lésungsmittellumireszenz L>Ll: Ky, 
Absorption der Lösungsmittellumineszenz durch 4 hv, + 4 At Kia 
Innere Löschung der Moleküle A . ..... 4r—A | Ku 
Lumineszenz der Moleküle A . . ...... | Kin, 
Selbstabsorption der Moleküle A . . . hv, +44 | Ku 
Emission der Lumineszenzvon d . ...... Emission 


hv — Lumineszenzquanten. 


Molekulare Lumineszenzquantenaushcute des Lésurgsmitt« 
=* 
K 
Yor 


Ky, + Ky, 
molekulare Lumineszenzquantenausbeute g., der gelösten Sulstarz: 


K ps 


die Übertragung der Anregungsenergie innerhalb einer Lösung der Substanz A 
im Lösungsmittel Z die in Abb. 1 und 2 dargestellten Schemata. Da sich 
„quantenmechanische Resonanz‘ von einer Emission und Reabsorption von 
Photonen allein bezüglich der Entfernung unterscheidet, über welche die Ener- 
gieübertragung erfolgt, hat Energieübertragungsschema Abb. 2 auch bei einer 
Energieübertragung durch quantenmechanische Resonanz Gültigkeit. Nach 
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beiden v bertragungsschemata ergibt sich für die Abhängigkeit der Lumines. Lösı 
zenzintensität J, von der Konzentration [A] einer im Lösungsmittel Z ge-| Polyvin 


lösten lumineszierenden praktis: 


Substanz A: Styrol- 

J P [A] und sin 
AT [AD (R+[Ap’ | Frage, 
Lösung 


(1) sich di 
wobei die Größe der} zwische 
Parameter P, Q und R| unmitte 
von der Wahrscheinlich- | übertra 
keit der in Tab. 1 und 2 | des Lö: 
aufgeführten Prozesse | quanteı 


abhängt (Tab. 3). Ein 
Kallmann und | sofern 

_ Furst’) konnten nach- | quante 
weisen, daß sich | Orienti 
die Konzentrationsab- Du: 


hängigkeit der Lumi- | wird, : 
neszenzintensität flüssi- | Funk ti 
ger Lösungen sehr gut | welche 
durch die Funktion (1) | zwisch 
darstellen läßt. Für Lö- Da: 


Lumineszenzemission L umineszenzemission I 3 
Abb. 1. Energieiibertra- Abb. 2.  Energieiibertra- x 
gungsschema nach Kall- gungsschema nach Birks®) mittel ein merkliches | Weise 
mann und Furst?) Lumineszenzvermögen Bei 


besitzt, kann (1) jedoch | kann ı 
keine Gültigkeit besitzen, da dann die Lumineszenzintensität J, der Lösung | wird, : 
für A — 0 nicht wie Funktion (1) gleich Null, sondern gleich der Lumines- | beobac 
zenzintensität des Lösungsmittels wird. Der Anteil der Lösungsmittellumines- | stark i 
zenz an der Energieübertragung kann jedoch, wie noch gezeigt wird, durch | die gel 
sinnvolle Erweiterung der Funktion (1) berücksichtigt werden. 


Tabelle 3 

Die Parameter P, Q und R der Funktion (1) als Funktion der Wahrscheinlichkeitskoeffi- | Für de 
zienten der in Tab. 1 und 2 aufgeführten Prozesse. (Durch 6 wird die Lumineszenz- | innereı 
quantenausbeute der Lösungen mit der in willkürlichen Einheiten (J) gemessenen Lu- (Tab. : 
mineszenzintensität J, in Beziehung gesetzt) 5 


Rsind 
Funktion nach Ubertragungs- | Funktion nach Übertragungs- | R ist. 
Parameter schema Abb.1 (Theorievon | schema Abb. 2(Theorie von 
Kallmann und Furst?)) | Birks°)) 
P ne Odor Goa Kea Nach « 
ß Kia (1 904) 
| sungen 
| bei ge 
(1 — Gx) | 
ergib 
t,°% gıbt 
K 


re 
ar. . ul 
Primaranregung 
7) 
LA 
hg - 
K A 
a (A) 
K 
FA \ Ky, ) 
(A) 
ti B A + hv, 
1 
4 / te Kea 


i H.K. Bothe: Die Energieübertragung in lumineszierenden Lösungen 343 


Lösungen lumineszierender organischer Substanzen in Polystyrol oder 
Polyvinyltoluol haben als Plastmasseszintillatoren für Szintillationszähler 
praktische Bedeutung gewonnen. Sie werden durch die Polymerisation von 
Styrol- bzw. Vinyltoluol-Lösungen der betreffenden Substanzen hergestellt 
und sind wie flüssige Lösungen völlig homogen und isotrop. Es ergibt sich die 
Frage, in welcher Weise sich die Energieübertragung in den polymerisierten 
Lösungen von der in flüssigen Lösungen unterscheidet. Am stärksten müßte 
sich die Polymerisation der Lösungen auf eine Energieübertragung durch 
zwischenmolekulare Stöße auswirken, da dieser Übertragungsmechanismus 
unmittelbar mit den intermolekularen Kräften verknüpft ist. Die Energie- 
übertragung durch Photonen dürfte dagegen kaum vom Aggregatzustand 
des Lösungsmittels abhängen, sofern sich nicht gleichzeitig die Lumineszenz- 
quantenausbeute des Lösungsmittels ändert. 

Eine Energieübertragung durch quantenmechanische Resonanz kann in- 
sofern vom Aggregatzustand des Lösungsmittels beeinflußt werden, als die 
quantenmechanische Resonanz zwischen zwei Molekülen auch von deren 
Orientierung zueinander abhängt. 

Durch die experimentellen Arbeiten, über welche im folgenden berichtet 
wird, sollte zunächst geklärt werden, unter welchen Voraussetzungen die 
Funk tion (1) auch auf die festen Polystyrol-Lösungen anzuwenden ist, und 
welcher charakteristische Unterschied bezüglich der Parameter P, Q und R 
zwischen flüssigen und polymerisierten Lösungen besteht. 

Das Lumineszenzvermögen von reinem Polystyrol ist wesentlich besser 
als das der meisten flüssigen Lösungsmittel, es läßt sich jedoch auf folgende 
Weise berücksichtigen: 

Bei der vom Lösungsmittel an die gelösten Moleküle abgegebenen Energie 
kann man unterscheiden zwischen einem Teil S, welcher nicht reabsorbiert 
wird, so daß er bei Abwesenheit der gelösten Moleküle als Lumineszenzlicht 
beobachtet werden kann, und einem Teil 1—S, für den die Reabsorption so 
stark ist, daß er nur in Erscheinung tritt, wenn eine Energieübertragung auf 


die gelösten Moleküle erfolgt. Die Gl. (1) bekommt dann die Form 
SPI), 
+ [AD + TAD + [AD + [AD 


Für den Teil S kann definitionsgemäß der Wahrscheinlichkeitskoeffizient der 
inneren Löschung 1/t, (Tab. 1) bzw. der Selbstabsorptionskoeffizie nt Ä,ı 
(Tab. 2) gleich Null gesetzt werden. Dann wird nach Tab. 3 Q’ = P. P und 
R sind für beide Anteile gleich, so daß Q’’= Q, P’= P” = Pund R’ = R” = 
Rist. (2) geht daher über in: 

SP (1-8) P[A] 


4=R+ [4] @+ 


Nach dieser Funktion strebt J, nunmehr in Ubereinstimmung mit den Mes- 
sungen von Kallmann und Furst’), Swank und Buck’), Krenz!) 
bei geringen Konzentrationen einem endlichen Grenzwert zu, für den sich 


| 12) F,H.Krenz, Trans. Faraday Soc. 51, 172 (1955). „ehe 
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Experimentelle Arbeiten 


Die experimentellen Arbeiten hatten zum Ziel, die Gültigkeit der Funk- 
tion (3) für die Konzentrationsabhängigkeit der Lumineszenzintensität von 
Polystyrollösungen nachzuweisen bzw. die Parameter P, Q, R und S dieser 
Funktion für eine Reihe 
lumineszierender Substanzen 
in Polystyrol zu ermitteln. 
Hierzu wurden die in Tab. 4 
aufgeführten Substanzen in 
Styrol gelöst und diese 
Lösungen ohne Zusatz eines 
Katalysators unter Vakuum 
bei 140° C polymerisiert. Die 
so gewonnenen festen Lö- 
sungen wurden auf Zylinder 
von gleicher Höhe (25 mm) 
und gleichem Durchmesser 
Abb.3. Plastmasseszintillatoren (etwa natürl. Größe) (13 mm) abgedreht und po- 

liert (Abb. 3). 

Zur Ermittlung der Konzentrationsabhängigkeit der Lumineszenzintensität 
genügten relative Messungen, wobei allerdings an die Genauigkeit der einzelnen 
Meßwerte hohe Anforderungen zu stellen waren. 

Es mußten vor allem alle Faktoren ausgeschaltet werden, welche bei der 
theoretischen Ableitung der Funktionen (1) und (3) nicht berücksichtigt 
worden waren. Um die Absorption des Lumineszenzlichtes im Szintillator mög- 
lichst gering zu halten, durfte die anregende Strahlung nur wenig in den Szin- 
tillator eindringen. Die Lumineszenzanregung erfolgte durch die 3-Strahlung 
eines TP?4-Priparates (max. Reichweite der Strahlung in Polystyrol etwa 

3 mm, Halbwertsschichtdicke 
0.4mm), und in einigen Fällen 
auch durch weiche Röntgenstrah- 
lung (CuK,-Linie = 1,5A. Halb- 

Photokathode  Wertschichtdicke 1,3 mm). Die an- 
regende Strahlung traf auf eine 
der beiden ebenen Stirnflächen 


Richtung der anregenden 


Chrahlınn 
g 


“ der zylindrischen Szintillatoren. 
Szintillator Gemessen wurde die Intensität 
Abb. 4. Geometrie der Meßanordnurg des aus dieser Fläche emittierten 


Lumineszenzlichtes (Geometrie der 
MeBancrdnung Abb. 4). Die Messung der relativen Lumineszenzintensität er- 
folgte mit einem Photosekundärelektronenvervielfacher RCA 5819, dessen 
Spannung elektronisch stabilisiert war. Die der Lumineszenzintensität 
proportionale Stromstärke am Ausgang des PSEV wurde an einem Skalen- 
galvanometer abgelesen. 


Meßergebnisse 

In den Diagrammen Abb. 5—12 sind die gemessenen Werte der Lumines- 
zenzintensität von Polystyrollösungen in Abhängigkeit von der Konzentration 


dargestellt. Da das Lumineszenzvermögen eine molekulare Eigenschaft ist, 
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war es ee: die Konzentration in Mol/Liter anzugeben. Wegen der 
Einfachheit bei der Herstellung der Lösungen wurden die Konzentrationen 


& 


Tabelle 4 N 
Verbindungen, deren Lumineszenzvermögen in Polystyrol untersucht wurde ry Fu 
Farbe der 
Verbindung Strukturformel Lumines- 
zenz 
| 
| 
Naphthalin L | | UV 
Anthrazen | | | blau 


a RY 


+ 


Stilben Din. UV 


p-Terphenyl «K > UV 


11.4.-Tetraphenylbutadien (TPB) blau 


Phenyl-biphenylyl-oxadiazol(PBD) 


Ann. Physik. 7. Folge.Bd.5 

way 
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mit ß- und Röntgenstrahlung. Erläute- 
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auf monomeres Styrol bezogen, 
Bei der Polymerisation findet 
eine geringe Erhöhung der 
Dichte und damit auch der Kon- 
zentration statt. Die Anregung 


100} 


oO 


Sn 
oO 


0’ 7 
Mol Anthracen in 1000cm? Styrol 


Abb. 5. Konzentrationsabhängigkeit der Lumi- 
neszenzintensität von Anthrazen in Polystyrol 
bei Anregung mit ß- und Röntgenstrahlung. Meß- 
werte: Anregung mit ß-Strahlung: O, Anregung 
mit Röntgenstrahlung: x. Die eingezeichneten 
Kurven ergeben sich durch Einsetzen der ent- 
sprechenden in Tab. 5 enthaltenen Parameter in 
Gl. (3). Angabe der relativen Lumineszenzinten- 


berechnen. Dieses ist zunächst nach QR, Q 


erfolgte in allen Fällen mit ß. 


Strahlung (T12%), bei Fluoran- 
then, Anthrazen, Stilben, PBD 
nd Tetraphenylbutadien auch 


mit weicher Röntgenstrahlung 


(CuK,-Linie). Die Einheiten der 
relativen Lumineszenzintensität 
sind bei den verschiedenen An- 
regungsarten nicht die gleichen. 

Aus den Meßwerten J,, J;, 
I, der Lumineszenzintensität 
bei vier verschiedenen Konzen- 
trationen [A,], [Aa]; [4g] und 
[A,] lassen sich die Parameter 
P,Q, Rund S auf Grund des 


. a sität in willkürlich festgelegten Einheiten Gleichungssystems 

QR+1,[4)-(Q+ R)— PQS— [A]: P=—h[AP 

“4 I, -QR + I, [Ag] -(Q + R)— PQS — [A,]- P = — 1, [4,7 (4) 

R + I; [Ag] -(Q + R)— PQS — [A,]- P = — I; [4,7 
I,-QR+ 1,[Aq]-(Q+ R)— PQS— [A,]: P=— 1, 


+ R, PQS und P aufzulösen. 


Aus QR und Q + R erhält man die quadratische Gleichung 
_Q+R 
Q bzw. + | 4 QR. 


ital 


‘ative Lumineszenzint 


Konz 
10°’ 7 
Mol Fluoranthen in 1000cm? Styrol 


Abb. 6. Konzentrationsabhangigkeit der Lumines- 
zenzintensität von Fluoranthen in Polystyrol bei 


lo”? 


rungen siehe Abb. 5 


- 


Wählt man für Q den kleineren 
der beiden möglichen Werte, 
so ergibt sich für S in allen 
Fällen ein Wert um 0,2. Dieser 
steht in guter Übereinstim- 
mung mit den Ergebnissen der 
Messungen von Swank und 
Buck’) bzw. Krenz!2), wel- 
che ergaben, daß etwa 20% 
der Lumineszenz von Poly- 
styrollösungen durch das aus 
reinem Polystyrol emittierte 
Lumineszenzlicht angeregt 
wird. Die aus den Meßwerten 
berechneten und in Tab. 5 zu- 
sammengestellten Werte für 
die Parameter P, Q, R und 8 
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Tabelle 5 


347 


Parameter P, Q, R und S der Konzentrationsabhängigkeit der Lumineszenzintensität 
von Polystyrollésungen 


| Parameter 
Gelöste lumines- | Anregende P Q R N 
zierende Substanz Strahlung [Z][Mol/Lit] | [Mol/Lit] [Mol/Lit] 
Fluoranthen 201 0,01818 0,4675 0,20 
Röntgen 201 0,01818 0,4675 0,20 
Anthrazen B 955 0,1178 2,231 0,19 
Röntgen 891 0.1178 2.231 0.19 
Stilben 62 0,01632 0,127 | 
Röntgen 57 0,01632 0,1727 0,26 
Tetraphenylbutadien ae 196 0,01237 0,2126 0,19 
Röntgen 140 0.01207 0,1618 | 0,20 
Phenyl-biphenylyl- 471 0,01286 0,3724 | 0,11 
oxadiazol Röntgen 314 0,01365 0,2582 0,11 
Naphthalin B 277 0,1237 0,6883 | 0,17 
Pyren 5 219 0,016 0,2374 | 0,17 
p-Terphenyl wo 346 0,02482 0,282 | 0,136 


dürfen, da sie definitionsgemäß Funktionen von Wahrscheinlichkeitskoeffi- 
zienten (Tab. 3) sind, weder negativ noch imaginär werden. Die Möglich- 


keit, durch 4 Meßwerte eine 
Kurve der Gl. (3) zu legen, ist 
hierdurch eng begrenzt, und 
schon geringe Fehler der MeBß- 
werte führen zu unzulässigen 
Werten für einen oder mehrere 
der berechneten Parameter. Die 
ausgezogenen Kurven in den Dia- 
gammen Abb. 5—12 wurden 
durch das Einsetzen der Para- 
meter in Gl. (3) ermittelt. 

Die Konzentrationsabhängig- 
keit der Lumineszenzintensität 
im Konzentrationsbereich von 
etwa 5 - 10% — 2- 10-1 Mol/Liter 
wird sehr gut durch die Gl. (3) 
wiedergegeben. Bei höheren Kon- 
zentrationen liegen die gemes- 
senen Werte teilweise über den 


Relative Lumineszenzintensität 


Sn 
S 


oO 


0° 10°’ 


Mol Stilben in 1000cm? Styrol 


Abb. 7. Konzentrationsabhängigkeit der Lumi- 
neszenzintensität von Stilben in Polystyrol bei 
Anregung mit ß- und Röntgenstrahlung. Er- 


läuterungen siehe Abb. 5 


sich aus der Gl. (3) ergebenden Werten. 


Dies mag daran liegen, daß bei diesen Konzentrationen ein Teil der anregenden 
Strahlung direkt von den gelösten Molekülen absorbiert wird. 


Die Einheit (I) der relativen Lumineszenzintensität war bei der Anregung 
durch Röntgenstrahlung so gewählt worden, daß bei Fluoranthen der Wert 


enen tibereinstimmte. Die 
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Anregungsart (Röntgen- oder 8-Strahlung) hat in allen Fälle 
nur geringen Einfluß auf die Größe der Parameter Q und R. 
Diskussion der MeBergebnisse 
Interessant ist ein Vergleich der sich fiir feste Polystyrollésungen er. 
gebenden Parameter mit den von Kallmann und Furst fiir die gleichen 


Tabelle 6 
Parameter Q und R der Konzentrationsabhängigkeit der Lumineszenzintensität einiger 
lumineszierender organischer Vi erbindungen i in Polystyrol und in flüssigen Lösungsmitteln, 
(Die Parameter der flüssigen Lösungen ergaben sich aus den von K. u. F3) ermittelten 
Werten nach Umrechnung von der Dimension g/Liter auf Mol/Liter) 


Lumineszierende Substanz | Lösungsmittel [Mol [Mol/Liter] 
Anthrazen Polystyrol 0,1178 2,231 
Xylol 0,00298 0,0208 
gr Benzol 0,00236 0,0208 
Toluol 0,00314 0.0205 
ae . | Polysty rol | 0,01918 0,4675 
Benzol 0,00208 1,09 
ER Toluol 0,00223 1,145 
e rphenyl Poly styrol | 0,0248 0,282 
| Xylol .0,00261 0,1925 
| Benzol | 0,00265 0,222 
Toluol | 0,00369 0,156 
Stilben Polystyrol | 0,01632 0,1727 
Xylol 0,00211 0,35 
Substanzen aus der Konzen- ‘ 
Lösungen berechneten (Tab.6). [ 2 
Zu diesem Vergleich mußten 800+ 


/ 
\ 
/ \ 


die Werte von Kallmann, 
welche sich auf eine Konzen- 
trationsangabe von g/Liter be- 
ziehen, auf die Dimension 
Mol/Liter umgerechnet wer- 
den. Verglichen werden kön- 
nen nur die Werte der Para- 
meter Q und R, da in P die 
jeweils willkürlich festgelegte 


> 
a. 


Relative Lumineszenzintensitat 


tinheit [J] der Lumineszenz- | 
intensität eingeht. @ ist bei 
0 0 10 
den festen Polystyrol-Lésun- 
gen etwa 10mal, im Falle An- 


Abb. 8. Konzentrationsabhängigkeit der Lumines- 
S0mal zenzintensität von Phenyl-biphenylyl-oxadiazol in 
größer als bei den flüssigen. Polystyrol bei Anregung mit ß- und Röntgen- 


Dies läßt sich nach den Defi- strahlung. Erläuterungen siehe Abb. 5 
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nitionen der Tab. 3 entweder auf eine Erhöhung der Wahrscheinlichkeit 1/t, 
der inneren Löschung im Lösungsmittel oder auf eine Verringerung der Wahr- 
scheinlichkeit & des Energieüberganges vom Lösungsmittel auf die gelösten 
Moleküle zurückführen. Die 

erste Aussage steht im ER 


sache, daB gerade reines 
Polystyrol ein besseres 
Lumineszenzvermögen be- 
sitzt als die betreffenden 
flüssigen Lösungsmittel. 
Ebenfalls ist nicht einzu- 
sehen, warum die Wahr- 
scheinlichkeit einer Ener- 
gieübertragung durch Stoß 


in festen Lösungen geringer 70” 07 0” 
sein sollte als = flüssigen Mol Tetraphenylbutadien in 1000 cm?Styrol 


da doch die Wechselwir- Abb. 9. Konzentrationsabhängigkeit der Lumineszenz- 
En eid d Mole intensität von 1.1.4.4.-Tetraphenylbutadien in Poly- 
ang zwischen den Mole-  styrol bei Anregung mit ß- und Röntgenstrahlung. 
külen in festen Lösungen Erläuterungen siehe Abb. 5 
wesentlich stärker ist als 
in flüssigen, vor allem dann, wenn die gelösten Moleküle direkt in die 
Lösungsmittelmoleküle eingebaut sind, wie es von Hinrichs!) im Falle 
Anthrazen in Polystyrol nachgewiesen werden konnte. 

Bei einer Energieübertragung durch Photonen oder quantenmechanische 
Resonanz hat der Parameter 
Q nach Tab. 3 die Form: 
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500+ 


Relative Lumineszenzintensitat 


700 
600| 
500 
400 
3001 
200| 


K 2 
Q= LL Vor (5) 


Hieraus und aus der Beob- 
achtung, daß@ bei den festen 
Polystyrollösungen etwa 10- 
bis 30mal größer ist als bei 
Lösungen mit flüssigem Lö- 


Relative Lumineszenzintensitat 


sungsmittel, geht hervor, daß vol 
entweder Konz 
3 2 
Kyıl ı/Ky.) fe fest 10 10 0 7 
>10 (6 Mol p-Terphenyl in 1000cm? Styrol 
Abb.10. Konzentrationsabhängigkeit der Lumi- 
oder aber I “N neszenzintensität von p-Terphenyl in Polystyrol 
bei Anregung mit f-Strahlung. Erläuterungen 
(1 — oz) fest > 10 m’ (7) siehe Abb. 5 


(1 — qp,) fliiss. — 
ist. Bei einer Energieübertragung durch Photonen dürfte das Verhältnis der 
Selbstabsorption (K,,) zur Fremdabsorption (K,,) der Lösungsmittellumi- 
neszenz nicht wesentlich vom Aggregatzustand des Lösungsmittels abhängen, 


13) H. Hinrichs, Z. Naturforsch. 9a, 617 (1954). 
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so daB ii ‘deat Art der Energieübertragung nur die RE (7) in Frage 
käme. Diese ist nur erfüllt, wenn die molekulare Lumineszenzquantenausbeute 
Yor der flüssigen Lösungsmittel nahezu 1 und die vom Polystyrol wesentlich 
kleiner ist. Ein derartiger Unterschied zwischen den molekularen Lumineszenz. 
quantenausbeuten ist jedoch wenig wahrscheinlich. 
Sehr gut läßt sich die Bedingung (6) bei einer Energieübertragung dure 
quantenmechanische Resonanz deuten. Die Erhöhung der Wahrscheinlich 
keit K,, der Energieüber. 
tragung vom Lösungsmittel. 
molekiil zu Lésungsmittelmole. 
kül gegenüber der Wahrschein. 
lichkeit K,, der Energieüber- 
_ tragung vom Lösungsmittel auf 
die gelösten Moleküle bei Poly. 
könnte bei dieser 
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fe Art von Energieübertragungs- 

60 mechanismus auf eine starke 

Resonanz zwischen den Glie. 
20 Konz. 

dern der Polystyrolketten und 
1? 10? 0 7 


damit auf eine bevorzugte 
Mol Naphthalın in 1000cm? Styrol Energieübertragung entlang der 
Abb. 11. Konzentrationsabhängigkeit der Lumi- Polystyrolketten zurückgeführt 
gy merianangg von Naphthalin in Polystyrol werden. Ein Einpolymerisieren 
ei Anregung mit ß-Strahlung. Erläuterungen 
siehe Abb. 5 der lumineszierenden Moleküle, 
wie es z. B. nach Hinrichs®) 
bei Anthrazen in Polystyrol 
_ der Fall ist, erschwert offenbar 
die Energieübertragung auf die 
betreffenden Moleküle; denn 
bei Anthrazen in Poly- 


AD 
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styrol ist Q relativ groB. 
Die Veränderung des Para- 
#4 meters R beim Übergang von 
flüssigen zu festen Lösungen ist 
nicht so charakteristisch wie 
— die des Parameters Q. Für Ter- 
Mol Pyren in 1000cm? Styrol Stilben ist B 
esten Lösungen von der glei- 
Abb. 12. Konzentrationsebhängigkeit der Lumi- chen Größenordnung wie bei 


größer, für Fluoranthen dagegen 
5mal kleiner. Dies liegt an der unterschiedlichen Größe von R bei flüssigen 
Lösungen. Bei den festen Polystyrollösungen hat R für alle untersuchten 
Substanzen mit Ausnahme von Anthrazen die gleiche Größenordnung. Bei 
flüssigen Lösungen hängt offenbar die Wahrscheinlichkeit eines der Energie- 
übertragungsprozesse, durch welche nach Tab. 3 die Größe von R bestimmt 
ist, wesentlich stärker von der Art der gelösten lumineszierenden Substanz ab 
als bei festen Polystyrollösungen. Man kann annehmen, daß durch die unter- 
schiedliche Assoziation der Moleküle in flüssigen Lösungen die Wahrschein- 
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lichkeit 1/t, bzw. K,, der inneren Löschung der gelösten Moleküle beeinfluß wird 
und sich dies auf die Größe von R auswirkt. Die innere Löschung von An- 
thrazen in flüssiger Lösung wäre demnach verhältnismäßig gering und nähme 


-}um den Faktor 100 zu, wenn Anthrazen in Polystyrol einpolymerisiert wird 


(Tab. 6). Bei Fluoranthen wäre das Gegenteil der Fall, nämlich in flüssiger 
Lösung etwa 5mal stärkere innere Löschung als in fester Lösung. Um hier- 
über exaktere Aussagen machen zu können, müssen die Untersuchungen auf 
noch mehr Substanzen ausgedehnt werden. 


Anthrazen und Naphthalin nehmen in Polystyrol gelöst insofern eine 
Sonderstellung ein, als bei den Lösungen dieser beiden Substanzen in Polysty- 
rol Q und R relativ groß sind. Von Anthrazen ist aus spektralen Unter- 
suchungen von Hinrichs!?) bekannt, daß es beim Polymerisieren von Styrol- 
lösungen in die Polysty rolketten eingebaut wird. Die Vermutung liegt nahe, 
daß bei Naphthalin ein ähnlicher Einbau erfolgt und sich dies auf “u Größe 
der Parameter Q und R auswirkt. 


Zusammenfassung der Ergebnisse 


Die drei in den letzten Jahren vorgeschlagenen Mechanismen für eine 
Energieübertragung in lumineszierenden Lösungen unterscheiden sich durch 
ihre Abhängigkeit vom Aggregatzustand der Lösungen. Auf der Grundlage 
der von Kallmann und Furst?) und von Birks®) abgeleiteten Funktionen 
für die Konzentrationsabhängigkeit der Lumineszenzintensität wurde die 
Energieübertragung in flüssigen Lösungen mit der in festen Polystyrollö- 
sungen verglichen. Hierzu war es erforderlich, die Funktion (1) so zu verall- 
gemeinern, daß sie auch für solche Lösungen gilt, deren Lösungsmittellu- 
mineszenz nicht zu vernachlässigen ist (Funktion 3). Die experimentell an 
Polystyrollösungen gefundene Konzentrationsabhängigkeit entspricht inner- 
halb eines großen Konzentrationsbereiches genau der Funktion (3), wobei 
sich in Übereinstimmung mit Messungen von Swank und Buck’) und von 
Krenz!?2) ergibt, daß die von reinem Polystyrol emittierte Lumineszenz etwa 
20% zur Ausbeute der Energieübertragung auf die gelösten Moleküle beiträgt. 
Die Anregungsart, ß-Strahlung oder Röntgenstrahlung, hat keinen wesent- 
lichen Einfluß auf die Energieübertragung in Polystyrollösungen. Die Kon- 
zentrationsabhängigkeit der Lumineszenzintensität von Polystyrollösungen 
unterscheidet sich von der flüssiger Lösungen vor allem dadurch, daß bei ihr 
der Parameter Q der Funktionen (1) bzw. (3) um etwa eine Zehnerpotenz 
größer ist. Dies läßt sich am besten bei einer Energieübertragung durch 
quantenmechanische Resonanz deuten. Zu einem eindeutigen Beweis, daß 
die Anregungsenergie in lumineszierenden Lösungen durch quantenmechanische 
Resonanz erfolgt, reichen die vorliegenden experimentellen Ergebnisse noch 
nicht aus. Aus diesem Grunde wurden die Untersuchungen auf Polystyrol- 
lösungen mit zwei gelösten lumineszierenden Substanzen ausgedehnt. Über 
die Ergebnisse dieser Untersuchungen soll in einer weiteren Veröffentlichung 
berichtet werden. 


Meiner verehrten Lehrerin, Frau Professor Dr. L. Herforth, bin ich für 
die Erteilung des äußerst interessanten Themas und die Förderung meiner 
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Arbeiten zu groBem Dank verpflichtet. Die Diskussionen mit ihr haben mir oft 
geholfen, auftretende Schwierigkeiten zu überwinden. Feinstri 


Herrn Prof. Dr. C. F. Weiss möchte ich danken, daß er als Direkter 
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des Institutes für angewandte Radioaktivität mir die Gelegenheit gab, die wi te 
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Feinstrukturuntersuchungen mittels 
Röntgenstrahlen haben in der natur- 
wissenschaftlichen Forschung und in 
der Technik besondere Bedeutung 
erlangt. 


FEINSTRUKTURUNTERSUCHUNGS -RONTGENAPPARAT 


„TUR“ M 60 
zur Untersuchung des atomaren Aufbaus von Werkstoffen "= = 
nach den üblichen Feinstrukturuntersuchungsmethoden, wie er. a 
Aufnahmeverfahren nach Debye-Scherrer, Laue, Schiebold, & : 
WeiBenberg u. a. 


BETRIEBSSTATTEN FUR RONTGENTECHNIK UND ELEKTROMEDIZIN 
BERLIN N 4 ReinhardtstraBe 18- Tel. 426671, 426672 
DRESDEN N 23 KleiststraBe 10: Tel. 54869, 53842 
LEIPZIG C 1 Roßstraße 17- Tel. 30200, 20215, 65265 


Besuchen Sie bitte unsere ständigen Ausstellungsräume in Berlin N 4 Reinhardtstraße 9 


VEB TRANSFORMATOREN- UND RÖNTGENWERK DRESDEN 


Prospektmaterial stellen wir auf 
Wunsch gern zur Verfügung. 
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Quantentheorie des Atoms 


Von Prof. Dr. A. RUBINOWICZ, Warschau 
Ubersetzung nach der 2. polnischen Auflage des Werkes ,,Kwantowa Teoria 
Atomu“. VII, 486 Seiten mit 38 Figuren im Text. 1959. gr. 8°. 
DM 29.80, Leinen DM 31.80 


Physikalische Blätter: Das vorliegende Buch über Quantentheorie ist das Ergeb- 
nis langjähriger Vorlesungen und erhält im wesentlichen das, was auch bei uns in den 
Vorlesungen gebracht wird. Es enthält seine eigene Note durch das außerordentliche 
didaktische Geschick des Verfassers. An Hand der historischen Entwicklung wird auf 
alle begrifflichen Schwierigkeiten eingegangen. Dabei wird immer wieder die große Rolle 
des Experiments bei der Entwicklung der Quantentheorie aufgezeigt. Das Buch spricht 
besonders auch den theoreihisch interessierten Experimentalphysiker an. 


Einführung in in die Theorie = Elektronenoptik 


Von Dr. phil. JOHANNES PICHT 


o. Prof. fiir theoretische Physik und Optik, Padagogische Hochschule Potsdam-Sanssouci 
(friiher Technische Hochschule Berlin-Charlottenburg) 


2., erweiterte Auflage 
VIII, 274 Seiten mit 70 Abbildungen im Text. 1957. gr. 8° 
DM 29.20, Leinen DM 31.— 


Experimentelle Technik der Physik: Das Wiedererscheinen der „Einfüh- 
rung in die Theorie der Elektronenoptik* dürfte wohl von einem großen Leserkreis dank- 
bar begrüßt worden sein. Sowohl von dem Fachmann, der es schon aus der ersten Auf- 
lage kannte und der es als Nachschlagewerk benutzen will, als auch von den Studenten, 
die in diesem Buch eine wirkliche Einführung in dieses nicht immer ganz einfache Ge- 
biet erhalten haben. Auch derjenige, der sich in der Technik mit elektronenoptischen 
Problemen zu beschäftigen hat, z.B. beim Bau von Fernsehbildröhren und Ikonoskopen, 
kann auf dieses Buch zurückgreifen, da es sehr klar und übersichtlich geschrieben ist, 
Für über den Rahmen hinausgehende Betrachtungen ist ein bis !956 gehendes Literatur- 
verzeichnis vorhanden. Die sorgfältige, übersichtliche Darstellung sowie die Vielfalt der 
beschriebenen Methoden, mit denen man bei elektronenoptischen Problemen zum Ziel 
kommen kann, macht das Werk zu einem guten Lehrbuch und Nachschlagewerk. 
Wisotzki 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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